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In  book  th<jire  is  given  an  anpli cation  of  normal  fundamental 
functions  and  integral  equations  for  solving  engineering  oroblems. 

Examples,  considered  in  work,  refer  to  problems  of  strength, 
stability  and  vibrations  of  elastic  systems,  however,  the  results 
can  be  used  also  in  other  fields  of  engineering. 
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Frequently  the  solution  of  engineering  problems  reduces  i.o  a  .soln^  on  of 
ordinary  differential  equations  or  their  systems  trith  boundary  con.-i j  *  j  ons  of  a 
general  form.  If  the  corresponding  equation  has  a  high  order  and  variAt] j  co¬ 
efficients,  then  the  problem  i*  found  *0  be  difficult,  since  the  fiactang  of  ii. 
accurate  solution  usually  is  not  successful.  This  refers  even  to  second 
equations,  if  they  do  not  reduce  to  knovm  equations  (for  example,  Bessel  eqi  f 
the  solutions  of  which  have  been  tabulated. 

Difficulties  arise  also  in  thor.e  cases,  when  an  accurate  solution  is  known 
(for  example,  for  differential  equations  with  constant  coefficients),  b.  t  w  >  t.^. 
limits  of  interval  of  changes  of  the  independent  variable,  the  sought  function 
its  derivatives  experience  discontinuities  (for  example,  problem  on  fixurc  oi 
rod  under  action  of  concentrated  forces  and  moinents,  problem  on  di;  ' ’•ibut  ion  of 
temperature  in  rods  with  branches). 

The  most  effective  way  of  solving  in  the  latter  case  is  lfi«  application 
normal  fundamental  functions,  as  was  demonstrated  in  works  of  the  outstanding 
Academician,  mechanics  and  mathematics  special. st,  A.  IJ.  Krylov. 

In  those  cases,  when  the  probie...  reduces  to  differential  equations  with 
variable  coefficients,  it  seems  expedient  to  proceed  io  integral  equations.  Th. 
idea  of  such  a  transition  is  intimately  connecteil  with  application  of  melhixl  of 


■uc9«itlvt  tpproxlMiions  for  the  solution  of  differential  equations,  however,  e 
trensltltn  to  integral  equations  makes  it  possible  to  use  more  general  and  more 
effective  solutions. 

In  this  work  there  is  considered  the  application  of  method  of  normal  fundamental 
functions  (Chapters  1  and  2),  and  also  there  are  investigated  boundary  and  normal 
integral  equations  (Chapters  3  and  U). 

Examples  of  application  refer  to  problems  of  Engineering  Mechanics,  however 
the  fairly  general  discussion  makes  it  possible  to  apply  the  results  also  in  solving 
other  engineering  problems. 
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CHAPTER  1 

NORMAL  PJNDAMEIJTAL  FUNCTIONS  OF  LINEAR  Dl  FFEliENTInL 
EQUATIONS  '.VITH  CONSTANT  COEFFICIENTS 

Effectiveness  of  use  of  normal  fundamental  functions  in  engineering  problems 
was  established  in  widely  kj.owr;  w.  rks  of  A.  N.  Krylov.  In  the  subsequent  works  of 
P.  F.  Papkovich,  Sh.  E.  Mikeladze,  N.  K.  -nitko  and  others,  these  functions  were 
anplied  in  solving  a  number  of  problems  of  structural  mechanics. 

The  special  advantages  of  normal  fundamental  functions  are  reflected  in 
constructing  discontinuous  solutions  of  differential  equations  with  constant 
coefficients  (of  solutions  with  given  discontinuities  of  the  derivatives). 

Comprehensive  experience  in  constructing  and  using  such  solutions  in  problems 
of  structural  mechanics  (beginning  with  the  known  nroblem  of  integrating  equatio;i 
of  elastic  line  of  rod)  made  it  possible  to  generalise  the  results  for  linear 
differential  equations  of  arbitrary  structure  (Works  of  oh.  E.  Mikeladze). 

In  this  work  there  are  established  general  formulas  for  determining  normal 
fundamental  runc..ions  and  simnle  differential  relationships  between  them. 

By  means  of  general  relationships  there  are  obtained  already  well-know-'i  systems 
of  normal  fundamental  functions  and  there  are  given  certain  applications  of  these 
functions.  In  particular,  they  are  used  for  an  approximate  integration  of  aiffer- 
ential  equations  with  variable  coefficients. 

The  solution  of  differential  equation  is  presented  in  matrix  form,  which  is  a 
mathematical  expression  of  a  known  method  in  structural  mechanics  of  initial 


pAr«iMt«rt  (Works  of  A.  A.  Umanskly,  N.  I.  Bezukhov,  N.  G.  Chudnovskiy  and  others). 


1.  Statement  of  Problem 

Suopose  there  is  given  a  linear  differential  equation  of  n-t.h  order 

.  .  +pA<)y(<)  =  f(x)  (1.1) 


or 


iPi(^)y-'‘(-')=/('«f) 

tmd) 


(1.^) 


the  solution  of  equation  (l.l)  in  a  certain  interval  of  change  x{a<x  ''b)j 
is  sought. 

The  totality  of  n  (linearly  independent)  solutions  of  the  homogenec.us  equation 
(1.1/  (A“0,  1),  satisfying  the  condition 

1  (=•  A 

Oirk  (j,  A=0,  I . n-l),  (1-3) 


is  called  the  normal  fundamental  system  of  solutions  of  equation  (l.l)  vdth  the 
initial  section  x  =  a. 

The  particular  solution  of  equation  (l.l),  corresnonding  to  zero  initial 
conditions,  is  designated  y»(Jf)- 
Thus , 


K(0(a)  =  o  (i-0,  1.  .  .  .  ,n- 1).  (1.4) 

(If  (K*(JC)}  and  K.  (x).  ,  are  known  then  the  solution  of  equation  (l.l) 

is  presented  as: 

y  M-Jy"H<>)  y,<’c)  +  y,  (a), 

•  -•0 


(1.5) 


where  y***  (o)  (*=0 . n — 1)  —  are  values  of  the  function  y(x)  and  its 

first  n — 1  derivatives  in  the  initial  section  x  -  a.  The  indicated  values  arc 
called  also  initial  parajneters.  In  solving  bo'indary  value  problems,  there  is  use 
usually  not  only  the  function  y(x),  but  also  its  derivative  un  to  n — 1  order 
inclusively. 

For  the  future  it  is  exnedient  to  introduce  "the  coluiriTi — solution" 


b(-0)  = 


y(^) 

y<n(x) 


(1.6) 


(x) 


From  equality  (1.5)  by  successive  differentiation  we  find 

l>(^)l  =  lK(x)|l>-(a)|+(F,(x)]. 


where 

"K,(x) 

(  V) 

.  .  y-iix)  - 

Kl'>  (X) 

■ 

F{— )(x)  . 

(3.7) 


(l.B) 


is  a  normal  fandamental  matrix  of  the  homogeneous  equation  (1.1),  and 


y(a) 

-  K.(x)“  - 

yO)(a) 

.  (K,(x)l  = 

y«-ii  (a) 

is  the  column  of  initial  parameters  and  the  column  of  a  particular  solution.  The 
normal  fundamental  matrix  in  initial  section  is  unitary  and  column. —  the  part¬ 
icular  solution  —  is  zero.  The  solution  (1.7)  corresponds  to  the  application  of 


the  method  of  initial  parameters,  vd.de ly  us«d  in  engineering  problems.  This 
solution  shall  call  the  solution  in  matrix  form. 


2.  The  Homogeneous  Hkjuation 

Let  us  consider  a  homogeneous  differential  equation  vdth  constant  coefficients 


{...!) 


.juopose  is  characteristic  nolynomial  of  equation  (2.1): 


1-0 


(2.2) 


the  roots  of  vrhich  v<e  shall  designate  X,  (s  =  0.  1....). 

If  vrt  subordinate  the  selection  of  arbitrary  constants  to  the  condition 

y<'>(a)  =  V  (^-0.  j 


v/here  is  a  certain  parameter,  then  the  solution  of  equation  (2.1)  vdll  be 

* 

such 


y{x)^-B  < '  -o'. 


v/here  B  is  the 
If  nov/  we 
ij,  ,  then,  as 
equation  (2.1) 


symbol  of  operation  of  a  cmnlete  integral  residue, 
exnand  expression  (2.4)  into  a  series  1/  degrees  of  the  rxarameter 
was  established  even  by  Cauchy,  the  riormal  fanda’nent.al  iunctio.-.s  of 
are  found  to  be  the  coefficients  of  the  exnansior; 


»'|(JC)  V  d-  ■  •  +^-i(.0V  '• 


A.  N.  Krylov,  On  Certain  Differential  hlquations  oi  Mathematical  Physics,  CITTL, 

M.— L.,  1950. 


After  discussing  this  result  A.  N.  Krylov  in  the  work  'On  Certain  Differential 
Equations  of  Mathematical  Physics"  turns  to  the  consideration  of  concrete  differer. 
equations  of  simole  structure,  for  which  he  also  makes  the  indicated  expansion. 
However,  it  is  possible  to  establish  certain  general  results,  valid  for  differenti 
equations  of  arbitrary  order  with  constant  coefficients. 

Vh  shall  assume  at  first  that  the  roots  of  the  characteristic  pclynomial 

X»,  Xi . X,-i  are  all  simple. 

In  introducing  the  vaiue 


A  IJ 

i*  0 


(i-l  K-l  —  k 

»-0  i-0 


into  equality  (2.4),  by  making  a  calculation  of  integral  residue  and  by  ccroiderif. 
the  exnansion  (2.5)>  we  obtain 

"-I  V 

^  - 

,..o  (2.7) 


(*  =  0,  1.  .  .  .  .  n-1; 


In  this  equality  X,  (s  =  0.  !,.  .,n  1)  are  roots  of  the  characteristic  nolynom. 

In  a  narticular  case  for  function  )',_i(.v)  there  is  obtained  the  following 
exoression; 


I 

I 


f  ^ 


V  — Oa, 
1^0 


1 


(2.8) 


*Thi3  result  also  can  be  obtained  by  methods  of  onerational  calculus. 


Wt  now  considtr  the  caae  of  multiple  roots.  Suopose  the  characteristic 


polynomial  (2.2)  has  m  different  roota  X*  («=0. 1,,  .,/n— 1)  with  a 

multiplicity 

Relationships  (2.4)  —  (2.6)  remain  in  force  also  for  the  considered  case  and 
therefore 

d-i-* 

- (2-9) 

*  ’  F{k) 


After  calculating  the  complete  integral  residue, we  will  find 


-I 


1-0 


JH-l 

— ^ —  n 

(1  _  i^y.  11' 


gX{x-a) 


U.io) 


Differentiation  in  this  equality  is  conducted  with  respect  to  X  and  into  final 
result  there  is  introduced  X  =  Xi. 

If  all  V,  1  (s  0,  ...,m-l),  then  m  -  n  and 

a-t  I 

»  —  f-i  I 

1-0  U-X, 


»-l 

■  P,'(n  ~  i) . 

I'.O 


formulas  (2.7)  and  (2.10)  coincide. 

At  k  n-1  from  equality  (2.10)  we  find 


y'.-.(^) 


- 1)'  dr--‘ 


( » •  <i> 

m  -  I 


(2.11) 


In  the  future  there  will  be  given  examnled  of  use  of  formulas  (2.7)  and  (2.10). 

We  note  that  in  the  formula  for  y^ix)  {k  =  0,...,n — 1)  the  value  ol 

the  coefficient  P»-  does  not  enter.  Its  magnitude  exerts  an  influence  only  on 
the  value  of  the  roots  of  the  characteristic  oolynomial. 


3.  Recurrent  Relationships  Between  Normal 
Fundamental  Functions 


From  formulas  (2.7)  and  (2.10)  it  is  oossible  to  establish  the  following  La:.r 
relationshin: 


Ak»(x)  =  K»-,  (A)-;>,_^K„_,(.r) 

dx 


(*  =  0.  1 .  «-l). 


(3.1) 


In  this  formula  it  is  necessary  to  ass’ome  H(.v)^0  at 
examnle  for  derivative  yo(^)  we  shall  have 


.  Tnus,  for 


dx 


K,(JC)=  -p„Yn-x(x). 


Equality  (3.1)  makes  it  oossible  to  seek  a  system  of  normal  fundamental 
functions  of  an  equation  with  constant  coefficients,  if  there  is  known,  for 
example,  function  y»-i(A)  .  "his  method  frequently/  is  found  to  be  in  nrac-.i.a. 

problems  very  effective,  since  the  determination  of  y.-t(’<)-  by  lormulas  (2.3)  an 
(2.11)  are  relatively  simple. 

On  the  basis  of  equality  (3*1) 

K,_,(x)  K,-.  (a)  (  ^). 

dx 

y’«-3(A)  — ^  K/I-J  (a)  +/’li  n-l  (a), 


what  results  in  a  subsequent  determination 


of  all  yil  v)  (fe=0,  1). 


Th«  recurrent  relationship  (3.1)  makes  it.  possible  to  express  derivatives  of 
the  f’inction  (i,  A“0,...,n— I)  by  a  linear  combination  of  the  normal  .''undamental 

functions. 


4.  Solution  of  Inhomoj^eneoua  Fxjuation 
Suppose  there  is  given  the  nonhomo^:eneou3  equation 


I  ^0 


(4.1) 


The  oarticular  solution  of  equation  (4.1)  satisfying  the  zero  initial  conditions 
may  be,  as  knovm,  presented  in  the  following  fonn; 

y^{x)  =  ]  + 


which  readily  is  verified  by  a  direct  substitution.  This  resul:,  can  be  estHtliohe  i 
also  by  means  of  theory  of  integral  residue,  if  equalities  (2.8)  and  (2.11)  are 
used . 

For  an  explanation  of  the  writing  in  the  form  (4.2)  we  shall  present  an 
illustrative  example. 

For  the  equation  y^^Hx) y  (x) /  {x) 


we  have 

K,  (t)=-.  cos(x--a), 

Kj  (x)  -  sin  (a  — 

J'.  (-<)*=  f  J'-t-i  ('*^  —  5  d  «)/(■')  ]  sln(A-  s)/{s)ds. 

€  a 

General  solution  of  equation  (4.1)  will  be  thus: 

j,(x)-=*Vy(*>  (u)  y'*(A)  i  )  )'„  ^  ,  .j)/(s)(/s.  (/,.3) 

*-o  “ 

The  solution-  in  matrix  form  has  the  form  (1.7). 


Jo 


5 .  Discontinuous  Solutions 


VJe  shall  seek  the  solution  of  equation  (4.1),  satisfying  the  given  initial 
conditions  and  having  given  discontinuities  (discontinuities  of  fi^st  order)  of 
derivatives  un  to  n — 1  order  inclusively. 

Such  a  type  of  oroblem  is  encountered  during  calculation  of  concentrated 
influences.  We  note  that  coefficients  of  equation  do  not  have  discontinuities 

’/C 

in  entire  interval  of  change  x  . 

Sunnose  the  function  >'‘’'’(-f)  ('  =  0 . n  — 1)  has  m  discontinuities, 

located  in  the  section  x  —  n,/ (1 'Jy -4^ '«.)•  The  discontinuity  of 

^(•»(jc)  in  section  x  =  n./  ^  will  re  designated  as: 

y ’>(a, ,  +  0)  -y<'>(a.  y  -  0)  =  ( 5  _  1 , 

(v=C  1,  .  .  .  ,  n—  1;  y  =  1.  .  .  «.), 

Function  nossessing  given  discontinuities  must  have  the  following 

structure : 

y)(jc)  =  ;p,(x)  +  Jv5(.v,  a.j)  (5.2) 

(v  =  0.  1 . rt-I). 

where  S{x,  a,j)  is  a  single  function,  determinated  by  the  equality 


S(x,  a.y)=  jj 


X  <:  n,y, 
x>a.j. 


(3.3) 


c. (jc)  is  a  continuous  function. 

Function  y(x)  which  is  a  solution  of  the  stated  nroblera  has  to  satisfy  equation 
(4.1),  the  initial  conditions  and  condition  of  discontinuities  (5.2). 


*  Presence  of  discontinuities  of  function  y(x)  and  its  first  n — 1  derivatives  is 
not  associated  generality  speaking,  with  the  continuity  of  coefficients  of  the 
differential  equation.  There  may  take  place  also  the  reverse  case,  i.e.,  continuous 
solution  in  the  presence  of  discontinuous  coefficients. 


It  e*n  b«  estAbllshad  that 


y  (JC)  *[2  (“)  t"')  +% ^  ('')+!".(*)  ( 5 ,  ^ ) 


corresDonds  to  the  indicated  conditions,  if  functions  of  y^/(x)  satisfy  the 
hoeiogeneous  differential  equation  (4*1)  and  the  relationship 

(1  v-ik. 


0  vv>A 

(»,  *  =  0.  1 . «— 1). 


(5.5) 


Functions  of  Yi^ix)  should  bo  continuous  together  vdth  the  n-1  derivatives. 
From  the  preceding  it  is  clear  that  function  K»y(A)  is  normal  fundamental 

function  with  the  initial  section  x  =  a^j-  It  is  sufficient  to  assume 

>'*/(*)  =  ^(^-«.y).  (5-6) 


so  that  all  the  Indicated  conditions  above  were  satisfied. 

Thus, 

>  (  0  =  s'y “ (fl)  (^)  +  "2  2  -S (X,  a,j)  \y^y, {X  - a,j)  + 


+  — s  +  a)/(s)ds. 


(5.7) 


The  solution  (5.7)  can  be  written  in  a  more  symmetric  form,  if  initial  values 
of  the  function  are  considered  as  the  given  discontinuities  after  assuming 

>r(a)  =  AW, .  .  .  ,  y<''-'>(<i)  =  Ai— ■> 

and  by  assusdng  all  (^  =  0, .  .  n  — I). 

Now,  equality  (5.4)  we  shall  write  out  as: 

>(x)="s  "S5(x,  a*,);i<*>K*(.v-a,y)  +  K,(.t).  (5.8) 

The  solution  of  type  (5.8)  long  ago  was  used  in  structural  mechanics  for 

equations  of 'particular  forms.  By  another  method  and  in  another  form  formula  (5.8) 

* 

was  established  by  Sh.  E.  Nlkeladze  .  However,  in  the  reasoning  of  the  author  an 


■*Sh.  E.  Mlkeladse,  Certain  Problems  of  Structural  Mechanics,  State  Engineering 
Publishing  House,  Noecow,  1948. 


IZ 


error  creot  in:  equality  (5.8)  is  correct  only  for  an  equation  with  constant 
coefficients,  since  only  in  this  case,  functions  ai.,)  satisfy  the 

corresDondinp  differential  equation. 

Let  us  turn  to  the  solution  in  matrix  form.  For  convenience  in  writing  we 
assume  that  in  section  x  =  aj there  is  in  a  general  case  a  discontinuity  of  all 
derivatives  which  characterizes  the  'column  -  discontinuity". 


yu. 


(5.9) 


Certain  elements  of  this  column,  of  course,  may  be  equal  to  zero. 

The  discontinuous  solution  of  equation  (4.1)  in  matrix  form  will  be: 

I>'(.»^)1  =  ^£5(a,  ay)lK(.v-a^)]lAj+(K,(x)),  (5.10) 

where  m  is  the  number  of  sections,  in  which  there  are  discontinuities. 

■jhe  initial  column-discontinuity  is 


IA,1  = 


>(<*)  I 

3»(')(a) 


y<"-'*(a) 


(5.11) 


The  matrix  [y(x — o/)]  has  the  form  (1.8>,  where  functions  >i(Jc)  are 

reolaced  by  y»(x— «>)•  If  all  elements  /  J  are  given,  then  the 

discontinuity  will  be  called  ii.dependent.  Frequently,  however,  there  are  en¬ 
countered  oroblems,  in  which  discontinuity  of  der  vative  y'*’ (Jt)  in  the  section 

derivatives/ 

x—Oi  deoends  on  values  of  function  y(x)  and  its  /  in  the  same  section 


(in  th«  presence  of  discontinuities  wa  shell  for  definiteness  consider  the  left-hand 
^•luee  of  functions  in  section  the  results  alaost  without  change  are 

applicable  for  the  case,  when  one  should  consider  right  side  values). 

Thus, 

(4yi-[CJI>(a^)l. 


(5.12) 


In  this  case  the  discontinuity  will  be  called  dependent,  and  i\Cj  is  the  matrix  of 
deoendent  discontinuity.  In  theory  of  rods  frequently  there  is  used  the  equation  of 
fourth  order,  in  which 


Aj»> 

Thm 


{Cj\ 


^tojy  i^j) 

C^jy  («>)  + 

~  0  0  0  0“ 

0  0  0  0 

_  ^uj  0  0  _ 


(5.13) 


(5.14) 


Discontinuous  solutions  can  be  anplied  both  in  case  of  concentrated  influences  on  a 
system,  and  in  the  presence  of  discontinuity- like  variation  of  the  parameters  in 
different  sectors  of  the  system,  being  described  by  one  and  the  same  differential 
equation  (for  example,  flexure  of  stepped  rod). 

6.  Rxamnles  of  Applying  Normal  Fundamental  Functions 
As  an  example  ws  shall  consider  the  problem  on  flexure  of  a  rod  of  constant 
section  under  action  of  given  external  forces  (Fig.  1). 

The  differential  equation  has  the  form 


y«>(x) 


El  ' 


(6.1) 


where  y(x)  is  the  sag  of  axis  of  rod;  EJ— the  strength  of  the  section  to  flexure, 
f(x)— the  distributed  load  per  unit  of  length  of  beam. 


Characteristic  polynomial  of  the  homogeneous  equation 

has  root  X»=0  of  fourth  multiplicity  ■'0  =  4. 


y 

rd) 

tlfUff 

» f  1 1 1  f  hfiT » » li  1 1 } 

44 

^  ■■  ■' 

- -1 

- - 

Fig.  1.  Flexure  of  Pod. 

From  formula  (2.10)  we  find  (n-4,  m-1,  a-0) 


3' 


From  equalities  (3.2)  we  determine 


K.(x)  =  ^K,(.v)  =  |f. 
K.(x)  =  ^^K.(x)  =  -^. 
K.(x)  =  £-K.(x)  =  l. 


Of  course,  this  system  could  have  been  written  by  not  resorting  to  total  results. 
The  function 


31 


At  points  of  applying  concentrated  bending  moments  M,  and  forces  P<  there  are 
discontinuities  of  the  derivatives 


J  FJ 


EJ 

Pl 

EJ 


lb 


On  the  basis  of  aqiiality  (5«7)  we  obtain  the  known  equation  of  the  elastic 
line  of  a  rod: 


jf  (X)  (0)  +y<‘>  (0)  X  (0)  ^  +>‘3) («)  ^  q. 

S  (x,  a^)  Aj®>  +  ^  5  (x,  o,y)  Aj' >  (x  —  a,y)  + 

/-I  i-l 

0 


(6.2) 

In  composing  the  eqxiation  it  was  assumed  that  the  elastic  line  has  in  the  sections 
0*1  discontinuities  of  sags,  and  in  the  sections  Oi;  — discontinuities  of  angles 
of  rotation. 

In  practice,  such  a  case  can  be  encountered  for  a  coopound  rod. 

Equation  (6.2)  will  be  valid  also  for  a  rod  with  a  graduated  change  of  the 
section,  if  one  were  to  introduce  in  corresponding  sections  the  discontinuities 


AJ»> »  M  (bj)  £y  (6 ,  _  0)  ]  • 

■=  Q  (bj)  J  • 

whex*e  M(bj)  •nd  bending  moment  and  the  transverse  force 

in  the  section 

The  normal  fundamental  matrix  of  equation  (6.1)  will  be  thus: 


1 

X 

21 

0 

1 

X 

0 

0 

1 

_0 

0 

0 

JG 


(K(x)l  = 


(6.3) 


Wb  consider  now  flexure  vibrations  of  a  weightless  rod  of  constant  section, 
carrying  masses  m)  and  moments  of  inertia  (Fig.  2).  In  these  sections  the-f> 

will  be  the  dependent  discontinuities. 


where  •-  is  the  angular  frequency  of  natural  oscillations. 

The  column  of  the  dependent  discontinuity  is  e:.nressed  as: 


1 

1 

1  / 

AJ*> 

0 
0 

0  - 

■’m  f 

El 


0 

0 

El 

0 


0  0 
0  0 

0  0 
0  0 


y(oj) 

y<‘>(ay) 

y(3)(a  )  J 


(6.5) 


In  sections,  wnere  supports  are  located  there  will  be  a  discontinuity  of  the  seconu 
and  third  derivative  in  accordance  with  magnitudes  of  reactive  moment  and  force . 

In  majority  of  real  cases  Ar*,/=0  and  K»t,  =  0;  for  elastic  supports /Cs*/<0.  K>t, 


Fig.  2.  Critical  Sneed  of  Shaft 


The  matrix  of  discontinuity  for  an  elastic  support  has  the  form 


0  0 
0  0 


^tlj 


0  0" 

0  0 
0  0  I' 
0  0 


A  rigid  support  may  bs  taksn  into  account  by  selecting  the  corresponding  cocTficients 
of  rigidity  of  the  elastic  support.  We  note  also  that  a  rigid  support  which 
eliminates  ths  sectior  of  linear  mobility  is  equivalent  to  an  application  of  m<-* **oc;a 
support  eliminating  angular  turns  corresponds  to  an  application  The 

solution  of  the  problem  in  matrix  form  is  exnressed  by  the  equality 


where  /”y(x)_7  i*  the  fundamental  matrix  of  equation  of  flexure  of  rod;  [C/]  is 
matrix  of  discontinuity,  corresponding  to  concentrated  mass  or  elastic  support. 

Equation  (6.6)  is  useful  also  for  calculating  for  critical  speed  of  a  weightless 
shaft,  loaded  with  separate  disks. 

■» 

For  the  case  of  forward  synchronous  precession  the  matrix  of  discontinuity  in 
section,  where  disk  is  located,  will  be  thus: 


“  0 
0 

0 

_  BJ 


0 

0 


0  0~ 

0  0 

0  0  . 


0  0  0 


where  iHi-  is  the  mass  of  disk;  !$■  —is  the  equatorial  moment  of  inertia.  Thus, 
for  «xaaq>le,  the  discontinuity  of  second  derivative 


*The  concept  of  forward  and  reverse  synchronous  processions  is  given  in  Chapter 
4,  potion  3. 

**lW.s  equality  in  work  of  A.  N.  Krylov,  "On  Detenaining  the  Critical  Speeds  of  a 
Revolving  Shaft"  (Academy  of  Sciences  of  USSR,  Moscow,  1931)  and  later  also  in  the 
work  of  Sh.  E.  Milwladse  "New  Methods  of  Integrating  Differential  Equations"  Moscow, 
State  Theoretical  Technical  Publ.  House  1951,  is  erroneously  used  with  the  minus  sign. 
It  follows  from  this,  that  the  calculation  of  gyroscopic  moment  of  disk  during  a 
forward  synchronous  precession  does  not  increase,  but  lowers  critical  speed  of  shaft; 
this  is  not  true. 


In  calculating  a  reverse  synchronous  precession  for  thin  disks  one  should  consider 

The  calculation  by  equation  (6.6)  must  be  made  in  sequential  order,  in 
determining 

^  1|  2«  3f  .  .  y  fft)m 

Suopose  there  is  a  certain  number  of  unknovm  initial  narameters  (in  considered 
problem  there  are  two);  then  the  same  number  of  homogeneous  boundary  conditions 
should  exist  at  x  =  b.  As  a  result  of  the  calculation  vre  obtain 

=  (6.8) 

In  this  equality,  the  column  [y(fl)l  contains  only  the  unknown  initial  parameters, 

and  the  column  [y(&)l  —  those  values  of  y*"*  (b)  ('*-0 n — I),  which  enter  inu 

boundary  conditions  at  x  =  b. 

In  considered  case  boundary  conditions  are  such: 

(y(fc)]=o, 

and  then  from  relationships  (6.8)  tnere  ensues 


delM«]=0. 


(6.9) 


Since  the  coefficients  contain  »*.  then  equality  (6.9)  is  the  characteristi 

equation  for  determining  the  eigen  values. 

Problem  of  calculating  consists,  essentially,  in  deteriming  of  elements  of  the 
matrix 

We  discuss  now  in  reference  to  the  considered  scheme  of  calculation  one  method, 
belonging  A.  N.  Krylov  (See  "On  Determining  the  Critical  Sneeds  of  Revolving  Shaft"), 
this  method  makes  it  possible  in  most  cases  to  simplify  the  calculation  considerably. 

Suppose  the  element,  standing  in  the  i-th  line  of  column  of  initial  parameters 
is  unknown  (for  definiteness,  we  assume  ^=2). 


Than  wa  introduce  the  unit  column 


!>(«)•) 


I  0"! 


(6.10) 


where  all  elements  — zeros,  excent  standing  in  line  i.  With  the  initial  column  of 
(6.10)  wa  conduct  antire  calculation,  which  determines  elements  of  matrix  [>!(»], 
standing  in  column  i.  After  having  made  so  many  such  calculations  of  so  manj'  un¬ 
known  initial  parameters,  we  determine  all  the  elements  of  matrix. 

We  note  that  in  a  majority  cases  vd.th  use  of  equalities  (6.7)  there  is  no 
necessity  to  calculate  all  elements  of  matrices  — Oi)].  which  also 

facilitates  the  conduct  of  calculation. 

In  calculating  the  mass  proper  of  rod  (or  shaft)  wo  oroceed  frcan  the  eqiiation 


><«>  (x) — •**«>  (x) = 0, 

tdiera  y(x)  —  amplitude  sag  of  axis  of  rod; 


(here  fF  —is  mass  of  unit  of  length  of  shaft). 

Characteristic  nolynamial  of  equation  (6.11) 

f  =  — oi*a‘ 

has  the  roots 


(6.11) 


)l,b  —  X.,=  — iaV^to. 

On  the  basis  of  forjaula  (2.7)  we  obtain  /unctions,  introduced  by  A.  N.  Krylov 


•  3 


=*  Yi»x  +  cos  a  V^x); 

a  ^ 

a  V^u>x  +  sin  a  «x); 


#-o 


a  y<x>x  —  COS  a  y  u)x); 


K,  (X)  -  2  T  • 


(6.12) 


The  functions 


yo(x),  K,(a).  K*(a)  also  can  be  obtained  by  means  of 

recursion  formulas  (3.1).  In  using  relationshio  (3.1),  we  will  find  fundamental 
matrja  of  equation  (6.11) 


~  ^(A) 

>",(A) 

^(A) 

>',(A)- 

I»'(A)1  = 

®*»*^'s(a) 

Ko(a) 

(a) 

^'a(A) 

m"a‘K,(x) 

^(A) 

(A) 

_»VK,(x) 

<u*a‘Kj  (x) 

cu*o«K,(.r) 

^(A) 

All  preceding  results  relative  to  vibrations  of  red  remain  in  force,  if  o:.ly 

* 

instead  of  matrix  (6.3)  we  use  matrix  (6.13)  .  To  the  same  degree  this  refers  also 
to  calculating  the  critical  numbers  of  revolution. 

Vie  nresent  an  examnle  ,  referring  to  determination  of  critical  speeds  of  shaf: 
with  one  disk  (Fig.  3)  in  calculating  the  mass  proper  of  shaft. 


*If  the  mass  of  rod  becomes  vanishingly  small  (a-*0).  then  matrix  (6.13) 
transforms  into  matrix  (6.3). 

This  example  is  found  in  work  of  A.  tV.  Krylov  "On  Determining  the  Critical  .3need 
of  Revolving  Shaft",  in  Sect.  10.  However,  in  tne  solution  errors  crept  in  and  made 
the  result  incorrect.  Besides,  the  already  mentioned  error  with  the  sign  of  the 
gyroscopic  moment,  in  equation  (58)  of  this  work  and  subsequent  equations,  the  ex¬ 
pression  for  y<'>  is  written  without  consideration  of  influence  of  discontinuity. 


The  coluBn  of  initial  values  has  form 


U.l  =  (C.I(3^(a)H 


The  column  of  first  discontinuity  (x  =  ) 


■  C 
yC)  (0) 
0 

y(3)  (0) 


0 

0 

0 

0“ 

~  yi^i) " 

t 

0 

0 

0 

0 

IA.l  =  IC,l(y(a,)l  = 

0 

Cll*/f 

EJ 

0 

0 

y>(a,) 

(a,) 

_  EJ 

0 

0 

0 

_y(3)  (a,) 

.L_ 


tj 


0 

0 

><■>(«,) 

y(ox) 


(6.14) 


Fig.  3-  Shaft  with  one  disk. 


_  EJ 

Since  boundary  conditions  at  x  =  1 


>(0“0.  y<»>(/)=0. 


then  there  must  be  calculated  in  matrixes  of  equation  (6.7)  only  these  two  lines, 


and  for  determining  the  discontinuities,  also  y(Oi)  and  y ’  (**1) 

have  from  relationships  (6.7)  and  (6.13) 


X 

X 

“  0  “ 

X 

^ i(®i) 

X 

>'.(«!) 

yC)  (0) 

y»(ai) 

X 

X 

X 

X 

0 

_  y»(a.)  _ 

_x 

X 

X 

X  _ 

1 

_ ! 

We  shall 


(6.15) 


Here,  elements  of  matrix,  not  participating  in  calculations  are  marked  x  . 


From  (6; 15)  it  is  evident 

Jf  (a,)  -y(«>  (0)  K,  (a,)  +  y*3>  (0)  K,  (u,). 
y«»  (a,)-y)  (0)  K,(a,)  +y)  (O)  K,(a,). 


(6.16) 


For  the  section 


a,  =  / 


-\y{0\- 

X 

(0 

X 

»',(0’ 

"  0  “ 

y»  (/) 

X 

X 

X 

X 

yi)(0) 

4. 

y»>  (/) 

X 

a>»»‘K,(/) 

X 

0 

1 

y»)  (/) 

_x 

X 

X 

X  _ 

y3)  (0) 

X 

X 

' 

0 

0 

X 

X 

X 

X 

X 

X 

YAl-a,) 

U>*/| 

-y»(a,) 

_x 

X 

X 

X  _ 

EJ 

y(ai) 

Hence 


j;(/)=y'»(0)K,(/)+y3,  (0)  K,(/-a.)X 

Xy>(a.)  +  -^3J->'3(/-a.)>'K). 
yo (/)=><“ (0) ..VK, (/)  +y^)  (0)  K.  (/)  4-  “^5- X 

X  K,(/-a,)y»(a,)  +  K,  (a,)- 


(6.1?) 


In  introducing  values  of  (6.16)  and  consxdering  condition  (6.14), 
a  system  of  homogeneous  equations 

y»  (0)  [k,  (/)  -f  Y,  (/-a.)  Y,  (fl.)  + 

+  -^  »'.(/-  a.)  K.  (a,)  j  +y>3)  (Q)  jr,  (/)  +  X 

X  K,  (/  -  a,)  K,  (a,)  +  Y,(i-a,)Y,  (a.)]  =  0. 

y*) (0)  {»VK,(/)  +-^5^  n (/-a,)  Y,  (a,)  + 

+  K.  (/-  fl.)  r,  (a.)J  +>.(3)  (0)  ( K,  (/)  +  ^  X 

X  K,  (/- a.)  K,  (a,)  +  K,  (/ -  a,)  K,  (a,) j  =  0. 


we  obtain 


(6.18) 


j^ll  ^tt 

y»(0) 

L^ii  ^is 

y3)(0) 

or  in  matrix  writing 


Th«  •quality  to  s«ro  det  Mai 


gives  the  characteristic  equation  for  detennining 


The  roots  of  equation  are  most  simply  found  graohically  by  means  of,  constructing 
the  function 


F  (••)  *■  ^ii^n  —  ^  1*^11 


(6.19) 


and  in  aoproaching  limits  by  means  of  a  linear  interpolation  between  the  ooints 
F>0  and  ^<0..  We  note  that  in  matrixes  there  are  calculated  elements, 
standing  in  columns,  corresponding  to  lines  in  columns  of  the  discontinuities.  For 
intermediate  sections  there  are  maintained  only  lines  necessary  for  calculating 
the  discontinuities;  for  the  last  section  —  lines  corresponding  to  boundary 
conditions. 


With  a  large  number  of  discontinuities  the  presented  method  becomes  cumbersome, 
since  all  the  time  it  is  necessary  to  determine  the  values  [>'(«*)]  by  the 
Initial  parasmters. 

In  this  case  it  is  expedient  to  use  method  of  A.  N.  Ki^lov,  the  application  of 
which,  MS  shall  illustrate  in  this  same  example.  We  introduce  the  column  of  initial 

values 


and  then  from  equalities  (6.15)  we  obtain 

>(«!)•“  »'i(ai). 

y'»(a,)*-K,(a,). 


Then  we  shall  have 


|7^ 

y)(0* 

y  *>(/)* 


X 

»',(/)  X 

X 

-  O' 

X 

X 

> 

\ 

X 

-VK^/)  X 

X 

0 

_  X 

X 

X 

X  _ 

0 

X 

X 

K,(/-a.) 

^,(/-a.)' 

X 

X 

X 

X 

u.*/, 

X 

X 

»',(/- a.) 

V',(/-u.) 

X 

X 

X 

X  _ 

~Tj 

+ 


0 

0 


hence 


y(/)*  =  ^.i“  YAI)  +  ^;  >',(/-fl.)Ko(a,)  + 
y«(e)*  =  =  <uVK, (/)  +  -^y  Ko  (/-  a.)  K, (a,)  + 


Analogously  we  find  and  and  then  function  (6.19). 

In  taking  into  consideration  the  mass  and  moment  of  inertia  of  sections  of 
rod  or  during  action  of  constant,  longitudinal  force,  the  differential  equation  of 
the  oroblem  will  have  the  form 

y«>(jc)+/>,y’>(x)+/>«v(.v)  =  o. 

(6.20) 

We  shall  determine  the  normal  fundamental  functions  of  this  equation,  by  using 
results  of  Sec.  2  and  Sec.  3- 

The  rootr  of  characteristic  nolynomial 


M»  dvsignAt* 


I 

V. 

% 

Further  on*  should  consldsr  ths  rslatlonship 


F***-F«- 


Frosi  squsUty  (2.8)  ws  obtain 


In  using,  now  rslstionship  (3.1)«  we  find 


Y  ^  A 

*'  '  ^  V^W  V-t-W 

,  »  [»  ^JC^±  vxl. 

!»  4i>*+2g*'»  ^  I  S  flM  ’  #<•«' 


nw-^i'.w- 
K.  w-j  n  w+i^f.  w 
nw-^nw- 


- - ! - I  f|»JC—  '•Jcj, 

p»— »*  Aan  c»«k 

_ L-(^  vxl. 

ii*— I*  Sen  '*  ws<i  ’ 

I**— ’i*  <e*S 


SqMslitF  (3«1)  adM  it  fMtlbla  t«  d«i«niins  all  alsnents  of  fuadsnantsl  matrix  of 

equation  (6.20). 

7.  Other  Applications  of  Normal  FupriyMnt.al  Functions 
Let  us  consider  an  STimatrle  deformation  of  a  closed  cylinderic  shell 

(Fig.  4).  The  radial  displacement  of  points  n5ddle  smrface  we  designate  as  y(x). 


1''^'^.  4.  iiidally  r^'^iinetric  defor::iation 

of  a  cylindrical  shell. 


Z6 


The  differential  eqiiation  for  the  function  y(x)  has  the  fonr. 


In  this  equality 


where  ^  I2(|  — |t*)  is  the  cylindrical  rigidity; 

is  the  coefficient  r**  linear  expansion; 

temperature  of  middle  surface  of  shell; 

A<(X)-  drop  in  temperatures  between  outer  and  inner  surfaces 
shell. 

Equation  (7.1)  conforms  with  the  equation  of  flexure  of  beam  on  an  elastic 

* 

foundation,  for  which  for  the  first  time  there  were  widely  used  by  A.  N.  Krylov  , 
normal  fundamental  functions.  We  shall  obtain  an  expression  for  these  functions, 
using  the  general  formula  (2.7). 

The  characteristic  polynomial 

=  + 

has  the  roots 


^-PCI+O;  ^=P(i-');  ^=-P(i-0;  ^--Pd+O- 


A.  N.  Krylov,  On  the  Calcuation  of  Beams  Lying  on  an  Elastic  Foundation. 
Academy  of  Sciences  of  USSR,  Moscow,  1931. 


Furthermore, 


y,  u)  -  2  IT”.— rTTr,'’"  ■  2  t  = 


1-9 


(1+0  ^  •'-')/  -fe”*  ?x  -os  Pa; 

4 


I 


3  lP>^]~‘  , 

t-9  ' 


J  ((+-/)  ^  --»(l-')jr  ^-»(«+0x 

+ 


MUO 


Kl-0  Ml-0 


NI+0 


—  -i-f***<*Px  sin  px  +j«<»Pa  cos  p.v); 

2P 


:  7  0 

\  •  rfC  / 


•-0  «mO  * 

_  J-  f-L  ^  (J+O  ^ - L  fp  (I  -0  jr - L  g-f  (1-0  -r  1  _L  ^p  (1+0  x]  ^ 

4p*  [  2/  21  21  21  \ 

-j^ypufcpxslnpx; 


-.i[- 


3  IpK' 

^(o-s^-24-4' 

^  *  #-0  * 

Jc+'O^  ,P(I-Ox  ^-P(l+Ox 

+  -^-7-: — 4 


1+/ 
1 


1+/ 


l+i 


l-( 


Px  sin  px -«<•#»  3a  cos  Px). 


Tho  normal  fundamental  raatrijc  of  equation  (7.1)  is  determined  on  tie  basis  of 
relatlonshios  (3.1).  It  has  the  form 


in^)i 


y.(^)  ' 

>'.(^) 

y,(^) 

y,(^) 

--HP*K,(x) 

y,(^) 

K.(x) 

K,(x) 

-4P*K,(x) 

-4P‘F,(x) 

I'o(^) 

K.(x) 

-4P*K,(x) 

-4?‘K,(a) 

-4?‘K.(a) 

K«(x) 

(7.3) 


The  oarticular  solution  of  equation  (7.1) 


*•(■*)  =  ]  5  +  a) /(*)'&• 


(7.4) 


If  onto  the  shell  in  the  section  x=ai  alon^,  a.  circular  contour  there  are 
aoplied  concentrated  bending  moments  M,  and  forces  P,.  then  in  this  section 
there  exists  a  discontinuity  of  second  and  third  derivative  y(x) 


0 

0 

D 

.  D 


Fip.  5.  Plane  problem. 


Solution  of  equation  (7.1)  is  exnrr-^sa  i 
by  equality  (f>.10),  where  matrix  [y'(x)] 
taken  from  relationship  (7.3). 

If  shell  in  the  section  x=fl<  ^‘as 
reinforcing  diaphragm  (disk),  then  the  discontinuity  will  be  dependent,  in  which 
the  matrix  of  the  discontinuity  will  correspond  to  the  matrix  of  elastic  support. 

The  following  example  refeis  to  one  problem  of  theory  of  elasticity. 

We  shall  consider  the  generalized  problem  on  a  plane  strained  state  .  Wc* 
assume  that  the  external  forces  act  in  the  middle  of  plane  of  a  thin  plate  of 
variable  thickness  (Fig.  5). 

Equations  of  equilibrium  of  element  of  plate  will  be  such: 


,  d(X,yh) 

— +— -+p*Ar-o, 

0jr  dx 


(7.5) 


The  more  general  statement  of  problem  will  be  used  later  on. 


£9 


In  assuming  that  mass  forces  are  absent  and  in  introducing  the  function  of  stresses 


1  a* 

h 

1  d*a' 

A  dx*  ’ 

1 

X **  '  • 

^  h  dxdy 


we  obtain  from  condition  of  continuity  of  the  deformations 


<«  f  I  /  ^ 


dy'  /  J  aylfAl 


+  2 


a*  / 14(«  O'*  N 

dx  dy  \  £ft  dx  dy  j 


&** 

dy* 

=  0. 


(7.6) 


For  a  nlate  of  constant  thickness  with  constant  oarameters  of  elasticity  and 
with  an  absence  of  heating  we  arrive  at  the  well  known  biharmonic  eqxoation 

V*4>(jr,  y) *=0.  (7.7) 

In  investigating  the  state  of  strain  in  beams-walls  the  solution  of  equation 
(7.7)  is  sought  in  the  form 


(x,y)  = 


tw. 

cos  oty 


In  introducing  the  value  ^(x,  y)  into  equality  (7.7),  we  obtain  the  following 
differential  equation  for  the  function  ♦(*): 

(x)  -  2«»+<*>  (x)  +«♦<;.  (x)  =  0.  (7.8) 

The  characteristic  polynomial  has  two  roots  of  second  multiplicity 

X,=  a.  li=— »' 

For  determinating  q‘s(x)  ^  shall  use  formula  (2.11). 

We  have  at  a  --  0 
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— - — flo 

(Jl— 4,)*  ■  ■  '  ' 

<-0 


-^1  4 


i-». 


*-— I — * - 1  ^  2—{ixtmkax—tmk*x). 


Further  on  basis  of  equalities  (3.2)  v«  find 


V,  (x)  =  (3*»**  a.V  —  3ixt»»kax), 

Vf  (x)  *=  (2i*»fc  ax — axiinhax). 


The  aoolication  of  normal  f\indamental  functions  introduces  into  the  considered 
oroblem  a  number  of  simolifications  in  satisfying  the  boundary  conditions.  The 
use  of  the  general  method  in  Sec .  5  nvakes  it  oossible  to  construct  a  solution  in 
the  nresence  of  a  discontinuity-like  variation  in  the  thickness  of  the  olate,  since 
for  each  of  the  sections,  equation  (7.8)  remains  in  force. 


8.  Normal  Fundamental  Functions  of  Euler's  Equation 
Vfe  now  consider  the  homogene-is  Euler  equation 


or 


Jt 

a 

(c,  =  l,  y®’(.v)=>(.'f)). 


(8.1) 


(8.2) 


idiere  ct —  constant  coefficients. 


3.7 


substitution 


equation  (8.2)  reduces  to  an  equation  with  constant  coefficients: 


(8.3) 


>(c)  =  0 


(8.4) 


The  characteristic  polyncmial  of  equation  (8.4) 

(8.5) 

The  coefficients  Pi  are  the  simplest  of  all  to  determine  by  usin#;  the  identity  of 
characteristic  polynosiials  of  equations  (8.4)  and  (8.2),  if  there  ij  introduced 
into  the  latter 

y(x)«x\ 


In  decomposing  the  characteristic  polynomial  of  equation  (8.2)  by  degrees  of  X,  we 
obtain  the  values  of  Pi- 

Suppose  there  is  given,  for  examole,  an  equation  of  foiurth  order 

+  ^y*)  (.K)  +  ^y»)(x)  +  (x)  +  (x)  =  0. 


The  characteristic  polynosdal  of  this  equation 

F(X)  - 1  (1- 1)  (X_2)  (X -3)  +  X(X- 1)  (X- 2)  f,  + 
+  ^(X— l)c,+Xr,+r4. 


After  decomposition  by  degrees  of  X  we  obtain 

/?  (X)  »  X*  +  >  *  (  -  6  +  O  +  0 1  —  3^1  +  + 

+  X  (  —  6  +  2C|  —  +  ^t)  "H  P«* 

Thus,  we  find 

2i»,y-'»(o=o. 


p.-l. 

y,=  _6  +  Cj. 

p%^  11  — 3c,  +  Cyf 

— 6  +  8c, — C|  +  Cj, 


P4-=^4. 
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Wb  shall  consider  the  solution  of  equation  (8.1)  at  the  values  (a>0 

since  x  =  0  is  a  sinpular  ooint  of  the  equation).  Suppose  the  value  x  =  a 
corresponds  to  5  =  a,  a=lna. 

For  equation  (8.4)  on  basis  of  indicated  earlier  formulas,  there  may  be  found 
the  normal  fundamental  functions  and  then 

>(E)=  j'y*»(>)j'.(o. 

in  which  m 

Then  arises  the  question,  will  the  functions  Fk(e),  if  we  assume  in  them  t=lnx. 
be  the  normal  fundamental  functions  of  equation  (8.1). 

It  is  possible  to  establish  that  this  will  take  place  for  equations  up  to  the 
second  order  inclusively  and  it  is  found  invalid  for  equations  of  higher  orders. 
Hovwver,  in  any  case 

and  therefore  even  for  second  order  equations  functions  ^'•(•’t)  and  ayi(Jt).  will 
form  the  normal  fundamental  system.  If  there  is  an  inhcmogeneous  Euler  equation 

(8.6) 

then  by  substitution  (8. 3), it  reduces  to  the  equation 

Vp  y"-')(t)=y  (^),  (g  y) 

1-0 

general  solution  of  which  has  the  form 

jf(?)=  2V*»(»)  K.(;)+J  fa)/(s)t/s 

*-0  ' 


1  i  =  k 
0  i*k. 


(8.8) 


Fig.  6.  Flexure  of  disk  with  a 
ccsinusoidal  load. 


Lat  US  consider  the  following  examnle  (Fig.  6).  The  differential  equation  of 
the  flexure  of  a  circular  nlate  of  constant  thickness  under  action  "of  cosinusoidal 
load  h 'S  the  form 

cos  6, 

where  =»  -j_  _L  ^  ,  2_  d»w  _ 

r*  ()9«  /■>  d/-»  dS*  ^  r  dr> 

_ 1,  _|_  J  <^u'  __  I  <?'t»  I  dw 

r*  drdi^r*  d9'  /■>  dr*  ^  ^  ’ 


D  —  cylindrical  rigidity. 

In  assuming  a  solution  in  the  form 

wir,  6)«=9(0  cos  9, 

we  obtain  for  9  (r)  the  !i)uler  equation 

f (0  +  -7  » -  4  ?”’  ('•)  +  1 1"  ’  (O  -  -  -  T  (r)  =  . 

"  r*  r*  D 


By  substitution 


it  reduces  to  following  ordinary  differential  equation  of  fourth  order: 


(8.9) 


The  roots  of  characteristic  nolynomial 

A^(i)  =  i«_4X*  +  2)i»  +  4*-3 


(^  =  2).  X,  =  3 


are ; 


(root  has  second  multiplicity). 
From  formula  (2.11)  we  obtain 


-LA 

1!  dk 

[  (1i  +  1)(a_3) 

A  (6—) 


I2X*  +  4X  +  4  li. 


+ 


-  =_±(t_a)e5— - 

4Ji»~I2X*+4X  +  4  4'  ^  16  ^16 


By  means  of  equalities  (3.2)  we  find 

***1  (0  =  +  —  (;  —  a)  €'-'  +  —  g-({-») _ \_g3 

^4  16  16  ’ 

<I»,  ({)  =  -1  <«-  4-L  (•-  _  a)  gi-. - L  _  ±  ^3  (£  -,) 

24  16  16  ’ 

(0  =  - ^  -|-  J_  ^  (£-•) . 


The  narticular  solution 

i 

(0  =  J*  (E  -  5  +  a)  -?|L  e*-ds. 


(8.10) 


(8.11) 


Thus, 

I«t  us  apply  now  the  obtained  solution  for  a  plate  (disk)  of  variable  thickness. 
For  this  purpose  we  shall  divide  the  plate  into  sectors  of  constant  thickness  with 

sections  6=a/  (/=0,  I . ).  We  shall  designate  the  cylindrical  rigidity 

D  («;— 0)  =“  Dj-\, 

D(aj+0)^Dj 

and  correspondingly 


?(«>— 0)  =  9/-i. 

^(«y  +  0)-i>>. 

From  condition  of  equality  of  the  bending  moments  per  unit  of  length  of  cylindrical 
section  we  obtain 


IfJ’*  —  ( I  —  h)  tJ* ’  -^PTyl  =  I?ji\  -  ( 1  —  «*)  ?)-!  —  1  •  (8.12) 
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Owing  to  th«  continuity  of  de.formation 


and  then  from  (8.12)  we  obtain 


“  (■—  - 1)  (9)1’,  -  (I  -  p)  ?;l’,  -Kf/-,!. 
From  condition  of  equality  of  i.ie  t.Dtal  transverse  force  we  find 

Dy  If )« - 2?"’  -  (J -  !■ )  f )' ’  +  (3 -  C)  9 J  - 

-  Dj-,  lfji\  -  2f‘V, -  (2- ,,)  9)'J,  +  (3 -f) 

It  follows  from  this 


4J.) _  ^  _  1^  I- _  (2  - ,C. ,  + 

+  (3-t‘)?/-.|+23'”  =  (-^2--l)x 
X  l9)L',  -  (4  -  )  ,JL',  +  3  (1  -  I 


(8.13) 


Thus,  we  shall  have 


I?  Wl  -  ,5:  -S  (*.  «/)  I®.  (=  -«y)l  (Cyl  l?  (<j,)  1  +  1®.  (!)(■ 

/  -^v 


(8.14) 


Here  [^*(^—tt/)j  —  is  the  normal  fundamental  matrix  of  the  homogeneous  equatic 

(8.9).  Its  elements  are  determined  from  the  relationship  (3.1)  or  by  direct 
differentiation  of  equalities  (8.10). 


The  avatrix  of  discontinuity  has  the  form 


0  0 
0  0 


\  Dj  }  —  J»  _(i_p)  I  f) 

_3(1-P)  -(4-3j0  0  1_ 

In  equality  (8.14)  into  lT(<^y)l  the  left-hand  values  enter 


l?(fly)l  =  l'f/-'I' 
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Mb  shall  dwell  on  determining  the  column  of  a  oarticular  solution.  At 
from  equality  (8.11)  we  obtain 

/-I  *. 

=  2  J  ‘***(=  —  5  +  “)  e*’ds^ 


»-l  S-l 
t 


■  f  4>,({-s  +  a)^(i^ 


O.-i 

e*^ds. 


(8.15) 


or./ 


In  determining  the  elements  i‘I»*(*)]  one  should  not  forget  that  all  connonentarr^f-  . 

the/ 

hand  side  of/equality  denend  on 


9.  Application  to  the  Integration  of  Equations 
with  Variable  Coefficients 


Suppose  there  is  a  linear  differential  equation  with  the  variable  coefficient" 


2  p,(x)y-')  (.«:)=/(•') 


(9.1) 


l/»o(^)  =  l.  y»»(A)=y(x)l 

and  there  is  sought  the  solution  of  equation  in  the  interval  d<x<b.  The  mean 
values  of  coefficients  can  be  determined  for  example,  in  the  following  manner: 

out/  • 

We  write/equation  (9.1)  in  the  following  form: 

-/(^)  +  i  (Acp-/»i  (jr))y"-“  (X) 


(9.2, 


(9.3) 

Suppose  y*(x)  (fc=0 . n— 1)  -..are  the  normal  fundamental  functions  of 

equation  with  constant  coefficients : 

,E/»icpy-'’(-<^)=o  (9.4) 

(with  the  initial  section  x  =  a). 

Solution  of  equation  (9.3)  with  an  arbitrary  right-hand  side  is  written  out  as 

y  (J«) = "s  y  *’  («)  (Jf)  +  j'*  ( ')  + 1  ^ 


«  ep  equals  averafe. 


X  i  (Pi cp-A»i(*))y (9.5) 


wh«r« 


(9.0) 


>'•('«)  •=  j  >'<.-1  (.v-s  +  fl)/(s)(ys. 


In  diff«r«nti»ting  equality  we  find 


y>(x)  -  (a)  Ki"  (A)  +  K;-'  (A)  +  jKl-j,  (  V  -  s  +  a)  X 

(—0. 1 . "-I) 


(V.v) 


R>'  introducing,  as  nreviously,  the  columr.-soiution 


y(^) 

yO>(x) 


(9.8) 


we  shall  write  out  system  (9.7)  in  following  fom: 

[yWI-(?(A))U(a)|  +  |r,(v)i+j(/c(A,  i)||>.(s)|*. 

a 

"Equation  (9.9)  represents  the  matrix  integral  Volterra  equaLion. 
The  matrix  is  the  nucleus  of  the  equation 


(a:(x.  s)j  = 

K.-,  (x-  j  +  a)  (p^ ~p^  (s)). . . K,_, (x  - s  -f  a)(/»,,p -Px{s)) 

*  * . .  ■  •  • 

LKllT‘’(A-,  +  a)(/,,„-/,,  (s)), ,  .?!,"--'’(x-j  +  a)(/,„, 


(9.10) 


Hlquation  (9.9)  is  solved  by  method  of  successive  anproximations,  in  which  the 
orocess  is  convergent,  if  all  coefficients  Pt(x)  are  limited  in  interval 
a<x<b.  Furthermore,  it  is  possible  to  show  that  solution  of  equation 

(9.9)  results  in  determining  the  normal  fundamental  functions  of  the  homogeneous 
aquation  (9.1)  and  the  narticnlar  solution,  satisfying  the  zero  initial  conditions 
Sunoose  it  is  necessary  to  determine  the  k-th  normal  fundamental  function  of 
aquation  (9.1)  Yk(x)- 

We  shall  assume  in  the  column  of  initial  values  all  ^^’’((i),  ('=^fr)  are 

equal  to  zero,  and  v<*>(a)=l.  From  equation  (9.9)  at.  y'*('-)="0 

cor, sequent  Iv,  /(x)— 0, 


we  obtain 


(9.11) 


M 


The  first  apnroxioation  is 

1>'*,„(ac)I=(k,(x)]:  (9.1^) 

second  approximation  is 

(K.,„(jr)l  =  lK*(x)l+Jl/f(x.  s)l[K,(s)lrfs 

A 

et  cetera. 

The  nresentation  of  solution  in  the  fonr.  of  series  leads  to  the  same  result-' 


!(«!•,  (x)l. 


(9.14) 


in  which 


(9.1i)) 


Calculation  is  terminated,  when  the  difference  between  two  successive  aoproximati''' 
can  be  assumed  to  be  negligibly  small  or  with  the  use  of  equalities  (9.14),  when 
considered  terms  of  the  series  is  small  in  comoarison  with  the  sum  of  preceding 
terms.  In  the  process  of  successive  approximations  there  is  no  necessity  to 
calculate  all  the  elements  of  the  column  [yk(x)];  onl^'  the  elements,  for  which 
Piep — Pi(s)  $k  0.  will  be  subject  to  calculation.  Others  required  t 

calculating  the  derivatives  1^»(-*).  are  comouted  from  equation  (9.11)  after  ^he 
indicated  elements  have  been  determined. 

Me  shall  consider  as  examnle  the  Bessel  equation  of  order 


>(X)  +  >(JC)  +  (»  -  (^)  =  0- 

We  shall  seek  an  anoroximate  solution  of  this  equation  in  the  interval 
The  mean  values  of  coefficients  by  formula  (9.2)  are  equal  to 


(9.16) 


Pi  cp  ■*=  -;p  In  6, 
5k 


(9.17) 


*,  f 


Th«  normal  fundamental  functions  of  the  equation  vdth  constant  coefficients 

(9.18) 

will  be  such: 


1 


~  Pjcp 


J'lw- 


(9.19) 


whore 


Php 

2 


/’itp 

I 

ruf 


Pitf  • 


/’icp  • 


Tn  the  considered  case 
In  considering  (9.17),  we  find 


a  = 


~  1  -  0.3  4?? 

K,(x)  =  — i-r 
•  0,8549 


0.854‘J  COS  0.8549  A  - 


0-:) 

V,  ^ 

0,3422  sin  0.o549  {x  -  y)] 


K,(a)  = - e 

0,8549 


1 


fie  shall  determine  the  normal  fundamenial  functifjn 
From  rolationshin  (9.11)  we  shall  have 

y.(^) 


sinO,8M9(A-yj 

yo(v) 


of  equation  (9.16), 


y.  (0 

J 

AJ'’W  . 

+ 


K,  (X  -  s  +  a)  (ricp  -  y).  Yi  I'  -  s  +  a)  +  ^t) 

[^Kl'>(A-54-a)^p,^— Y  j.  Kl'>(A-s  +  o)^p„p- I 


K,(s) 


{/s. 


0 


H«nce  the  first  approximation  for  the  function  ^ofA) 


Second  anproximation 


^'o(i)("r)=  PoW- 


Ko  (2 )  (JC)  =  K,  (X)  +  J  K,  (X  -  s  +  a)  X 


(9.20) 


(9.21) 


In  Table  1  the  calculation  of  second  annroxiiration  is  exolained.  In  the 
calculation  the  interval  '  is  divided  into  ten  equal  sectors.  For  comnilinp, 

tables  of  values  of  the  functions  y',(x— s  +  a)  is  calculated  as  soon  as  first 

column. 

Since  s  =  a,  then  into  the  column  tne  values  y'i(x).  are  entered.  All  othei 
columns  are  filled  in  such  a  manner  so  that  the  elements  on  secondary  diagonals  are 
identical. 

This  ensues  from  the  circumstance  that  in  the  division  into  equal  sectors  for 
these  diagonals  the  magnitude  .v— s+a  remains  constant. 

In  column  12  there  is  introduced  an  exnression,  standing  in  equality  (9.21)  ir 
brackets.  Furthermore  each  element  in  column  12  is  multiplied  by  the  corresnondin, 

value  _ s+a)  result  is  summarized  by  the  traoezoidal  rule 

for  all  values  in  a  given  data  line.  This  sum  is  contained  in  column  13.  In 
column  14  there  is  given  the  value 

Y»  »n  (•*■)  =  Y*{x), 


then  —  value  of  the  magnitude  (x),  F«(,»(x)  accurate 
value  V*(x).  In  the  given  case,  a  third  approximation  gives  an  accuracy, 
adequate  for  a  majority  of  engineering  problems  (Fig.  7). 

SooBtimes  it  is  expedient  to  apply  another  variant  of  the  presented  method. 


which  consists  of  the  circumstance  that  equations  vd.th  constant  coefficients  strive 
to  obtain  a  form  as  simple  as  possible. 


Zh 


Table  1.  Solution  of  Equation  (9.16) 


Wb  shall  turn  to  the  generalized  problem  of  a  olane  stress  condition  and  shall 


assume  that  h  — h(*)  (3ee  Fig.  5).  Then  for  <•(')  we  obtaij’  the  equation 

(jc)  (x)  =0. 

iirtiere  ,  .  _  rf  i+|i 

,,(x)_2-(£A)— 

By  transferring  the  term,  containing  <{»('> (a),  into  the  right-hand  side  of  equalitj  . 
ws  arrive  at  the  equation 

f  (X)  =  (0)  (A)  - J  V-,  (X  -  s)  p,  (s)  4.C. ) (5)  ds. 


where  __  is  the  normal  fundamental  functions  of  equation  (7.6). 


Fig.  7.  The  Comparison  of  accurate 
and  approximate  solutions 


The  first  approximation  in  this 
oroblem  expresses  the  stress  condition  in 
a  olate  of  constant  thickness. 

In  conclusion,  we  shall  say  several 
words  on  another  method  of  using  equations 
with  constant  coefficients  for  solving 
equation  (9.1),  which  sometimes  is  annlied 
in  engineering  problems.  In  this  method, 
the  total  interval  of  variation  of  x  is 


divided  into  sectors  within  limits  of 

which  the  coefficients  Pi(x)  are  assumed  constant.  And  here  the  use  of  the 
solution  in  form  (1.7)  gives  an  essential  advantage. 

Thus,  determination  of  [y(^)l  reduces  to  the  multiplication  of  matrices, 

corresponding  to  individual  sectors. 


*  To  accurate  value. 


CHAPTER 


NORMAL  FUNDAMENTAL  FUNCTIONS  OF  LINEAR  DIFFERENTIAL 
EQUATIONS  WITH  VARIABLE  COEFFICIEfJTS 


There  are  considered  two  methods  of  approximate  determination  of  normal 
fundamental  functions:  method  of  successive  approximations  and  method  of  linear 
approximation. 

Both  methods  are  used  for  solving  the  normal  integral  equation,  to  which  a 
differential  equation  under  initial  Cauchy  conditions  reduces. 


1.  Statement  of  Problem 

There  is  given  a  linear  differential  equation  of  n-th  order  with  variable 
coefficients 

It  is  required  to  find  the  solution  of  this  equation  in  certain  interval  of 
variation  of  x(a(x(b). 

The  functions  /'i(A)(<“l . n)  and  / (x)  are  assumed  limited  in  the 

indicated  interval. 

The  set  of  n( linearly  independent)  solutions  of  the  homogeneous  equation  (1.1) 
{Fj(.<)),  A  =  0.  I . a— 1.  satisfying  the  condition 


1  =  k 
i’^k 

{i.  k^O.  1. 


•  .  1), 


(1.^) 


13  called  the  normal  fundament.al  3  imu  1  lar.e  Tn.)  eqiial:^)’i  (1.1)  wilt,  the  ir.itlal 

sect. inn  x  a. 

If  it  is  known  that  (a)  —  i;  a  riartiriiar  3>ilution  ">*'  equatio;,  (1.1), 

sat  isfyin^t  zero  init  ial  condii  inns 

K;-'(u)=0,  /-  0.  1.  .  .  .  n  1. 

then  the  solution  of  the  equatior  is  nresentel  as: 

where  >'**’ (a)  values  of  function  ,v(x)  and  its  first  n-1  derivatives  In 

section  x  -  a. 

We  shall  now  nrnve  the  converse  assertion.  If  t.he  solution  of  equation  (l.j 
under  arbitrary  initial  conditicr.s  a.nd  if  the  arLi'rary  function  j'(x )  can  te 

presented  in  the  form 


(1.0 

(1..0 


(1.5) 

where  -  0  at  /  (x)  0,  then  functior.  fbe  normal 

fundamental  functions  (homogeneniis )  of  equation  (1.1  ) 

Z*(0  (1.0) 

and  function  Z,  (x)  is  a  particular  solution  of  equation  (1.1),  satisfying  c.er 
initial  conditions 


^*(-0  .  (.V). 


(1.7) 


For  proof,  let  us  assume  at  first  /(v)=0 
conditions  in  sucn  a  form: 


yu)(a) 


I  0  /  /  A: 


(i.-O 


Tnen,  from  re lationsnin  ( 1 .  i> ) 

y  (x)- 


By  differentiating,  we  fi.nd 

^o'(A)  :y*’f(j)Z;''(v). 


In  view  of  tho  dopondonci  of  (1.8)  at  x 


a  we  now  obtain 


Zi-.  (.)-{; 


which  orovee  equality  (1.6). 

Supooae  now  the  initial  conditions  are  zero: 

yo(a)-o  (/-O . s-l).  (1.9) 

Then  fron  relationship  (1.5)  it  follows 

y(x)-Z*(x). 

but  in  view  of  dependence  (1.9)  equality  (1.7)  proves  to  be  valid. 

We  note  still  a  subsequent  result.  If  under  the  above-indicated  conditions 

y (x) - "s  y *’(«)  «»*  ( Jf)  +  «•  (^) .  (1.10) 

then 

(1.11) 

2.  Determination  of  Normal  Fundamental  Functions 
and  of  a  Particular  Solution  by  the  Method 
of  Successive  Approximations 


In  designating 

(2.1) 

and  in  considering  relatlonshlo,  being  obtained  after  repeated  integration  of  the 
equality  (2.1), 

y*’  (x)  -  2  y (")  + 

•  • _ m 

•-»  times 

(*—0.1,.  .  .,«-l), 

e  shall  from  equation  (1.1) 

f-A/f +'r  y*>(a)A(x)+/(x).  (2.2) 

*-e 


^6 


where 


/vt-- 


i-i 


*i-\ 

.  .  f  (  t^)  dx,...  dxi, 

_ m 

‘  tijnee 


(2.3) 


IrnH-M 


(/-/!+*)! 


(2.4) 


At  x->a  equation  (2.2)  reverts  to  an  identity  by  virtue  of  equality  (1.1). 

/  \  ^ 

Equation  (2.2)  reoresents  the  normal  integral  equation  . 

Cursa  established  it  in  another  form,  the  solution  of  equation  (2.2)  can  bu 
nresented  as : 

T(->f)“"2V*(a)«»*(*^)  +  <I>*(^).  /P 


In  this  equality 

®.W“/.+AV.+JV*/,+  .  .  .  -S/V'A.  (2.6) 

where  /V'A-  signifies  s  times  the  reneated  application  of  operator  N,  where 

(/v/.-/.). 

Correspondingly 

♦.W“/+Ay+«’/+  -  2A?>/.  (^  .7) 

The  series  (2.6)  and  (2.7)  converge  uniformly  and  absolutely. 

By  virtue  of  equalities  (1.11),  formulas  (2.6)  and  (2.7)  solve  the  nosed 
nroblem. 

For  a  determination  of  function  K*(.<)  ®nd  its  derivatives  one  should 
use  equality 


♦Theory  of  normal  integral  equations  is  considered  in  Chanter  3. 

♦♦E.  Cursa,  Course  of  Mathematical  Analysis,  Vol.  3,  3tate  Theoret.  Technical 
Publ.  House.  Moscow- Leningrad,  1934 


J  .  .  .  J  ,  .  .dxt  (<>*) 

••ttimii 

l  +  J  -  •  •  J 

«-*  tUMS  (^■■^) 

M  0.|.| 

<*"*!*-  + f  .  .  •  f  d>,ix,.,)dx^.i .  .  .dxt 

:  ;  («A) 

(<"■0.1,.  .  .,#1—1). 


(2.8) 


Equal! tl«s  (2.8)  make  It  possible  to  determine  all  the  elements  of  e  normal 
fundamental  satrix  of  equation  (1.1). 

Now  we  preaent  an  example.  The  equation  of  stability  of  rod  of  constant 
section  on  two  end  knuckle  bearings  has  the  form 

y<*>(x)  +  aV(Jc)«0.  (2.9) 

The  initial  section  a  »  0; 

y*>(x)-T(x). 


Prom  equalities  (2.3)  and  (2.4)  we  obtain 


yVf--  ypy(x)J  J  9iXi)dx,.  .  .(/jr,= 


l-l  '  i  i 


I 


00 


/•w— S'-wW 


_a» 


a*. 


t 


Equation  (2.2)  will  be  such 


<-i 


•  jit 


9  (.<)—-»•  J  J  f  (-‘fi)  dx^  dx^  -y  (0)  a* — jf* '  >  (0)  a-x. 

U  0 


According  to  equality  (2.6) 


♦t  ( jt) = n=’(^) = — ’ + +  ■ 

I  /l  I  \ 

('-5r+-47-  •  •)• 


/  a*jcJ  ,  a»-»* 

—  a  a.t - + 

\  3! 


5! 


For  equation  (2.9)  the  accurate  solution  is 


and  furthermore, 


^0  (■*)  =  cos  Yi  (x)  «=  -  sin  a.v 


—  a*  cos  a.v,  K',^>(a)=  —  asInoA. 


(2.11) 


By  comparing  formulas  (2.11)  and  (2.10),  we  readily  note  that  in  given  car,-:, 
each  term  of  the  series  (2.6)  represents  a  corresoonding  term  of  the  expansion  of 
accurate  functions  into  a  oower  series. 

We  shall  make  several  remarks  of  a  nractical  nature. 

The  magnitudes  of  A/*/*  should  be  determined  by  using  anproximate  methods  '^f 
calculating  the  integrals  of  which  the  simplest  is  trapezoidal  rule. 

In  calculating  the  iterated  integrals  with  a  variable  upper  limit,  an  essential 
decrease  in  computing  work  is  obtained  by  subdividing  the  interval  into  sectors  of 
equal  length  and  use  of  the  'Ring  rule'  (Table  2). 

Into  the  column  are  entered  three  numbers,  being  encompassed  ly  the  arrow;  for 
obtaining  the  true  magnitude  of  integral  the  values  in  the  column  must  be  multinlied 
by  ij.  where  ^ -is  the  length  of  sector,  n-is  the  number  of  integration 

operations . 

For  determining  "  I.  2,  3,...)  the  integral  oneration  N  is  completed  or 

function  A/’  '/*,  the  values  of  which  already  are  in  corresnonding  column  of 

{N»/,  =  fk). 


the  calculating  table 


Ttbli  2.  DlagrAB  of  InUgratlon  on  basis  of  ths  Ring 

rult 


f(*) 

BH 

7^ 

(i^)’ 

(UY 

Xg  m  a 

'•r 

"  V 

^0^ 

■X  0 

( 

EE 

<  □ 

T 

Eg 

<  o 

^  a 

KEY:  (a)  factor. 

3.  Discontinuous  Solutions 


Lst  us  assuBS  that  in  the  solution  of  equation  (1.1),  there  are  given,  in 
addition  to  the  initial  conditions,  discontinuities  of  first  order  of  function  y(x) 
and  its  n— 1  first  derivatives  in  the  sections  (/'*“  part  of  given 

discontinuities  can  have  zero  values. 

Vie  shall  designate 

y’>(fl>  +  0)— y<”(Oy-0)»= 

/»-0.1 . «-l,\  (3.1) 

\y-15 . m  /* 

The  initial  values  also  may  be  considered  as  the  given  discontinuities  after 
assuming 

y(a)«A(0) 

••••••  (3.2) 

y  •-!)(«).  A 

We  shall  introduce  a  single  discontinuous  function  S(x,c),  determinate  by  the 
equality 


So 


(3.3) 


x<c, 
X  >C. 


If  ^  (J^)  is  the  arbitrary  integrand,  then 

J  *5 (-^i.  f)  /  (-*^1)  ^S(x.c)j/  (x,) 

0  t 


and  furthermore 


//(*  ■■  j  .  .  .rf.t,  = 

mass 


•  .  y/(x.)dx,.  .  .dx^. 

rtimee 


(3.4) 


13.5) 


In  particular,  at  /('^)“‘' 

'  •'.-i 

J  J.  •  •  J  .  .  .  dx^^S{x,c)  ~  ^  - .  (3.6) 

If  f(Jt)-"  is  an  arbitrary  differentiable  function,  then 

^(5(^.0?('«)l=«^('>f.0  “O':)-  (3.7) 

For  a  discontinuous  function  we  shall  have 

y{x)=  2  AJ'»5(jf,a;)  +  Jj,{i)(;c,)(/^,  (3.8) 

0 

and  further 

Jf<'>(x)»  5Ay>5(x.a,)  +  yy2»(x,)(/x,.  (3.9) 

/»o  0 

By  introducing  (3.9)  into  equality  (3.f’)  and  by  using  dependence  of  (3.4), 


obtain 


>(x)-  2  4(0)5 (.V, a,)  +  i:A(«»5Ua^)(;c-a.)4- 

/-o 

+  Jj>>'(x,)dA:,dx.. 

•  0 


(3.10) 


Successively  by  anolyinp  this  metnod  we  shall  find 

«-  1  m 


*-o  y-4 

X  ■'*-1 


+J  J  •  •  •  J 


(3.11) 


i/ 


we 


*  tijsss 


In  (llfr«r«ntl4tinf  •quality  (3*ll)f  m«  •stablish 


+  J  J  .  .  .  I  ,) 

•  •  *  ■ 


(3.12) 


•  xttMt 


(vwO,lt  .  .  .  irt  1). 


In  daaignatlnc 


and  by  introducing  equality  (3.12)  into  equation  (1.1), we  shall  obtain 

f  -  AT, +^5  *2^  (X)  +  /  (a). 

where  Nf  is  given  by  equality  (2.3),  and 

Equation  (3.13)  is  a  nornal  integral  equation,  equivalent  to  the  differential 
equation  (1.1)  together  with  initial  conditions  and  given  discontinuities  of  the 
function  itself  and  its  n—l  first  derivatives  in  the  sections  (/*1,.»,^)- 
If  the  Indicated  discontinuities  are  absent,  then  equation  (3.13)  and  (2.2) 
agree  by  virtue  of  equalities  (3.2). 

The  solution  of  equation  (3.13)  will  be  such: 

tW-  S  +  (3-15) 

y.o*-o 

Khere  the  fiinction  are  determined  by  the  equality 

•wW-/./+A^/„+/V’/„+ .  .  (3.16) 

and  the  function  ♦•(•’f)  is  given  by  the  orevious  formula  (2.7). 

It  is  nossible  to  establish  that 

/il-=0 . «-l\ 

ly-0.  .  .  .  .m  j- 

where  K*y(.t)  — is  the  k-th  normal  fundamental  function  of  the  homogeneous  equation 
(1.1)  %d.th  the  initial  section  .v«=a^. 

For  an  equation  with  constant  coefficients 


(3.13) 

(3.U) 


SZ 


(3.18) 


Normal  fundamental  functions  Y„(\)  an  initial  section  x  a  can  be 

designated  now  (ik-0 . /,-!). 

All  other  functions  J'tyC.v)  can  be  represented,  as  linear  combinations  of 
the  functions 


where  the  constant  coefficients  qi,j  are  determined  from  n  equations  n 


i=k, 

ii-k. 


The  exnediency  of  the  method  of  determining  y  (a)  j'^\  Z”^y  means 

of  the  series  (3.16)  or  the  equality  (3.19)_7  ^  determined  by  peculiarities 

the  nroblem 

Solution  of  equation  (1.1)  with  given  discontinuities  of  function  y(x)  and  its 
n-1  first  derivatives  in  m  sections  =  has  form 

>(•<:)=  V  (3-20) 

/•iO  *.»  0 


U.  Annlication  of  the  Linear  Approximation  Method  . 


V/e  discuss  now  another  method  of  determining  normal  fundamental  functions, 
based  also  on  the  solution  of  equation  (2.2). 

For  determining  <tj^(A)  =  y*"’ (''^)  there  is  solved  the  equation 


(4.1) 


<K)ther  methods  of  linear  approximation  will  be  considered  in  Chanter  3. 


5? 


whtr«  thf  solution  of  •quatlon  ^(.v)agr«e8  with  the  function  (In  determining 

•  ptrtieuUr  solution,  tttisfying  ssro  initial  conditions,  wo  procssd  from  the 
•quation 


Wi  shall  divide  interval  of  variation  of  x  into  a  number  of  small  sectors  and 

will  designate  the  boundary  sections  . .  .  .  Within 

limits  each  seetor  ws  shall  assume  the  function  (.v)  as  linear. 

For  first  section  (a‘,<jc<a'}) 


where 


A|-4r,  * 


In  introducing  values  of  ?(a)  into  equality  (4.1)f  we  shall  obtain  at  a'>-.V} 


M 


<-l 


or,  by  introducing  the  value  k,, 


Ti- 


\ 


Pn 


hml 


(■>1  -<«<>)' 


•  (4.2) 


Her.,  4nd  htncforth  th.  following  abbreviated  designations  are  used; 


/*  (•*))  "  Ay 

The  equation  (polygonal)  of  the  function  ©(x),  valid  within  the  limits 
can  be  exoressed  in  the  following  manner: 


where 


t(J<)“t,+ j'  |S(jr,Ary)*y(v-X;)- 

/-• 

— S(x,Ay+i)*y(A  — (4.3) 

— — —  p 


and  the  single  discontinuous  functions  are  determined  by  the  equality  (3.3). 

In  Introducing  values  of  T(x)  from  equality  (4.3)  in  equation  (4.1),  we  shall 


obtain  at  jca*x. 


f. 


(>“1.2. 


•  .r). 


(4.4) 


where 


«•/- 


M 


i-\ 


In  the  latter  equality 


(/+!)![ 


For  j  -  we  shall  have  other  formulas: 


U.5) 


(4.6) 


PiKO) 


,TnF[('+  J- 


(4.7) 

(4.3) 


Let  us  note  that  coefficients  a,/  remain  identical  during  determination  of 
all  fundamental  functions  and  of  the  particular  solution. 

Of  prime  practical  interest  is  the  subdivision  of  interval  into  equal  sectors. 
Here 


X,-Xy  =  A(v-y), 


*irhore  4  is  the  length  of  a  section. 

Equalities  (4.6),  (4.7),  and  (4.8)  now  will  acquire  the  form 


(/+>)» 
Pi(^O) 


o+»v 

a* 

‘  (<+»)! 


(4.9) 


(4.10) 

U.ll) 

(4.12) 


The  latter  formulas  are  conv'eriientlv  nresented  as 


wher*  the  coefficients  '[l(^y)  remain  one  end  the  seme  for  any  (linear) 

differential  equations. 

Wi  now  preeenti  as  an  example,  these  coefficients  for  an  equation  of  fourth 
order  end  for  subdivision  of  the  interval  into  ten  sectors  (Table  3).  The  matrix 
of  coefficients  hi^J)  is  triangular.  All  elements  of  the  main  diagonal  are 
separated  by  a  heavy  line,  equal  to  unity. 

Elements,  standing  in  secondary  diagonals,  with  the  exception  of  those  belonging 
to  tho  first  columii,  are  identical;  therefore,  calculation  will  subject  only  the 
elements  of  first  two  columns.  The  tables  remain  valid  also  in  the  subdivision  into 
a  smaller  number  of  sectors. 

Below  there  is  presented  a  diagram  of  the  calculation  (Table  4).  The  first 
oart  of  the  table  contains  the  values  a,y,  in  which  the  coefficierits  a„  are 

differentiated  by  a  heavy  line.  The  values  /*.  are  known.  At  >*=0  ?o ;=/»«• 
Furthermore  the  magnitude  is  multiplied  by  elements  of  column  0  in  the  table 
and  is  entered  into  the  column  0  of  table  3./7y. 

For  obtaining  there  are  summarized  all  the  terms,  standing  in  the  line 

and  they  are  divided  by  the  magnitude  1 —  After  obtaining  'f’  column  is 

filled  in  the  table  et  cetera. 

The  values  (**=0.1.  •  •  •  ,s— 1)  are  determined  from  relationships  (2.8). 

3ince  here  the  integral  operations  are  ccmoleted,  then  the  accuracy  of  calculation 
increases. 

In  the  considered  methods  of  approximate  integration  of  differential  equations 
there  exists  an  effective  method  of  verification. 

Tho  aoproximate  value  9  (x)  is  introduced  into  the  integral  of  equation  (2.2) 
and  there  is  determined  the  difference  between  the  left  and  right  sides  which 
constitutes  the  error  of  the  solution  •(■»). 

By  integrating  tCx)  in  accordance  with  equalities  (2.8),  we  find  error  in 
•'unction  Vktx)  and  its  derivatives.  From  nhysical  considerations  there  is 


Table  3.  Table  of  coefficients  7/ (».;)  (for  brevity  the  J  coIijtjis  have  n  5t  been  extended) 


established  the  admissibility  of  auy  one  error. 


errors/ 

V/e  note  that  a  majority  of  engineering?  problems  /  of  an  order  of  2  to  yfj 
are  entirely  admissible,  since  they  correspond  to  the  acciu'acy  of  the  given  ir.itial 
magnitudes . 


5 .  Met  hod  of  Moving  Origin 

The  effectiveness  of  nrcviously  nresentel  methods  decreases  with  an  increase  i 
the  length  of  interval  of  x  variatior.,  in  the  elongat.ion  of  which  the  solution  is 
sought. 

This  circumstance  is  neculiar  to  almost  all  methods  of  approximate  integrat.i')' 
of  differential  equations. 

In  using  the  method  of  successive  approximations  the  convergence  fir  large  x 
values  deteriorates. 

Thus,  for  example,  in  equalit.ies  (2.10)  there  is  obtained  an  expansio.'  by 
deg”  !es  of  x,  which  converges  at  any  x  values,  but  a  small  number  of  approximations 
gives  good  accuracy  only  at  aJC<^V^. 

For  the  method  of  linear  aonroximaiior!  the  accuracy  increases  with  a  decrease 
in  the  length  of  sector  of  interpolation,  hut  an  increase  in  the  num.ber  of  secti'Tis 
results  in  a  large  increase  of  computing  work,  the  indicated  reasons  make  expedient 
ohe  application  of  a  special  method  of  calculati or,,  to  the  discussion  of  which  w^ 
now  turn . 

We  shall  consider  at  first  the  met,hod  of  successive  approximations . 

Suppose,  for  example,  there  is  determined  the  k-th  normal  fundamental  furctior. 
i.e.,  equatio”:  (4.1)  is  solved. 

Furthermore  from  the  calculatio:  it  is  asceriu .  .ned  tfiat  v.-ith  three  to  four 
approximations  the  values  9  for  x<ai  agree  with  the  required  accuracy.  Then, 
by  means  of  formulas  (2.8)  all  y**’  (^i)  -\) 


are  found.  Now  it  is 


nosslble  to  transfer  the  initial  section  into  x  a^^  and  to  make  the  calculation 
according  to  the  equation 


vrtiere 


*-• 


(5.1) 


■  -  1  Jt 

1  pM\  .  •  f  ^(Xi)dx,.  .dxt, 

•i  *1 

— rtiMfc 


I— I 


/•  (jf)  —  2  ~  ~  — ~  • 


thus  equality  (2.2)  Is  valid  at  an  arbitrary  value  a. 

Briefly  written,  equation  (5.1)  has  the  form 

f- A/jT  +  /„. 

where 

y •)(«») /*(Jc) 

t  —  0 

is  a  known  function.  The  subsequent  stage  of  the  calculation  reoeats  the  nrecedin 
For  determining  the  particular  solution  there  is  solved  the  equation 

In  second  section 

T  -  +*  J  ><*'  (a,)/.  + 

/here  N|9  and  f^{x)  '-he  sane  n-i  in  equality  (5.1). 

We  now  establish  an  evaluation  making  it  nossible  to  determine  a^^  nrior  to 
omoletion  of  calculation.  The  matter  reduces  to  an  evaluation  of  terms  in  the 
cries  (2.6)  or  (2.7). 

We  shall  give  a  very  "rigid''  evaluation,  which  will  assure  the  condition  of 
anid  convergence  in  the  orocess  of  successive  anproximt tions . 

Sunpose 

M-max 

P^-max|/»^(x)l. 


(5.2) 


Then 


■  •  '4-1  ■  j 


l-l  m 


l-l 


(5.3) 


whore 


(,>0,-0, 


|W*/J<  i  P, J  .  ,  .  Tl^7.!rf^,  .  .  .dx, 

l»l  «  « 


<.(s«4 


(5.4) 


and  in  general 


t)  • 


Under  the  condition 


<  I 


(5.5) 


(-.1 


series  (2.6)  and  (2.7)  will  be  ranidlj.'  convergent  which  can  be  established  by 
comnaring  an  evaluation  made  more  strictly. 

From  the  latter  it  follows,  as  already  was  indicated, that  the  series  (2.6)  ar.i 
(2.7)  are  absolutely  and  unifonnl.v  cor.ve''gent  for  the  finite  values  A  and  P.,  ,  : 
which  one  is  readily’  convinced  after  nresenl inr  N  in  the  canonical  form  of  a 
Volterra  operator. 

V.’e  now  nresent  some  examples,  relating  to  equation  (2.9).  For  determining 
solve  the  equation 

a*. 


0  u 


we  shall  have 


/  =  P, 

Evaluations  of  (5.3)  and  (o.4)  give 


W.  y. 


lAPA  • 


(5.0) 


4 


condition  of  convorgence  (5.5) 


(5.7) 

The  more  accurate  evaluation 

(5.8) 

established  convergence  of  nrocess  of  successive  apnroximations  in  an  arbitrary, 
but  finite  interval  of  variation  of  x. 

At  -r**/,  ,  real  value  is  6  times  less  than  the  evaluation  (5.6). 

Condition  (5.7)  assures  a  raoid  convergence  of  the  nrocess  of  successive 
aonroximati ons  in  the  interval  a<x</,. 

However,  in  nractice,  an  accurate  evaluation  of  the  magnitude  a^^  (limits  of 
ranid  convergence),  is  not  required  since  it  is  ascertained  in  the  process  of 
calculation  with  an  arbitrary  value  a^;  it  may  be  found  onl;/  that  seme  of  calculations 
at  large  x  values,  will  anpear  to  be  sunerfluous. 

This  is  virtually  established  after  the  first  two  anproximations . 

We  note  that  selection  of  magnitude  a^  and  "lijiuting''  number  of  utilized 
annroximations  is  determined  by  peculiarities  of  the  operator  Obviously,  the 
more  simple  the  structure  of  the  operator  is  the  greater  number  of  anproximations 
can  be  applied  and  the  rarer  can  the  transfer  of  origin  be  used. 

Method  of  moving  origin  withe  t  any  changes  is  extended  to  the  method  of  linc'ir 
approximation,  where  the  use  of  more  than  ten  sectors  becomes  unwieldy. 


6.  Quasinormal  Fundamental  Functions 


The  differential  equation  of  the  type 


dx"' 


Pi  ('<)-;..>'(■<) 


(-«) 


(6.1) 


we  shall  call  an  inhomogeneous  binomial  equation;  correspondingly  at  f(x)  -  0,  it  is 
homogeneous.  To  equations  of  this  type  belong  a  large  number  of  equations,  en¬ 
countered  in  engineering  problems. 


We  shall  assume  that. 


p,(x)f  0 


»t  \<i<l,a<x<b. 

The  sum  +  + determines  order  of  equation  (b.l). 

We  call  the  magnitude 

a  quasiderivative  of  the  function  of  order  m,  whe,’ e 

/n  =  Vg  f  Xj  -}  V 

The  designation  of  a  quasiderivative  is  supplied  by  the  superscrint  in  braci-.e 
In  writing  out  equation  (6.1)  in  the  form 

e 


and  integrating  both  sides  of  equality  n  times  (n  =  <,  f  .  .  ^y) 

from  a  to  X,  we  obtain 


where 


a-  I 


jf  =  A(y  +  2  >'«*!  (a)F^  (.v)  +  F(x), 

t-o 


^.,-1 


•t  +  u-l 


'"-W  ! 


iinss' 


with  ttu  limi 


(6.3) 


•  •  •  I 


\  iiXMS 


(6.4) 


(x-fl)* 


■  i  J.. 


*! 

'..-1 


.)  • 


•  ■ '  a 


'a-.,-  I 


•  i 


•»a-  I 


(6.3) 


(6.6) 


Quasinormal  fundamental  functions  are  determined  by  the  series 

(ik-0.1 . n-1). 


(6.7) 


These  functions  satisfy  the  equality 

[1  l  =  *. 
0  ii^k. 


(6.8) 


The  narticular  solution  of  equation  (6.1),  satisfying  zero  initial  conditions, 

is  exoressed  by  the  series 


Y^{x)^FJcNF  +  N*F+  .  .  .  (6.9) 


As  an  examnle  we  shall  consider  equation  of  flexure  of  a  beam  of  variable 
section  on  elastic  foundation: 


^  f  (Jf)  7^1  +  *  (^).V  "  /  (^)- 

djr  [  w*. 

Here  y(x)i5  the  sag  of  axis  of  beam; 

£7(x)is  the  strengtn  of  beam  to  flexure  bend; 

l«(x)is  the  coefficient  of  elasticity  of  foundation; 

^  (x)is  the  distributed  load  oer  unit  of  length  of  beam. 

From  equation  (6.3)  we  shall  have 


^  — II  *^11* 


(6.10) 


where 


+  2y'(0)f.w+fw. 


.  .dx,. 

0  0  0  0 


;6.ii, 


Qxiasiderivatives  have  following  physical  meaning: 


y 01  (x)  =>»(-«);  >-"i(0=y”(-0; 


where  M(x)  and  Q(x)  are  the  bending  moment  and  transverse  force  in  section  x. 
The  solution  of  equation  (6.1^)  is  written  out  as: 

y(^)«=  I  y'(0)^(x)  +  K.(x),  (6.12) 

*-o 

where  the  f'onctions  Ynix)  and  K.C^f)  are  exoressed  by  the  converging  series 


(6.7)  and  (6.9). 

Equation  (6.11)  may  also  be  solved  by  the  method  of  linear  approximation. 


As  previjusly,  in  a  number  of  cases  it  is  expedient  lj  use  the  method  of  the 
mobile  orifrin  section. 

Integral  equation  (6.3)  may  be  annlicable  for  the  solution  of  nonlinear 
equations  of  the  form 


dx'i 


Pi  (^) 


dx''- 


(6.13) 


To  an  equation  of  this  tyne,  belongs  the  well-known  equation  of  M.  V. 
Ostrogradskiy 

^.+»i'  +  Py’  =  o,  (6.U) 

which  is  a  subject  of  analysis  in  works  of  a  nuEiuer  of  outstanding  mathematicians. 
Equation  (6.3)  remains  in  force,  if  only  in  equality  (6./*)  vs'e  reolace 

7(Jt)y(jf)  by  ^(Jf'y)- 

We  note  *hat  the  anplication  of  method  of  succensive  approximations  together 
with  method  of  moving  origin  gave  a  solution  of  equation  (6.14),  entirely 
satisfactory  for  engineering  aopli cations . 


C  H  A  P  T  K  R  3 


BOUNDARY  AND  NORMAL  INTEGRAL  EQUATIONS 

Modified  Fredholm  and  Volterra  integral  equations,  which  are  called  boundary 
and  normal  integral  equations  are  considered. 

Origin  of  these  terms  vd.ll  be  clear  from  the  discussion  later  on. 

It  is  oossible  to  show  that  the  boundary  integral  operator  reduces  to  a 

Fredholm  ooerator,  i.e.,  it  is  expressed  in  form 

)  G(jr,  s)y{s)ds. 
a 

and  normal  operator  is  equivalent  to  the  Volterra  operator,  but  the  presentation 
of  considered  integral  equations  in  form  of  classical  integral  equations  frequently 
is  difficult,  and  the  principal — completely  unnecessary  from  the  point  of  view  of 
practical  use. 

It  is  necessary  also  to  consider  that  the  origination  of  Fredholm  or  Volterra 
equations  frequently  dense  is  found  to  be  a  very  complicated  matter,  whereas 
boundary  and  normal  integral  equations  naturally  ersue  from  differential  equations. 

For  an  illustration  of  statement  above  it  suffices  to  turn  to  the  problem  about 
stability  of  a  rod,  for  which  construction  of  Fredholm  equation  requires  a  number 
of  artificial  reasonings  and  computations.  It  becomes  intelligible,  why  the 
classical  Fredholm  and  Volterra  integral  equations  which  proved  to  be  a  very 
effect  ive  apparat.us  for  a  general  and  qualitative  ir'vestigation,  v#ere  not  vddeLv 
used  in  solving  engineering  problems. 


Boundary  integral  equations  were  used  also  earlier  for  solution  of  theoretical 
and  apnlied  questions. 

Usually  these  equations  ensued  as  the  repult  of  anolication  of  method  of 
successive  aDproximations  during  the  solution  of  differential  equations.  However, 
with  such  an  anproach  there  was  lost  the  generality,  neculiar  to  the  apparatus  of 
integral  equations. 

As  an  examnle,  it  is  possible  to  noinl  ti;  the  meth^.d  of  successive  approximations 
in  problems  of  stability  (Vianello  method),  which  is  a  combination  of  particular 
methods  including,  graphical-analj’-tic  operatioiis. 
shall  present  another  example. 

For  calculating  a  beam  on  an  elastic  foundation  A.  N.  Krylov  used  process 
successive  anproximations,  which  can  be  represented  as  the  solution  of  a  boundary 
integral  equation 

y-'lKy+f 

by  the  method  of  simple  iteration.  Parameter  of  equation  x-— I. 

It  is  easy  to  establish  that  the  homogeneous  equation  y’^XKy 


corresponds  to  the  problem  on  the  vibration  of  a  beam  with  a  certain  distributed 

"masses''  and  has  all  positive  eigenvalue  values  Vi.  tt.. .  ,  The  process  of  simple 

!  I 


iteration  is  convergent  at 


<1. 


/ .  N.  Krylov  detected  the  divergence  of  the  process  only  by  a  direct  analysis 
of  the  obtained  series.  For  the  case  in  the  work  "On  Calculating  Beams, 

Lying  on  an  Elastic  Foundation'  the  process  of  successive  approximations,  is  proposed ; 
it  did  not  give,  however,  satisfactory  results. 

Meanvhile,  the  use  of  theory  of  integral  equations  makes  it  possible  to 
construct  effective  convergent  processes,  to  establish  a  comprehensive  generality 
between  problems  on  strength,  vibration  and  stabilit,;/  of  rods. 

Works  of  the  outstanding  scientist  A.  N.  Krylov  have  promoted  development  of 


method  of  boundary  integral  equations. 


In  article  by  P.  F.  Papkovich  ,  which  continues  the  work  of  A.  N.  Krylov  , 
there  is  indicated  the  nrocess  of  determining  the  eigenfunctions  and  eigenvalues 
now  widely  used,  somewhat  earlier  a  similar  method  was  used  ly  P.  Vetnhinkin 

The  method  of  boundary  integral  equations  in  the  works  of  E.  P.  Grossman,  D.  Yu. 
Panov,  P.  M.  Reese,  and  3.  A.  Tumarkin,  is  further  developed. 

It  is  necessary  to  note  the  works  of  Sh.  3.  Mikeladze,  in  which  there  are  widely 
used  the  Volterra  equation  and  in  individual  cases,  a  transition  to  nrirmal  i.ntegral 
equations  is  observed. 

A  consideration  of  boundary  and  standard  equations  as  a  general  mathematical 

■K-a-c-a- 

device  for  the  first  time  was  done  by  the  Soviet  scientist  Yu.  V.  Henman 

Equations,  similar  to  the  considered  equations,  Yu.  V.  Henman  called  equations 
in  indefinite  integrals". 

In  the  nresent  chanter  there  are  considered  elements  of  theory  of  boundary 
and  nor.nal  integral  equations,  there  are  indicated  methods  of  solving  homogen'^  s 
and  in  homogeneous  boundary  and  normal  integral  equations.  Much  attention  will 
riven  to  systems  of  int.egral  equations,  which  are  oresented  in  the  form,  of  matrix 
integral  equations. 

The  considered  methods  can  be  apnlied  to  any  engineering  nroblems,  which  red  ice 
to  ordinary  differential  equations  or  their  systems,  and  also  to  nartial  differential 
equatioris,  which  reduce  to  ordinary  after  a  aeoaration  nf  variables. 


*P.  F.  Pankovich,  Concerning  the  itUestion  of  Anolicati lity  of  the  Process 
luccessive  Anproximations  for  the  Flexure  of  lieams  on  an  elastic  Foundatior,,  'AnnJ.ei 
'athematics  and  Mechanics",  Vol.  1,  No.  P,  lV'33 . 

**A .  N.  Krylov,  On  Calculating  Beams,  living  on  an  Elastic  Fo'indation,  Academy/  o;' 
Iciences  of  USSR,  Moscow,  1931. 

***v.  P.  Vetchinkin,  Theory  of  Screw  Propellers,  Moscow,  Zhukovskiy,  V.  VIA  Pubi. 
louse,  1926. 

****Yu.  V.  Repman,  On  Determining  Critical  Forces  by  Equations  of  Stability, 
ransactions  of  Laboratory  of  Engineering  Mechanics,  Engineering  PuLl .  House,  Moscow, 
942. 


1.  Claasificatlon  of  Equations 


Equations  of  the  fora 

where  y  is  an  unknown  function  x;  Ky  is  a  linear  integral  operator;  is  a 

linear  functional;  ft  and  f  are  the  functions  x  (in  the  interval  a<x<b)-,  X 
is  a  oaraneter  of  the  equation,  we  shall  call  on  one-narameter  integral  equation. 
Here  and  henceforth  there  are  considered  only  real  values  of  the  functions  and  of 
the  independent  variable. 

In  a  general  case  the  operator  Ky  can  be  nresented  in  the  fora  of  a  table 


^>*=^•11  i  +  J  ]' (luty  (x,)dx^dx^+  .  .  . 

«lll 

J^,  ^  jg 

«.lt  «iit 


9  MM 

•  •  -+^#11  I  *litty  (X^ydx^-Jfqan  )  qxtt  /  +  .  .  . 

*"*  "iM  0|». 

•  •  •  Piny{Xx)dxi-\-p^.^  j  .  .  . 

(1.2) 


In  this  equality  and  P'**  are  given  functions  of  x,  Om*  and  are 

constant  munbers.  The  subscript  "r"  is  connected  with  olace  of  function  or  oara- 
meter  in  the  integral  exoressions;  the  subscript  "s"  indicates  the  line  in  tabular 
writing  (1.2);  the  subscript  'k"  is  equal  to  the  number  of  column.  After  the  line, 
containing  the  function  q,  comes  the  analogous  line  (with  the  same  number),  con¬ 
taining  function  p. 

In  a  brief  fora,  first  fora  of  nresentation  of  integral  operator  /"equality 

(1.2)  7  will  be:  -  -  ,  , 

iCy-  2  2  [quk  J  quM-  •  J  qk,ky(xt)dx„ .  .  .dxi+ 

* 

+  •  •Pk-l.$.k  J  Ptak}' {Xt)  dx„  .  .dXx). 


(1.3) 


Formulas  (1.^)  and  (1.3)  exnress  a  boundary  integral  operator.  If  all  the 


constant  limits  of  integrations  are  identical 

then  the  integr^vl  operator  (l.j)  we  agree  to  call  normal 

m,  m,  M  _  I 

1  <7i,*  ]  <7ii.  ■  •  -V*-!...*  i 

*_  I  a  a 

the/ 

In  accordance  with  this  we  shall  differentiate  /  boundary 
(equality  1.1)  and  normal  integral  equations 

i~o 


dx, 

l.U 


integral  equation 


(1.3) 


From  the  nrecedinr  it  is  clear  that  the  normal  equations  <ura  a  narticular  case 
of  the  boundary,  similarly  as  the  integral  Volterra  equations  are  a  narticular  case 
of  Fredholm  equations;  however,  essential  peculiarities  of  the  normal  equat:  -no  ma.e 
a  senarate  examination  of  them  expedient. 

Equations  (1.1)  and  (1.3)  contain  linear  functionals  L*y,  i.e.,  narameters, 
depending  on  y.  As  L*y  usually  there  are  used  values  of  the  function  y(x)  or  •.  . 
derivatives  at  certain  noints  (x  =  a,)  or  values  in  fixed  sections  of  the  integral 
expressions,  entering  into  Ky. 

The  selection  of  functionals  for  a  boundary  equation  is  not  obligatory,  since 
they  in  essence  already  are  contained  in  boundary  integral  operator.  Therefoi-e,  a; 
the  basic  form  of  the  boundary  integral  equation  it  is  possible  to  adant  the 
following : 


y=>^Ky+l. 

Solution  of  the  equation  satisfies  all  boundary  conditions  of  the  problem. 

In  rarer  cases,  the  boundary  integral  equations  are  used  in  general  form  (I.l) 
For  normal  integral  equation,  the  general  form  is  given  by  the  equality  (1.3). 

Let  us  present  an  example. 

The  differential  equation  for  the  stability  of  rod  of  variable  section, 
supported  at  ends  on  hinges,  has  the  form  (Fig.  8). 

-w+-Ti7r 

rfjr*'  EJ{J^) 


EJ{x) 


(1.7) 


Fig.  8.  Stability  of  rod. 


where  y(x)  ie  the  sag  of  axis  of  rod; 
EJ(x) — strength  of  section  of  rod  to 
f  lexiire . 

From  equation  (1.7)  we  obtain 


>(*i) 


dx,  +  ^iO). 


0 

By  repeating  the  operation  of  integration, 


we  find 


J’W  -  -  '’IJ  ^  <®> 


EJ(x^ 


(1.8) 


The  obtained  equation  is  a  normal  integral  equation  for  the  stability  of  the  rod. 


If  one  %#ere  to  determine  —(0) 

dx 


from  a  boundary  condition  >(0  then, 


we  arrive  at  the  boundary  integral  equation 

\  •  0  0  0  ' 

Operator  Ky,  entering  into  equation  (1.1),  is  linear,  i.e.,  bounded  operator, 

possessing  oroperty  of  additivity: 

^(yi+y*)=^yi+^yi.  (i*io) 

where  y^^  and  y,  are  arbitrary  integrands. 

Boundedness  (and,  consequently,  the  continuity)  of  the  ooerator  ensues  from  the 
fact  that  all  the  functions  <7frt  and  in  equality  (1.2)  are  assumed  to 

be  bounded. 


The  linear  operator  is  also  homogeneous 

^(i*y) =t*^y» 


where  ^  is  an  arbitrary  parameter. 

Also  the  functionals,  entering  into  equation  (1.1)  also  possess  analogous 
properties . 

A  boundary  or  normal  integral  equation  we  call  homogeneous,  if  it  admits  a 


t.rivi4l  solution 


y(x)  5  0. 


(1.11) 


Thus,  for  •XAflinle,  equation  (1.8)  and  (1.9)  are  homogeneous.  By  virtue  of  the 

homogeneity  of  operators  and  functionals,  entering  into  equation  (1.5)  and  (1.6), 

the  latter  will  be  homogeneous  only  in  the  case,  if 

1=0. 

In  solving  homogeneous  boundary,  equations 

y-l/Cy 

is  determined  spectrum  of  eigenvalues  {^•'  and  their  corresponding  eigenfunctions 

{yi}.  A  homogeneous  normal  equation  of  the  form 

y^kNy 

does  not  nossess  any  other  solutions,  excent  an  identity  equal  to  zero. 

In  solving  inhomogeneous  equations  (bounda.-'y  M.d  normal)  the  parameter  A  is 
gi ven. 

In  a  number  of  cases  it  is  convenient  to  use  the  second  form  of  integral 
operator : 


^y  -  +  (i.Lc^) 

whore  Q,,  and  /»,  are  given  functions  of  x, 

and  ^s'  are  constant  numbers. 

With  identical  constant  limits  of  integration  ve  shall  have  a  second  form  of  a 
normal  integral  operator: 


Ny^  1  Os]q,y('^)dx, 

a 


(1.13) 


As  a  rule,  more  simply  the  integral  equation  is  obtained  vd.th  the  first  form  of 
operator;  this  form  is  more  convenient  also  in  solving  the  integral  equation  by 
method  of  successive  approximations.  The  second  form  has  the  advantag  i  with  the 
use  of  certain  other  methods  of  solution  (for  example,  methods  of  approximation). 

In  practical  problems  there  may  be  encountered  systems  of  integral  equations, 
which  expediently  are  presented  in  matrix  form.  Thus,  a  matrix  boundary  value 
integral  equation  has  the  form 


12 


M=MA:JbH-[/). 


(1.14) 


whex*6  th«  inAtric68-coluinns  of  the  unknown  and  given  functions 


■>r 

vr 

1>1  = 

yt 

.  1/1  = 

ft 

.yo. 

and  boundary  value  operator 


(1.15) 


(1.16) 


Finally,  the  matrix  integral  equations  of  the  following  structure  are  of 
interest  : 


1 - 

: 

\K,,.  . 

a:...  . 

.  Ko, 

.  Ku 

-  y  - 
yO) 

r/.i 

/. 

.  K„  . 

.  y<’>. 

-A- 

(1.17) 


or  in  short  form 


Here  matrix-column  /”y_7  *^11  be  formed  by  function  y(x)  and  its  derivatives  uo  to 
order  inclusively. 

Equations  of  the  form  (1.17)  are  encountered,  for  example,  during  calculation 
of  shafts  for  the  critical  speed  with  a  calculation  of  the  gyroscopic  effect  of  the 
distributed  nkasses. 

♦Equations  of  this  form  may  be  called  integro-differential.  However,  their 
distinction  from  integral  equations  is  immaterial. 


7^ 


2  .  Formation  of  Integral  Equations 
from  Diffarantial  Eauationa. 


Suppose  we  ha/e  a  linear  differential  equation  of  n-th  order  with  variable 
coefficients,  given  in  the  closed  interval  a<x<b\ 

.  .  .  -i-p,(>t)y'»(A:)-/(x)  (^.i) 

with  linear  boundary  valtie  conditions  of  general  form 


if 


(a) +  (*))= 

(*-0.1.  . 


0,* 


V  ^*, 


P.*  =»  0 
(^*-0.l. 


.A-  I). 


.  .n-l). 


{'2.2) 


then  conditions  (2.2)  are  Cauchy  conditions  (at  x  -  a  there  is  given  the  value  of  the 
function  and  its  n-l  first  derivatives). 

As  the  fundamental  variable  in  composing  the  integral  equation  we  shall  take 

><*»(x)  =  T(x). 

In  considering  equalities 

y-‘>(x)  —  y<*-'>(a)-f  (x,)  dxi  et  cetera  (2.3) 

m 

we  obtain  from  equation  (2.1) 

t-  jii*'  (o)/.  W  +/ W,  *2.4) 

where 

.  .  .  j''9(x,)dx,  .  dxy, 

AW  “  —  i  Pi  (^)  •  ^2.6) 

PIquation  (2.4)  is  a  normal  integral  equation. 

In  another  form  (viz.,  in  the  form  a  Volterra  equation^  it  was  encountered 


♦  ... 
earlier  •  Squivalence  of  normal  integral  equation  and  Volterra  equation  in  a  given 

case  is  readily  established  by  means  of  a  Dirichlet  identity. 


jj  .  .  .  J  ?(x,)dx, 


•  “(TUT)! 


Owing  to  the  linearity  of  operator  Nf  the  solution  of  equation  (2.4)  in  a 


general  case  can  be  nresented  in  following  form; 


0 

vrhere  the  function  is  the  solution  of  eq'iation 

9  =A^<P  +/*, 

(*=0.1.....  n-1). 


(2.7) 


(2.8) 


and  function  ®*(-*)  satisfies  the  equation 

9=/V9+/.  (2.9) 

Suonose  }  (*=0.I . ,n — I)  is  the  sequence  of  normal  fundamental  functions  of 

equations  (2.1)  and  !'•(■*)  is  the  particular  solution  of  this  equation  at  zero 
initial  conditions.  It  is  possible  to  show  the  validity  of  the  equalities 

®,(x)  =  K<->(x). 

0,(x)=K<«>(x). 

Thus,  the  solution  of  the  normal  integral  equation  results  in  a  determination  of 
normal  fundamental  functions  of  the  corresponding  linear  differential  equation. 

If  the  initial  Cauchy  conditions  are  given,  the  function  on  right  side  of 
equation  (2.4)  is  known,  then, the  solution  of  the  normal  integral  equation  (2.4) 
determines  the  solution  of  differential  equation  (2.1),  satisfying  the  indicated 
conditions.  If  there  are  given  boundary  conditions  of  a  general  form  ^condition 
(2.2)_7,  then  by  using  equality 


♦E.  Gursa,  Course  of  Mathematical  Analysis,  Vol.  Ill,  Moscow-Leningrad.  State 
Theor  Tech  Publ.  House,  1934;  Sh.  S.  Mikeladze,  Certain  Problems  of  Structxiral 
Mechanics,  Moscow,  State  Engineer .Publ.  House,  1948. 


75' 


{b  “  n) 

><’*(6)  =»  ^  /'+•’  (<»)  Y,  + 

Iw4 

*  •»!  . 

+  H  ^  f{x,-.)dx»-.  .  .  .dx, 

(v»0, 1.  .  . ,«  —  1). 


we  shall  obtain  on  the  basis  of  condition  (2.2)  the  system  of  n-equat-ions  relative 

to  the  n-unknown 

In  solving  this  system,  we  find 

•  *  -*1 

a  a  a 


a  t  s,  Mf-i 

I  ...  i  9Md.Kf 

aa _ a 


t  tiJMa 


,  dx^ 

(»  =  0. 1.  .  .  .  .n  -  1). 


(2.11) 


The  coefficients  are  dotermined  by  coefficients,  entering  into  boundary 

conditions  (2.2). 

By  introducing,  now  t.he  relationship  (2.B)  into  equation  (2.U),  we  obtain 


+  V  f,(x)e„]  . 

*.0  a 


•'p-l 

tlMS 


.  +  (2.12) 


where 


•  —  I 


f -/(•«)+  I/,(.v)c. 

*-  0 

Equation  (2.12)  is  a  boundary  integral  equation 

9-/C9+^.  (2.13) 

where  operator  is  expressed  in  first  form  /"equality  (1.2  )_7*  The  boundary 

integral  equation  (2.13)  is  equivalent  to  the  differential  equation  (2.2)  under 
boundary  conditions  of  a  general  form. 

Vve  note  also  that  the  boundary  integral  operator  is  expressed  in  as  the  sum  of 


% 


normal  integral  operator  and  linear  series  of  functionals  and  given  functions. 

This  result  can  aleo  be  established  directly  from  equality  (1.2). 

In  practical  problems  of  boundary  conditions  frequently  they  have  a  more  simple 
structure.  Suopose,  for  example,  the  value  y(*)  of  the  problem  at 

=  (2.14) 

(*  =  0.1 . n-1). 

the  sections  a*  usually  coincide  with  ends  of  interval. 

Boundary  conditions  of  the  form  (2.14)  we  call  simple.  For  obtaining  a 
boundary  integral  equation  it  is  sufficient  in  equalities  (2.3)  to  select  each  time 
a  lovrer  limit  of  integration  in  such  a  way  that  condition  (2.14),  is  satisfied. 

For  example: 


«(*“•> (jc)  =  1  +  j  ?(•*!) 

«ii- 1 

y(.-2)  (JC)  =  _  2  + 1,  _  I  (jc  -  a,  -  2)  -1-  ^  J  ^  1  9  (-^t)  ‘^'*1- 

(2.15) 

If  conditions  (2.14)  for  certain  derivatives  are  not  given,  for  example,  for 
y<’*,  then  the  integration  is  made  by  assuming  a. —a,  and  initial  value  y'’H®) 
is  determined  similarly  to  that,  as  was  shown  for  a  general  case. 

In  certain  cases,  equation  (2.1)  is  conveniently  reduced  to  a  matrix  integral 
eqxiation  (1.17). 

In  making  the  integration  in  the  intervals  from  o.-i  t,o  x,  wo  find 

J  pi(xi)y— -  j  /»,(x»)y--»)(x,)</x,_ 

■  -  *  - 1 

-  i  ^•(-«i)>(JCi)<^^i+ >-i  +  f  /(JC,)</JC,.  (2.16) 
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Furtheraore  on«  should  consider  the  dependence 

y  0  (jf) .  jy •+•)  (X,)  rfx*  +  T/  0  -0. 1 . «-?>). 

•t 

(T,-y'>(o,).  (2.17) 

The  system  of  equations  (2.16)  and  (2.17)  can  be  presented  in  matrix  form: 


where 


y 

yi) 

y«-j) 

y*-i) 


"0 

/f..  0 

0 

—  — . 

0  . 

0  /r« 

0 

y 

yi) 

0 

0 

0  /f,-,..-, 

y*-*) 

K,-^.u 

Ka-\,  /i_| 

y.-i) 

’t. 

0 

Ti 

0 

_  • 

Y.-, 

0 

T._,_ 

•»-  1 

+ 


(2.18) 


•l-l 

.  .  ./I  -  1). 

1 

(4=0,  1,  .  .  .  ,«  —  !). 


In  engineerinf  oroblems  frequently  there  are  encountered  binomial  differential 
equations 

-<l(*)y(x)^f{x).  (2.19) 

The  fonhation  of  integral  equation  reduces  in  this  case  to  a  successive 
integration  with  proper  selection  of  constants  of  the  limits  of  integration. 


Usually  in  practical  nro'jlems 


P»  (<«)  ^  0, 

5  =  0.1,.  .  .  ,y:  a<x<fr; 

in  a  converse  case  the  coefficient  in  the  orior  derivative  equation  (2.12)  vanishes 
at  a  certain  noint  and  the  solution  must  contain  a  singular  ooint. 

Sunpose,  for  example,  there  is  given  a  differential  equation  for  the  vibration 
of  a  rod 


(2.20) 


where  y(x)  is  the  amplitude  sag  of  axis  of  rod; 

SJ(x)  is  the  strength  of  section  to  flexure; 
f  is  the  density  of  material  of  rod, 

F(x)  is  the  az*ea  of  cross  section, 

p  is  the  angular  frequency  of  the  natural  oscillations. 

Wb  now  consider  a  cantilever  rod  (Fig.  9),  for  vdiich  the  boundary  conditions 
have  the  form 


jf(o)-a  y’(0)=o,  y’>(0  =0.  —  (fyy‘’>(<))  -o. 

dx 


(2.21) 


By  integrating  both  sides  of  equality  (2.20)  from  x  to  ^  we  obtain 

/ 

^  S]  J  c-'*)  y 

s 

By  repeating  the  operation,  we  obtain 

t  i 

S  j*  j*  ^  ^ 

*  *\ 

By  transposing  EJ(x)  to  the  right  side  of  equality  and  by  integrating  twice  from 
0  to  X,  we  find  xj^,  i  i 

(2.22) 

0  0  X,x, 


-7^ 


Equ*  lion  (2.22)  is  widely  used  in  engineering  crmnutations ,  beginning  with  the 

works  of  P.  F.  Peokovieh,  E.  P.  Groesmsn  end  othere.  It  is  s  homogeneous  boundary 
integral  equation. 

There  exist  also  other  methods  of  formation  integral  equations  frt«n  differential 
equations;  they  are  reviewed  later  on  in  connection  with  applications. 


3.  The  Solution  of  Homogeneous  Boundary  Integral  Equations 


Let  us  consider  the  solution  of  the  equation 

(3.1) 

Operator  Ky  is  assumed  to  be  symmetric,  positive Lv  determined  and  the  equation 
nossesses  the  real  eigenvalues  ii,  Xi.  Ij,....  and  the  corresponding  eigenfunction- 

y>«  y*> 


The  sequence  of  the  eigenfunctions  will  form  an  orthonormal  system  i 
interval 


i  I 


(yry/)rjyj(^)yj(<)A(x)  dx 


0  i  tj. 


(3.2) 


where  h(x)  is  a  given  positive  f'lnction. 

The  indicated  sequence  is  not,  in  general,  complete ,  but  if  f(x)  is  an 
arbitrary  function  with  a  square  being  integrated,  then  on  basis  of  theorem  of  the 
Gilbert — Schmidt  function. 


(3.3) 


is  expanded  into  a  uniformly  and  absolutely  converging  Fourier  series 

^(Jc)*=  5  (3.4) 

where 

^♦“(^.yi). 

Condition  (3. 3), in  essence,  denotes  that  the  function  g(x)  may  be  an  arbitrary 
ooundaj*}'  continuous  function,  satisfying  the  boundary  conditions  of  the  problem. 

For  solution  of  equation  (3*1)  there  can  bo  used  with  necessary  changes, 
-nethods  of  solving  homogeneous  Fredholm  integral  equations.  The  most  effective  in 


practlcAl  problems  is  found  to  be  in  most  cases  the  metnod  of  successive  approxi¬ 
mations.  For  determining  the  first  (minmum)  eigenvalue  and  first  eigenfunction  the 


calculation  is  made  according  to  the  scheme 

(3.5) 

%«hex^  ^(0  and  ym  are  the  i-th  approximation  for  eigenvalue  and  eigenfunction. 

The  magnitude  X(i)  is  determineci  from  condition  of  the  very  best  "proximity"  ol 
initial  and  subsequent  approximation. 

In  equating  norms  of  the  functions 

lly(/)!l=ll>'(/-i)l'.  (3.6) 


Po-uW . 


(3.7) 


we  obtain  from  equality  (3.5) 

More  accurate  results  (for  a  given  approximation)  are  given  with  the  use  of  a  scalar 
norm  of  function 

^  “  (3.8) 

but  more  simple  are  the  calculations  neculiar  to  determining  the  norm  of  the 
function  on  basis  of  the  maximum 

11^11  =  max  Igl. 

(3.9) 

If  we  present  the  initial  apnroximation  Xto)  expanded  into  a  series  according 
to  the  eigenfunctions 

(5.10) 


then  with  the  application  of  scalar  norm  we  shall  have 

^1) 


Nj)' 


‘V  ' 


(3.11) 


At  i  Qo  X(i)  — > 


where  enproximations  give  for  ,  an  evaluation  from  above. 


t<'ith  the  use  of  norm  on  basis  of  maximum  obtain 


_ y« 

'1  XI 

•  ~~  — ,  . 

1 

r 


i-h 

*1  Xi  h  ■  ■  Cl  vi 

<1 

-1 

Xi  ’ 

[t] 

«l  XI  ' 

1  h  }  i'-'ail 

(3.12) 

In  these  equalities,  signifies  the  abscissa  of  section,  corresnonding  to 

the  maximum  Ixio)!,  is  the  same  for  the  first  approximation  y*'* 

et  cetera. 

If  function  does  not  vanish  at  one  of  the  ooints  .  et  cetera, 

then  >«i)  *♦^1  with  an  increase  of  i. 

In  practical  calculations  the  indicated  limitation  is  Immaterial,  since  the 
point  Xml  tend  to  the  point,  where  yo)  has  a  maximum  value,  and  the  first 
approximation  always  can  be  selected  in  such  a  manner  that  yi(Xa,)«^0.  This 

method  aM.v  be  called  the  method  of  cmparing  ordinates.  The  values  ho  may  be 
larger  or  smaller  than  depending  foraxample,  on  the  selection  of  the  initial 
approximation. 

During  the  calculation  it  is  convenient  to  assume 

Then  by  virtue  of  equalities  (3.6)  and  (3.7) 

Suppose,  for  example,  there  is  determined  the  frequency  of  flexural  vibrations 
of  a  rod  and  equation  (2.22)  is  used. 

In  selecting  initial  approximation  1a  the  foni 


we  obtain  with  the  anolication  of  rjori;i  on  basis  of  naximum 


*(i) 


I 


J  J j  j 

0  0  X. 


Usually  determinate  from  tnis  formula,  differs  from  the  accurate  value  by 

3  to  b%. 

With  the  use  of  the  method  of  successive  annroximatj.on3  ^equality  (3-5)_7 
eigenvalue  can  be  found  frcm  condition  of  The  minimum  of  square  deviation  with  a 
"weight"  h(x): 

«  =  JCK(i)— >'(<  !»)’ Arfxr-.  j  (X,„/v.Vo-n  -y(i-t))*A  J-v. 


From  the  condition 


we  obtain 


St 


(O 


1(0=  -  - 


(3.14) 


In  using  equality  (3.10),  we  find 


1„)  = 


■: 

•Hmr^ 

■  ■  •  (vl’ItF 

-(vM 

X(i)  -♦ 

+ 

I 

-f 


(3.15) 


giving  always  an 

evaluation  from  above. 

A  rapid  convergence,  peculiar  to  the  above-indicated  methods,  is  explained 
also  by  the  fact  that  usually  in  equalities  (3.12)  and  (3.15)  coefficient 
significant  larger  than  the  remaining. 

Let  us  turn  to  determining  the  second  and  highest  eigenvalues. 


It  is  possible  as  previously  to  proceei  froro  the  equation  (3.1)  in  solving 
it  by  the  sMthod  of  successive  approximations,  but,  as  is  known,  the  process  of 
orthogonalisation  must  be  used  both  for  the  initial,  and  also  for  ‘  no  subsequent 
approximation. 

It  is  more  convenient  to  use  the  transformed  equation 

y-y.Kiy.  (3.16) 

for  which  first  eigenvalue  is  equal  to  the  second  eigenvalue  of  equation  (3.1).  In 
the  theory  of  integral  equations  it  appears  that  the  equation  possesses  these 
nronertios  ^ 

y{x) "  ^  J  *)  ”  )  A  (5)  y  (s)  ds 


in  relation  to  equation 

In  anplication  to  boundary  value  integral  equations  we  shall  have 

.  (3.17) 

]jf]{x)k{x)dx 

m 

However,  this  operator  ensures  the  orthogonality  for  function  y^  only  of  the  initial 
approximation  and  therefore  it  is  not  useful  for  practical  calculations. 

Ws  shall  indicate  the  structure  of  operator  K^y,  which  in  solving  equation 
(3.16)  by  the  method  of  successive  aoproximations 

>(()“ 

always  assures  orthogonality  of  a  subsequent  approximation  for  the  first 

eigenfunction  independently  of  the  selection  of  initial  approximation. 

This  operator  iias  the  form 

IKyyihdx  (3.10) 

f(ty^Ky—y^  — 


T 


It  is  readily  established 


*  * 

lifiny^hdx  =  j  K^y^|-^)y^hdx'^0. 

•  a 

We  note  also  that  function  is  not  assumed  to  be  normalized. 

Equation  (3.16)  is  solved  by  the  same  methods  as  equation  (3.1). 

In  determining  the  j-th  eigenvalue  (and  eigenfunction)  there  is  solved  the 
equation 

y  = 


where 


KjV^Ky-J  y, 


(-1 


We  will  now  consider  the  solution  by  the  method  successive  approximations  of  the 
matrix  equation 

I>)=MA:i(yl. 

(3.19) 

The  eigen  "functions"  of  this  equation  will  be  designated 


>».2 


(*=».2.3 . ). 


They  satisfy  condition  of  orthogonality  (3.2)  and  normalization 
((y.I.  I>J)  =  J  y,J  (•<)  VrJ  {x)  h,  (x)^  =  {j  I 


(3.20) 


Equation  (3.19)  is  solved  on  basis  of  the  scheme 


where 


Ly(i-i).<i  _ 


.  [K][yo-u\  = 


m 

m 


In  d«t«minin|;  th«  i-th  approxiinatlon  for  the  first  eigenvalue  fron  the  condition 

II  |y(oIII-||(j'(«-.)lll. 


V  obtain 


'./-I 


In  annlying  the  norm  on  the  basis  of  maxinum  (method  of  comparing  ordinates) 


^(()“ 


(3.^1) 


or 


^(0 


(3.22) 


where  Xmi  is  the  value  x,  with  which  magnitude 
maximum  value. 

With  the  application  of  the  scalar  norm 


I  _  V 


(3.23) 


dx 


has  a 


Less  accurate  results  (for  a  given  aoproximation)  are  given  by  a  more  simole  method 

of  determining  X(.),  based  on  comnaring  one  of  the  components  lyiol.  for  examnle, 
yd)*'. 

Ws  shall  have  y(0. (^d3'(<-i).i  •  •  •  +^/»y(/-i).*) 

.  / 

and  further  - 

^2  ^'>>(1-1)./ 


(3.2i*) 


where  jc^  is  the  section,  in  which  |3'(,_n. /■!  has  a  mayi mum  value. 

It  is  Dossible  to  show  that  orocesses  of  solution  by  formulas  /”(3.22)  — 
(3.24)_7  converge  to  a  minimum  (in  absolute  value). 

For  determining  the  second  eigenvalue,  there  is  solved  the  equation 


(3.25) 


where 


In  these  equalities 


-l/Clly|-[y,l- 

1''  M  *  I* 

’  11  lyil 

'-"‘'I 

^>'11* 

I'. 

I'll 

.  l^il^ 

.yin- 

In  a  developed  form,  the  matrix  equation  (3.25)  is  equivalent  to  the  system  cf 
equations 


(3.26) 


Coefficient  ft 


is  equal  to 


II  lyil  11* 


!is( 

>-i  / 

I 

a 

Above-stated  solution  can  be  aoplied  for  calculating  the  flexiire  torsional 
vibrations  of  rods  and  in  other  problems. 


Ltt  ua  conaidtr,  now  th«  solution  of  honogeneous  integral  equations  (1.17); 

>  (0 


r  ^ 

rAfn 

.  .  Afo.  1 

><'»(x) 

•  x 

/f.. 

.  ACi. 

-  y’>(j^) . 

.  .  Af„  . 

L  y’>(jr)  J 


(3.27) 


or  more  briefly  written  as 

We  shall  call  a  scalar,  nroduce  of  an  order  of  v  functions  /  and  £  with  a  "weight' 


h,  the  following  magnitude: 


*  « 


•/-I  ‘ 


In  the  considered  case,  tne  functions  hiix),  (/-O..,.,  v)  ^  also  negative. 

We  shall  assume  that  matrix  operator  of  equation  {3.27)  is  symmetric: 


(3.29) 


For  elastic  systems  this  ensues  from  condition  of  reciprocity, 
If  operator  /”Ky_7  is  positively  determinate, 

(/))'’>>  0. 


(3.30) 


then  equation  (3.27)  possesses  real  and  positive  eigenvalues. 

Eigenfunctions  of  equation  (3.27)  yt  and  y<  by  virtue  of  conditions  (3-29) 
are  orthogonal  and  normalized 


(b/|.  [*J' 

•  r—\  ' 


Equation  (3.28)  is  solved  by  the  method  of  successive  annroximations : 

(3.31) 

Let  as  consider  first  line  of  this  equality; 

(3.32) 


where 


By  the  method  of  comparing  ordinates 


we  obtain 


>(0  =>'(/-!)  1^-^, 


1  _  >(!-») _ 


(3.33) 


The  orocess  thus  constructed  converges  to  the  least  eigenvalue  in  absolute  value. 
In  determining  the  second  eigenvalue  the  calculation  is  made  by  equation 

(3.34) 


where 


l/Q.)lI_Vl  =  [y]  _  ly^  1  1(5':’ . 


In  solving  equation  (3.34)  previously  indicated  methods  are  wised. 

4.  The  Solution  of  Homogeneous  Normal  Integral  Equations 

The  indicated  equations  have  the  form 

y  =  )./Vy  h  (^-1) 

The  number  m,  entering  into  this  equation,  corresponds  to  the  number  (of  homogeneous) 
boundary  conditions,  which  it  must  satisfy  in  the  considered  problem.  It  is  necessar. 
to  remember  that  part  of  the  boundary  conditions  (at  x  =  a)  is  satisfied  already  in 
constructing  the  operator  Ny. 

Suppose,  for  example,  there  is  considered  the  oroblem  of  flexure  vibrations  of 
a  rod  /"equation  (2.20)J7  with  boundary  conditions  (2.21).  In  taking  a  constant 
limit  of  integration  '  we  obtain  from  equality  (2.20)  by  successive  integration 

II  It 

4)  y  (-'^4)  ^•'^4  d  dXi  + 

(L2) 

+3'(0/o(Jc)+y<'»(/)/i(A).  . 

*Wb  recall  that  by  our  definition  the  equation  i?  called  homogeneous,  if  it 
admits  a  trivial  solution,  y  =  0. 


wh«r6 


AW-l,  /iM-x-l- 


(4.3) 


Th»^  ‘tr  two  conditioni  (2.21)  already  were  taken  into  account  in  the  formation 
of  equation  (4.2). 

Equation  (4.2)  is  a  normal  integral  equation. 

L,y-yiO)\  L,y^y<'H0) 

The  eolution  of  equation  (4.1)  by  virtue  of  the  linearity  of  the  integral 
ooerator  is  presented  in  the  form 

y=  lLty%(^), 

'-0  /  .  ,  N 


where  the  function 

seriea 


(x)  are  exoressed  by  absolutely  and  uniformly  converging 

(/-O.  1,  .  .  .  .  m). 


By  introducing  the  relationship  (4.4)  into  in+l  boundary  value  conditions,  ^ 
shall  obtain  bh-1  equations  relative  to  the  same  number  of  unknown  parameters 

Liy  (/  =  0 . m);  by  equating  to  zero,  the  determinant  of  the  system,  we 

find  characteristic  equation  for  determining  the  eigenvalues. 

The  indicated  scheme  in  its  basic  features  was  used  by  Sh.  E.  Mikeladze. 

We  shall  consider  a  practically  important  case,  when  equation  (4.1)  contains 
two  functionals.  This  makes  it  possible  to  formulate  also  certain  more  general 
results.  For  definitiveness  we  shall  consider  equation  (4.2). 

The  boundary  value  conditions  at  x  =  0  (Fig.  9) 

y{0)^G,  /(0)=0 


results  in  a  system  of  equations 

y(0  iu)+y<'>(0«>,(0)-0. 

and  to  the  characteristic  equation 


♦.  (0)  «,  (0) 


o'. 


(4.6) 


In  determining  the  functions  ‘t>,(jic)  /"equality  (4.5)_7  vdll  retain  terms, 
corresponding  to  the  "k-th  approximation" . 

*  k 

In  this  case  (0)  =  2  (0)  -  v  ^ 

(4.7) 


k 


* 


(0)  =  I «>('> (0)  =  2 b;K\ 


where 


uX 


(4.8) 


We  note  that  for  determining  the  derivatives  d>,(x)  special  calculations  are  not 
required,  since  their  values  already  are  contained  in  calculation  tables  for 
determining 

From  equalities  (4.6)  and  (4.7)  it  is  evident 

10,1*  if.).’  f  d,X’=0  (4.9) 

f— 0  V-^ 

or,  by  expanding  into  a  series  by  degrees  of  x,  obtain 

f»(X)=0.  (4.10) 


where  the  characteristic  polynomial 

F,(^)=  I 

•-0 


(4.11) 


1  (On-ib,  —  c,-id,)  0  <  ra  <  4 , 

/-o 

ik-m 

2  — c*-irf/i-*+i)  *  +  (4,12) 


In  a  limiting  case  at  k-*-oo  we  shall  have 

(4.13) 

Expression  (4.13)  is  a  Fredholm  series  for  the  considered  problem  on  eigenvalues. 
The  roots  of  the  equation 

fW^O  (/..U) 

Xi,  Xi.— ®11  real  and  positive  which  ensues  frcan  symmetry  and  positive 


definitiveness  of  the  corresponding  boundary  operator. 


Lat  us  &ssune,  as  usual,  to  number  these  roots  (eigenvalues)  in  increasing 
order.  Coefficients  vd.th  even  n  are  positive,  with  odd,  are  negative;  | 
From  the  equality 

the  known  relationships  follow 


et  cetera 

The  roots  of  the  characteristic  nolynomial  (/*.li)  are  approximate  values 

(‘“K2 . )  g^rid  Dossibl;/,  in  general,  both  real,  and  alsc  con.jugate  comnlex. 

Tnere  exists  an  imnortant  relationship,  valid  for  equation  with  an  arbitrary 
number  of  functionals;  the  characteristic  polynomia]  of  the  k-th  anproximation  has 
k-i  1  first  coefficients,  in  accuracy  conforming  with  corresponding  coefficients  of 
the  Fredholm  series. 

Thus,  for  examole,  at  k  2 

(/♦.16) 

but  et  cetera. 

With  the  use  of  following  degrees  of  operator  Ny,  variations  occur  only  in  the 
terms  with  n>k'\-\. 

This  circumstance  is  important  in  practical  problems,  where  it  is  usually 
required  to  determine  several  first  eigenvalues,  depending  essentiallv'  only  on 
first  coefficients  of  the  Fredholm  series. 

For  seeking  eigenvalues,  i.e.,  roots  of  equation  (4. 14),  it  is  expedient  to 

^  V 

apply  Lobacevskii  met.iod;  by  assuming  a  lack  of  multiple  roots,  we  obtain  for  a 
first  evaluation  the  approximate  equalities 


^2 


14.17) 


•nsuing  from  relationship  (4.13).  For  the  second  evaluation 


(4.18) 


Coefficients  /1<^>  are  coefficients  of  the  exriansion 

These  coefficients  are  determined  by  equality 

For  seeking  the  accurate  value  ■  there  is  required  the  use  of  the  approximation 

of  an  order  k  =  2n.  Ws  note  that  equality  (4.17)  and  (4.18)  give  in  practical 
problems  sufficient  accuracy  for  the  determination  of  the  first  two  values. 

It  is  possible  to  use  in  the  calculation  also  some  of  the  ’’inaccurate  values" 
(follo»d.nF,  directly  after  accurate  values  in  the  characteristic  nolynosiials )  and 
then  the  matter  reduces  to  seeking  the  first  roots  of  the  indicated  polyncmiala . 
Here,  there  may  be  used  different  methods,  of  which,  in  addition  to  the 

w 

Lobacevskii  method,  we  shall  give  attention  to  Newton's  simple  method  and  the 


iBSthod  of  graphic  construction  of  function  F'k(X). 


In  order  to  Judge  about  accuracy  of  the  calculation  there  should  bo  made 


analogous  calculations  for  the  polynomial 

It  is  Dossible  to  use  upper  and  lower  evaluations  for  the  first  and  second 

« 

eigenvalues  . 


5 .  The  Solution  of  Inhomogeneous  Boundary  Integral  Equations 


Let  us  consider  methods  of  solving  the  equation 

y-l/Cy  +  f.  (5.1) 

where  1  is  a  given  narameter;  f  is  a.  given  function 

Method  of  Simple  Iteration 

By  anplying  usual  nrocess  of  successive  annroximations ,  we  shall  have 


(t  e=  1,  2,  3,  .  .  .  ). 


(5.2) 


For  the  n-th  annroximation 

>^<.)-/+>A7+>*A'V+  .  .  .  (5.3) 

where  K‘/  is  the  i-th  degree  of  the  operator  Kf: 


AC/-=AC(A-7) 


(5.1) 


(/-I.  2.  3.  .  .  .  ). 

The  nrocess  of  a  simple  iteration  vd.ll  be  convergent,  if  the  sequence  of  y^m) 


converges  to  an  accurate  solution 


lim  {y(,)} 


VJe  shall  explain  adequate  conditio.. s  for  the  convergence  of  the  series 


+  W/+  ...»  5).'A'/. 

(-0 


(5.5) 


Let  us  consider  norm  of  the  integral  operator  K.  For  a  continuous  of  rator 

ii*'/ii<cii/i:. 

where 


||A/|!-max|  A/I,  ll/ll*  max  |/1, 


♦3.  A.  Bei-nstein;  Fundamentals  of  the  Dynamics  of  Structures.  Civil  Engineer. 
Publ.  House,  Moscow,  1911;  A.  F.  Smirnov,  The  Static  and  Dynamic  Stability  of 
Structures,  riailv#ay  Transport  Publ.  House,  Moscow,  1917. 


The  miniauffl  value  of  the  constant  C,  assuring  inequality  (5<6),  is  called  the  nom 
of  the  integral  operator 


Thus, 


11^/ II  <11  ^11  ll/ll. 


(5.7) 


In  accordance  with  equality  (1.2) 


||/C||<max  j|^,„lj  +|Voiil  J‘ kii*r/ km|rf-«Ct</-«f| 

I  •••  1  “lit 


+ . .  .J. 


+ 

(5.8) 


Let  us  give  an  example. 
Suppose 


/w  =  -f  • 


Then 


6  18 


|/f/l<  Xdx^dx, 


max|/(x)|<  —-X 


0<.r<t 


||lf/ll-=max|/r/|<inax 

0<«<i  0<«<i  2  I  3 


||/r||<max  - — X 

o<jt<i  2 


In  turning  to  a  general  case,  we  assume  that 

I>^MI^II  =  ^<»- 

i.e. 


IM< 


I 

\\K\\ 


In  accordance  with  this 


nax 


(5.9) 

(5.10) 


By  virtue  of  the  latter  equalities  the  majorant  series  for  y  converges,  in  which 


(5.11) 


ll/ll 

>  -IMiKh' 


» 

This  ijnnortant  result  is  a  consequence  of  Banach's  theorem  established  for 

functional  equations  vdth  linear  operators. 

Thus,  under  condition  v5.9)  the  series  converges  uniformly  and  absolutely; 

function,  being  expressed  by  the  series  (5.5),  satisfies  the  integral  equation  (5.1), 

in  which  one  may  be  convinced  by  direct  substitution. 

Further  we  shall  establish  that 

lim  =0 


for  an  arbitrary  (bounded )  function  y  )  thus,  the  solutim  does  not  daoer.d  o.i 
selection  of  the  initial  aoproximatirn . 

In  connection  with  this,  hencefore,  we  shall  assume  usually  y  ,-i>  3  0  and  then 


.h  now  turn  to  an  evaluation  of  the  error  of  successive  anrroxima’ io;  s  >'•>  • 

t.*=y— y(.«),  (5.12) 


where  v  is  the  accurate  solution  of  equation  (5.1). 

From  equality  (5.3)  at  y^o)^  0  and  exnression  (5-5)  we  obtain 

By  virtue  of  (5.9)  we  shall  have 


max  I «.  I  =  (Is  IK-. 


1  —  <7 


(•^.:i) 


It  is  of  interest  to  establish  also  another  evaluation  for  magnitude  «•. 
Ws  shall  designate  difference  between  two  successive  anfjroximations 

In  the  basis  of  equality  (5.2) 

(5.11) 


the  magnitude  is  the  error  in  satisfying  the  main  integral  equation. 


*L.  V.  Kantorovich.  Functional  Analysis  and  Aoplied  Mathematics,  'Advances 
n  Mathematical  Sciences",  No.  6,  1948. 


^<0 


From  relationshiDS  (5.12)  and  (5.14)  we  establish  that  error 


itself 


satisfies  the  integral  equation 

Now,  by  virtue  of  relationship  (5.11) 


(5.15) 


llsll  <1  _  i_,  •  (5.16) 

irtiich  establishes  connection  between  difference  of  two  successive  approxijnations 
and  the  error  of  the  solution. 

We  consider  now  the  nractically  imnortant  case,  when  the  operator  Ky  has  real 
eigenvalues  Xi(i='1.2 . <»)  and  their  corresponding  eigenfunctions  yt. 

If  error  of  the  initial  aporoximation  is  presented  in  the  form 

(5.17) 


and  there  is  considered  the  relationship 

ii 

ensuing  from  equalities  (5.2)  and  (5.12),  then  we  obtain 


(5.18) 


(5.19) 


From  the  latter  equality,  there  ensues  a  well-known  result:  the  process  of  simple 
Iteration  converges,  if 

1X|<IX»1.  (5.20) 


%diere  Xi 

equation 


smallest  (in  absolate  value)  eigenvalue  of  the  homogeneous  integral 


y^XKy. 


(5.21) 


We  note  that  in  considered  case  (the  homogeneous  equation  possesses  infinite 
spectrum  of  eigenvalues)  the  solution  of  equation  (5.1)  may  be  obtained  by  the  well 
known  expansion  into  series  by  eigenfunctions 


The  use  of  this  solution  is  very  imnortant  for  a  theoretical  analysis,  practically 
it  Is  less  effective,  than  the  aoplication  of  method  of  successive  annroximations, 
since  a  determination  of  at  least  several  first  eigenfunctions  is  required. 


Method  of  Comnlex  Iteration 
Let  us  now  consider  the  case 

(3.^3) 

Of  great  nractical  value  are  eqiations,  for  which  but  all  the  eigenvalues 

Xi  are  nositive  (flexure  of  blades  and  turb'xr.achine  disks  in  the  field  of 
centrifugal  forces,  the  flexure  nf  beams  on  elastic  foundation  et  cetera). 

Thus , 

y^-?Ky  +  /-.  Oft  (^  20 

y,^KKy,i  K>0  (<-l,  2:t..  .  .  .ck)  (j-'S) 


Process  of  complex  iteration  was  shown  by  Viarda  "Integral  equations'  )  in 
connection  with  oroblem  on  longitudinall^.’-transverse  flexure  of  rcxis .  According  ’ 
this  method 

>(/+!)  =  «>'(/) +  P(/-P^3'(.)).  /C. 


where  a  and  0  are  narameters,  identical  for  all  anproximations . 

In  accordance  with  equality  (3 •■^6)  each  subsequent  approximation  is  a  linear 
combination  of  the  two  preceding,  obtained  by  the  method  of  simple  iteration. 

Proce*'s  can  be  generalized  also  for  the  case,  when  there  is  used  a  linear 
combination  of  several  preceding  anproximations . 

The  error  of  the  i-th  anproximation 

— pKv  —  — P  if  —  {^U7) 

can  be  presented  in  the  following  form: 

Pl*^*i-i  +  (1  —  a  —  p)_y.  (3..B) 

From  this  equality  it  follows  that  ths  necessary  condition  for 

•.^0  (5.  >9) 


9H 


will  be  such: 


*  +  (5.30) 

Now  %<e  obtain 

In  presenting  the  error  of  initial  anproxiniation  in  the  form  (5.17)»  we  find 

(«  — P-^)V  (5.31) 

Equality  (5.29)  is  found  to  be  valid,  if 


or  by  taking  into  account  the  relationship  (5-30) 


-?-<l 


(5.32) 


(i«=  1,  2  3 . 00). 


It  remains  to  show  that  there  can  be  found  such  a  value  9.  that  all  inequalities 
(5.32)  will  be  satisfied.  This  condition  is  essential,  since  at  u<0  and 
such  a  value  is  impossible  to  find. 

It  is  possible  to  establish  that  inequality  (5*32)  at  will  be  satisfied. 


if 


0<p< 


(5.33) 


where  X|-  is  the  minimum  eigenvalue . 

The  value  9  expediently  is  selected  with  such  a  calculation  that  it 
corresponds  with  minimum  values  y>.  i.e.,  most  rapid  decrease  of  error.  The 
solution  of  this  question  depends  on  relationships  of  the  coefficients  Ci.  i.e.,  on 
the  character  of  error,  but  also  on  the  magnitude  .  If  it  is  assumed  that  the 
coefficient  of  expansion  of  error  by  first  form  ci  is  considerably  larger  than  the 
remaining,  then  one  should  select 


P 


^1 _ 

!*+■*, 


(5.34) 


Then  q^  =  0,  also  at  fairly  large 


(*+*! 


(5.35) 


Viarda  recomraends  the  value 


k.2a»*  (^.^6) 

with  which  for  i  ^  1  also  for  fairly  large  i 

(5.37) 

In  accordance  wiT,h  equalities  (5.^6)  and  (5.30)  the  calculation  is  made 
accordir.f  to  equation 


>'(1  +  1)— (l-p)y,,,  +  P  (/-- 


In  assuming  y(0)=0,  we  find 


(3.38) 


or. 


by  usi-.g  (y.34), 


In  solving  the  oroblem  or. 


y(i)' 


*1 


longit  udinaxly-trar. averse  flexure 


(3.;.9) 

of  rods,  we  sr.all 


;iave 

Pi  (i.U-') 

where  is  the  sag  from  effect  of  a  transverse  load;  N-is  tensile  force,  acting,  j 
the  rod;  —  critical  force  according  to  Suler.  Formula  (5.40),  possessing 
great  accuracy,  is  used  in  engineering  calculations.  It  could  have  beer  obtained 
also  frcom  equality  (5.-^^),  if  it  is  assumed  that  functions  "[  ar.d  agree  with  an 
accuracy  un  to  factor  (cul've  of  flexure  from  lateral  load  is  similar  to  the  f^rst 
form  for  the  loss  of  stability). 

Process  of  complex  iteration  for  the  equation  (5.2.i)  is  convergent  at  any 
values  ji  but  at  u'  Xi  the  convergence  is  obtained  more  gradually  since 
is  close  to  unity  /':See  (5.35  and  (5.37)^. 


Method  of  Complex  Iteration  With  a  Variable  Parameter 
A  gradual  convergence  with  large  tn  in  the  preceding  method  is  associated 
•/ith  the  fact  that  parameter  fl  was  taken  as  constant  for  all  approximations, 
rocess  of  coonlex  iteration  /"equality  (5.38)_7  can  be  written  out  in  the  following 


orm: 


>(•+!)  ■=  >(0 + PAy<*>* 


(5.41) 


where 


(5.42) 


Magnitude  ^y(t}  is  the  difference  between  two  successive  approxinations  or  error 
in  satisfying  function  ^(0  of  equation  (5.24). 

Process  of  ccmplex  iteration  with  variable  narameter  is  expressed  by  the 


equality 


(i-a  1. 2. .  .  .). 


(5.43) 


where  the  parameter  can  be  determined  from  the  condition,  so  that  the  function 
>(!+»  in  the  very  best  manner  satisfies  the  main  integral  of  equation  (5.24). 

For  generality  we  shall  return  again  to  equation  (5.1)  and  then 

^>(0  “  /  +  o  ~  !'( <  )• 

From  the  condition  of  minimum  of  the  square  deviation  with  the  "weight"  h(x) 
/~the  function  h(x)  enters  into  the  condition  of  orthogonality  of  the  eigenfunctions/. 

(5.44) 

vrtiere 


we  obtain 


(5.45) 

(5.46) 


From  equalities  (5.45)  and  (5.43)  we  find 

u  =  ^>(0  ~  Pi  (^y(n  —  ^^y(i>) 

and  by  virtue  of  (5.46) 

P  —J! _ 

S 

It  is  possible  to  establish  that 


(5.47) 


(5.48) 


a>ei 

‘J?*! 


<0 


loj 


and  thar«foi*a  conditions  (5 >46)  determines  minimimi  of  the  error. 

We  note  that  in  determining  it  is  found  necessary  to  calculate  K&ytOt 
which  insediately  is  used  in  the  following  approximation: 

J ) = >(<  f  I )  +  Ph- I  ^  1  >• 

The  magnitude  &y(i+o  is  determined  from  equality  (5.47).  The  proof  of  the 
convergence  of  this  method,  as  also  of  several  subsequent  methods  is  made  difficult, 
however,  it  is  obvious  that  if  orocess  of  the  solution  converges,  then  it  converges 
to  an  accurate  solution. 

The  latter  immediately  ensues  from  equality  (5.43),  since,  if 

>(o|  <8, 

then 

I  *>'<'>  I  <T^. 


and  the  function  >(i)  satisfies  the  integral  equation  with  an  error  not  exceeding 

> 

IM  ‘ 

Ws  shall  show  also  that  if  function  on  the  right  side  ^  can  be  expressed  in 

form 

/-c*y». 

then  already  the  first  approxijnation  on  basis  of  equality  (5.43)  results  in  an 
accurate  solution  (initial  approximation  is,  as  usual  y«t')  =  0)- 


Ws  shall 


Pi» 


This  result  ensues  also  from  equality  (5.22). 

Ws  note  that  it  is  valid  with  the  arbitrary  function  h.  In  oractical 
calculations  (for  simplicity)  it  is  possible  to  assume 

4(x)-l. 


In  oroblem  about  the  flexure  of  turbomachine  blades,  propellers,  the  indicated 

|A|  ^ 

process  gives  entirely  satisfactory  results  at 

Method  of  Similar  Iteration 

Let  us  assume  that  an  approximation  can  be  improved  by  multiplying  be  a 

certain  coefficient  can  be  assumed  as  the  initial  appr  ^ximation 

>1.)  =  ?<!'(. )  (5./»9) 

The  subsequent  approximation  is  determined  from  "be  equality 

y(i+i)  =  /  +  P/^>(0  (5.50) 

{1  =  0,  1.  2.  3,  .  .  .  )• 

The  coefficient  p,  can  be  determined  from  different  considerations.  In  asstaning 
that  the  initial  and  subsequent  approximations  coincide  at  the  point  (a 
where  yo)  has  a  maximum  value,  we  obtain 

>(/+!)  =  (5.51) 

whence  x 

B,  B  -  . 

Relationship  (5.51)  is  equivalent  to  the  condition  of  equality  to  zero  of  the 
error  of  integral  equation  at  x  =  x,^,  if  into  this  equation  is  introduced 

"an  improved"  approximation: 

The  process  of  successive  approximation  constructed  in  this  manner  converges 

A  I  . 

usually  at  "xi"  [  ^  ^  ^  similar  method  was  used  by  V.  P.  Vetchinkin. 

The  significantly  best  results  are  obtained  in  the  case,  if  is  determined 
from  the  condition  ,  *  /  *  \ 

♦In  this  section  of  instructions  on  convergence  of  processes  are  given  on  the 
basis  of  an  experiment  of  their  application  during  calculation  of  turbomachine 
blades . 


which  gives 


5/^ 

4 


(5.52) 


The  process  of  solution,  given  by  equalities  (5*50)  and  (5.52),  we  shall  call 
a  sinilAT  iteration  by  the  equality  of  areas.  Practically  this  process  converges 

<5. 


<  10,  where  it  is  estjecially  effective  at 
The  ranidity  of  convergence  increases,  if  as  an  anproximation  of  (*  >  3) 
there  is  taken  the  half  sum  of  two  nreceding  anproxixaations : 

Finally,  one  can  determine  and  from  the  condition  of  the  minimum  of 

the  square  deviation  of  initial  and  subsequent  approximation 

^<  -](>{<+  »  -^Ji,)*  (  or  ^4  “  I  ^ • 

i.e.  for  equality 

which  gives 


ft- 


}O(0- 


(5.53) 


One  of  the  Important  variants  of  the  method  of  similar  iteration  (on  basis  of 


equality  of  functions)  ootnted  out  by  S.  A.  Tumarkin  .  Let  the  functions  >'«'+»  and 

Th; 

(5.5A) 


coincide  in  all  sections.  This  is  possible  at 

/ 


y<*>  conforsis  with  the  accurate  solution,  then  the  parameter  «,  will  be 

constant;  in  reality,  equality  (5.54)  determines  the  magnitude,  depending  on  x. 

In  using  in  an  approximate  solution,  equality  (5.50)  established  for 
we  shall  obtain 


•This  method  was  comnunicated  to  author  in  a  Tjersonal  visit. 


(5.55) 


it 


>(.) 


=  0.  1.2,.  .  .) 


The  nrocess  of  successive  epproximations  converges  in  practical  problems  at 

Relationship  (5.55)  loses  meaning,  when  the  denominator  of  the  expression 
tends  to  zero  %dtich  virtually  is  encountered  fairly  rarely. 

We  now  shall  point  out  that  for  all  variants  of  method  of  sisdlar  iteration 
the  following  is  valid: 


A.  If  a  process  of  subsequent  approximations  cotjverges,  then  converges  to  an 
accurate  solution: 

B.  If  function  ^  is  similar  to  one  of  eigenfunctions,  then  already  first 
aonroximation  results  in  an  accurate  solution. 

Suopose,  therefore. 


>('+!)-♦>(/)• 


(5.56) 


In  order  tliat  tends  toward  an  accurate  solution,  there  siust  be  in 

accordance  with  equality  (5.50) 

(5.57) 

Sunnose,  for  example,  the  process  of  similar  iteration  is  made  on  the  basis  of 
equality  of  areas:  • 

I/djt 

><*+»  ■“  ^yo)  i — - - + /. 

(5.58) 


By  integrating  both  sides  of  equality  from  a  to  b,  we  obtain 

5/d* 


f  .  f  h 

J  J 


\y^^.^\)4x-*  ]  >(()  dx. 


lo6 


By  virtue  of  (5.56) 


whence  also  there  ensues  the  reletionshin  (5.57). 

If  the  solution  is  constructed  on  basis  of  equation  (5.55),  then  from  (5.56) 

there  ensues 


1  .e . , 


the  function  Lends  toward  an  accurate  solution  of  the  problem. 
We  shall  now  indicate  the  validity  of  the  secotvl  assertion. 

Sunnose 


where  is  the  n-th  eigenfunction  of  onerator  Ky. 

Then,  by  assuming,  as  usual,  y(0)*»/,  we  obtain  from  equality  (5.5B) 


;  y<i 


yiKdx 


_/ 

I- 


An  analogous  result  is  obtained  from  equality  (>.55). 

I'M  note  that  "quality''  of  the  aoproximation  can  be  evaluated  also  by  the 
magnitude  0i.  With  a  good  convergence  of  the  nrocess  already  for  second  or 
third  annroximation 

'We  shall  dwell  briefly  on  the  solution  of  matrix  equations 

'We  shall  assume  that  the  corresnonding  homogeneous  equation  has  the  eigenvalues 
Xi,  X|,  Xj,... 


The  nrocess  cf  simple  iteration  for  equation  (5.59) 

l,V|HI)l  =  ^  +/ 

convorpes,  if 

ixKU.l 

At  |X|>|  Xi|  there  are  constructed  nrocesses  of  iteration,  analogous  to  the 

nreviouslv  Indcated  nrocesses. 

Wb  shall  consider  as  example  the  application  of  the  method  of  similar  iteration. 

In  thl,  cast  U,H..|-f,AlAriljr,„|  +  |/|. 


loG 


Th#  HAgnituds  Is  dstsmined  from  following  condition:  ths  corrsctsd 


aDproxiostion 


Kol-My(ol 


in  ths  vsry  best  mannsr  has  to  satisfy  equation  (5.59). 

The  error  of  equation  (5.59) 

U>Ji)I  *  P/  l^l(>(ol  +  (/I  —Pi  (5.60) 

Thus, 

'  'P/  y  +  /i  —  Pil'c).  1 

Pil  2  1-^  f\  —  P/l'jO.  J 

/—I  I . 


+  /«— Piy(<).« 


The  absolute  value  of  the  error 

By  detersdning  J,  frosi  the  condition  of  minimum  of  squai*e  deviation  J|  (Ay*,,]  l*dx, 

Oti  • 

i.e.  from  the  equality  obtain 

Pi  “  srr?  7. - s —  • 

2  2^./>’(0.y)  ** 

»— I « '  •  y— I  / 


If  we  take 


Pll  3^(0. 1 
P<ay(0.  ■ 


then  for  the  method  of  similar  iteration  on  the  basis  of  the  equality  of  functions 

we  obtain  the  system  of  equations  (at  /.4  0) 

Aj'io. , 


J'(i)..  ^^£,^.y>(/).y 

(''*=1.2 . «). 


Ic7 


(5.61) 


There  exist  also  different  methods  of  anproximation,  reducing  the  nroblom  to 
solution  of  a  system  of  linear  algebraic  equations  (Fredholm  method,  method  of 
ffloaents ,  method  of  collocation  and  others).  In  nractical  problems,  these  methods 
usually  are  inferior  to  the  nreviously  indicated  metnods . 


0.  The  Solution  of  Inhomogeneous  Normal  Integral  Fquations 


In  thie  solution  of  equation  (1.5)  its  comnonent  funct'onals  are  consLderdd  as 


naramoters  and  for  brevity  are  designated 


(6.1) 


In  solving  a  normal  equation  by  the  methcxi  of  successive  aoproxiraations 

+  X  C/j 4  /  (b..2) 

we  obtain 

2  O,  (b.3) 

where 

(6.4) 


.  .  . . 

0=/+)iV/4- W/+  . 


(6.5) 


The  narameters  Cj  are  detennined  from  boundary  conditions  of  the  oroblem.  oeries 
(6.4)  and  (6.5)  converge  uniformly  and  absolutely  with  an  arbitrary  value  1.  Let 
us  note  that,  if  calculation  is  made  on  basis  of  equation  (6.2)  and  the  narameters 
are  determined  after  each  annroximation,  then  nrocess  converges  only  at 

ixKU.i. 

In  "decinhering''  the  Ci  values  in  accordance  with  equality  (6.1)  the  normal 
equation  acquires  all  the  orone'-ties  of  a  boundary  equation. 

Sometimes  it  is  convenient  for  the  series  (6.4)  and  (6.5)  to  be  used  in  another 
form,  for  example, 

-h  ...  (L  t\ 


where 


=  (s  =  l,  2,  3.  .  .  .).  = 


In  above  mentioned  series  each  term  of  the  series  is  the  difference  between 
two  successive'  approxiaations .  The  calculations  are  stopped,  when  a  new  term  of 
the  series  is  small  in  comparison  with  sum  of  all  the  obtained  terms  of  the  series. 
In  evaluating  the  convergence  of  the  process  of  successive  approximations  it  can  be 
established  that  the  terms  entering  into  the  series  (6.4)-(6.6)  diminish  with  an 
accuracy  up  to  a  factor,  as  ^)1.  _  large  values  of  the  parameter  X 

ft . 

and  of  the  interval  of  determining  the  function,  the  convergence  deteriorates  and 
in  a  number  of  practical  problems  even  at  becomes  gradual.  For 

improving  the  convergence,  it  is  possible  to  aoolj'  method  of "mobile  origin",  with 
which  the  function  is  found  at  first  in  the  fairlj'  small  section  a<x<ai,  which 
assures  a  rapid  convergence.  Fiu'ther,  the  solution  is  constructed  for  the 
follo%cing  section  (in  normal  integral  equation  it  is  assumed  and  initial 

parameters  in  section  Jf=Oi  are  determined  from  the  preceding  solution.  The 
condition  of  "rapid  convergence"  of  the  simple  iteration  process  for  a  normal 
equation  can  be  recooiasnded  in  such  a  form: 

iXlllAniO-  (6.7) 

In  a  number  of  problems  in  structural  mechanics  under  condition  (6.7)  there  is 
required  not  more  than  three  to  four  approximations. 

Analogous  results  aire  obtained  in  solving  matr^  equations 

lyl^Xl/VIIyH- |c,|/,l  +  (/l 

b2^  the  method  of  successive  approximations;  we  shall 

(>)-  2C,(0,1  +  (4»I. 

.  .  . 

the  fur.ctions  (Oj)  and  (®J  can  be  presented  in  the  form  of  a  series,  for 
example ,  _ 

I®I- 


lo^ 


in  which 

I®,o)l  -  I/l. 

Ml  shall  now  consider  ths  application  of  the  linear  approximation  method,  in 
which  there  is  used  an  aoproximate  integration  by  the  trapesoidal  rule. 

The  indicated  nsthod  for  boundary  integral  equations  results  in  the  necessity 
of  solving  a  system  of  algebraic  eqxiatlons;  for  normal  integral  equations  the  matter 
is  greatly  sisplifled  because  the  values  y(*i)  (i^O,  1,  2,...)  can  be 

determined  successively  one  after  another. 

Suppose  the  normal  equation  is  given  in  form 

yW-  £Q,WW,(x,)y(x,)rfx,+/+  ICJ„  (6.8) 

•—la 

/'parameter  x  is  contained  in  the  function  Q<(x)]. 

Solution  of  equation  (6.8)  will  have  form  (6.3),  and,  for  example,  for 
determining  function  ^(x)  there  must  be  solved  the  equation 

IQ,  (JC)  !</,  (Xi)y  (Jf.)  dx,  +/.  (6.9) 

•—I  a 

Now  'im  consider  first  variant  of  method  of  linear  approximation.  Mb  shall 
subdivide  the  interval  of  x  variation  into  k  sectors  »#ith  the  sections 

The  length  of  a  sector 

X/-|»Ay. 

we  shall  designate 

Q*(Xy)-<?,y.  /(Xy)-/^ 

I,  2, .  .  . ,  k, 

s li  •  .  t  e. 

On  the  basis  of  equality  (6.9)  we  shall  establish 

,  •  r 

t/-l.  Z. 


J:lo 


Hance , 


#w| 


/-I 


^  q,  ,y.  (A.  +  A.  -i) 


+  //  • 


In  final  form 


where 


yi= 


1 

1  — d/y 


(/-I.  2.  .  .  .  ,*). 


.  (a,  +  'i.+  l) 

«— I 

(/“I.  2..  .  ..*;  v=l,  2, 


« 

*>>  **4~  yi 

•-1 

«»1 


y-1). 


(6.11) 


Calculation  by  equation  (6.11)  is  convenient  to  make  according  to  scheme,  shown  in 
Table  5.  Magnitude  y«“/*  is  multiplied  by  the  column  "0"  of  the  table 
and  the  result  is  filled  in  column  "0"  of  table  5y  summarizing  the  numbers 

in  first  line  and  by  dividing  by  1  —an.  we  find  y^.  By  multiplying  y^  by  colunn 

"1"  of  table  a/,  et  cetera,  we  find  all  the  values  yi- 

A  simplification  of  the  calculation  is  attained  by  subdividing  the  section 
into  equal  intervals. 

With  equal  intervals,  in  addition  to  a  linear  approximation  (integrating  by 
trapezodial  rule),  there  may  be  used  integration  by  Simpson's  rule  et  cetera  which 


1:11 


•/u»  ixtfo  I  •?» 


results  in  the  epnearance  of  corresnondinf  factors  for  the  coefficients  3)*' 

The  internolation  formulas  by  Chebysh.ev  and  Gauss  are  not  successfully  annlied  ir. 
a  similar  manner,  since  1’ )r  different  sectors  t:.e  ooints  of  internolation  will  not 
be  common. 

iVe  now  discuss  the  second  variant  ol'  the  linear  annroximatlon  which  leads  tc 
another  scheme  of  calculation.  V/e  snail  write  equality  (6.10)  in  such  a  form: 


=  .  !  -■  i  j,. 

#*•1 


(6.12) 


where 


A. 


t 


0-1.  2.  .  .  .  k~\). 

(-'*.0=  Y  ‘7icyo-^i  )• 


(0.13) 


From  relationship  (6.12)  we  shall  obtain  the  calculation,  formula 


yi 


(y=i. 


o 


.  *):  CVo  =  /o). 


9% 

$-\ 


(h.U) 


Thus , 


/ 


With  a  snail  numbsr  of  calculation  ssctions  ths  mors  convsnient  Is  the 


first  variant  of  linear  approxination  nsthod,  with  a  large  number  of  sections, 
the  second  variant. 

W<s  now  consider  the  application  of  tho  linear  approximation  method  for  the 
solution  of  matrix  normal  Integral  equations 

bJ-M/vilj-I+l/l. 


As  an  example  we  take  the  equation 


The  mr;^’nitude  of  the  known  paraxneter  X  is  included  in  the  operators. 

If  normal  operators  are  given  in  the  second  form  /"equation  (6.8)_7>  then  in 
an  exnanded  writrng  we  shall  have 

>1  (x)  -=  i  Qi,  (^)  1  ^1,  (Jfi) yi  (^i)  dx^  4- 

+  ITu  (x)  f  tu  (^i)  yt  (-^i) 

yt  (•<)  -  2  Qt,  (•*)  f  (  v,) j'l  (  t,)  rf.t,  4- 

**■>  *. 

+  2  W  I  Myt  (-^i)  dx^  4-/,. 

m 


As  earlier,  we  shall  briefly  designate 

yi(^j)"*yijl  Qi,  “» Qi,/  ^1,  (■*/)=  ^1,/ •  •  •  et  cetera 

In  using  the  second  variant  of  the  method,  similar  to  relationship  (6.12)  we  obtain 

•—I  ^ 

•-I 

a. 

yv  “  2  1  + 

I 

+ 2  ^  ^tijytAi^ + /i/. 


(6.16) 


•-I 


IIH 


wh*r« 


"^15/  1 + Y  i^j + A/+i) 


(y-i.  2..  .  ..ik-i), 


and  t.  p. 


Equality  (6.16)  vrill  be  written  as: 


where 


^ujyij + *nyy  V 
Hijyx!  ^ttjytj  “  ^tjt 


1 

*ii;  =  ~  ~  ^ 

« 

®ti;  “  — ^  2  Qttfluj' 


*My=>  “ 


«-i 


Fw-/./  +  io,/,';.’/-,  +  ,-,• 

F„-/„+'iQuii'J!,-,  +  s',  Vi!!/-.- 


(fc.l7) 


From  equations  (6.17)  we  obtain  the  calculation  formulas 

yi/“ - I 

*iiy"*y  —  *iiy**iy 

•  (/=  1,  2,  .  .  ,*). 

As  nrevicusly,  the  values  y,(x)  and  y»(*)  in  each  calculation  are  determined 
in  sequence. 


CHAPTER  U 


APPUCATION  OF  BOUNDARY  AND  NORMAL  INTEGRAL  EQUATIONS 
TO  PROBLEMS  OF  STRUCTURAL  MECHANICS 

1.  Fl»xur*  of  Rods  In  4  Fltld  of  Ctntrifugal  Forces 
This  problsB  hss  s  nuabsr  of  tnglntsrlng  sppUcmtions  in  calculating  for  the 
strength  of  blades  of  stssD  and  gas  turbines,  ccmpressors,  blades  of  propellers  and 
helicopters.  The  application  of  method  of  successive  approximations  and  boundary 
Integral  equations  In  calculating  for  the  strength  of  propellers  Is  given  In  works 
by  V.  P.  Vetchinkln,  D.  Yu.  Panov,  P.  M.  Rlsa,  S.  A.  Tumarkln. 


Fig.  10.  Flexure  of  rod  In  a  field  of  centrifugal  forces. 

In  Fig.  10  there  Is  shown  the  system  of  coordinates  being  used.  The  y  axis 
coincides  with  axis  of  rotation,  the  f  axis  Is  directed  radially  and  passes 


through  the  center  of  gravity  of  the  root  section  of  rod.  Axis  of  rod  is  assumed 
in  the  form  of  space  curve,  deviating  little  from  a  radial  direction. 


(1.1) 


The  origin  of  the  local  axes  x^y^  is  placed  at  the  center  of  gravity  of  cross 
section;  principal  axes  of  the  section  I.  are  turned  by  an  angle  a  in 
reference  to  the  local  system. 

Equations  of  the  flexure  of  a  naturally  twisted  rod  have  the  form: 

\bJ^  B^^|  ’  3  \ej,  EJf  I  ‘ 

Af.,  +  -Lf  '—.- L\sln2.A)„. 

^  \£y,  f/j/  2  U./,  £Ji)  ^ 

where  u  and  v  are  the  elastic  disnlacements  of  axis  of  rod  along  the  x  and  y  axes 

respectively;  BJ%  and  Eh  — are  the  orimary  rigidities  of  section  during 

flexure.  The  bending  moments  Al,i  and  during  flexure  of  rod  in  a  fi^Ad  of 

f: 

centrifugal  forces  are  equal  to  ^ 

Af,i  (r)  =  (r)  +po)»J  (v  (rj  - 1»  (r))  (r,)  dr^, 

M,X  W  +(»•  lV.“  w-™  W  )/'(<•.) 

where  p-  is  the  density  of  the  material;  «>  is  the  angular  velocity;  F  —  across 
section  area.  .  Mj,i{r)  and  designate  bending  moments  from  a  transverse 

load  Pm  and  Py  and  the  initial  distortion  of  axis  /”its  coordinates  orior  to 
deformation  >’•('’^1 

(r)^  -  II  Py  (''*)  +  P^>’ J  CV*  (r,)  -y,  (r))  r^F  (r,)  dr^.  ^  ^  ^  ^ 

iM-i  (0  •=  f  IPm  ('»)  i  (O  -  '’•^0  ('•.))  f  (/'»)  dr^. 

"  r  fi  ' 

Equation  (1.1)  which  takes  into  account  relationships  (1.2)  vdll  furn>  a  system  of 
integro-different.ial  equations.  We  shall  transform  it  into  a  matrix  boundary 
integral  equation. 

As  the  chief  unknowns  we  shall  take  the  comnonent:  ;he  curvature 


dr* 


ii.U 


"^3.  A.  Tumarkin,  Equilibrium  and  Oscillation  of  Twisted  Rods,  Transactions  of 
the  Central  Aero-Hydrodynamic  Institute,  No.  3/a,  193'a 


Qjr  MAns  of  IntogrAtion  bgr  ports  froM  rolttlonohlpi  (1.2)  m  find 

(15) 

t 

f 

C%  I  J 

C(r)  !•  tht  etntrlfucAl  fcreti  oetinf  in  tho  tootlon  r.  For  «  rigidly  flxtd  blade 
in  tho  root  ooetion  , 


whort 


5-1 


(1.6) 


By  introducing  equality  (1.5)  with  a  oonoideration  of  the  dependence  (1.6)  in 
equation  (1.1),  wo  obtain 


ffl  ,  Ifairt 
1 1 J  k..  Ktn\[  t 


f 

♦ 


(1.7) 


whoro  tho  paraoMter  of  equation 


X-— P. 

and  integral  operators  are  detersdned  by  the  equalities 


(1.8) 


-y  (s —  •*"*“  f  I 

'  i  u 

M  f, 

^•if-y(-^~;^J*in2e.rj'c,  (rj)  ^  J  r,? (r,) dr, dr, ;  (1.9) 


The  Magnitudes  f  and  ^  signify  the  conponents  of  an  elastic  curvature 
cAuMd  by  tht  momtnts  And  M^x  ;  thty  art  obtaintd  from  tquAlities  (1.1) 


tfitr  rtplkclng  ^mx  and  M^x  by  aI^,  and  calculation  by 

aquation  (1.7)  la  ainpla,  since  it  contains  only  two  integral  operators  (%dth  an 
accuracy  up  to  a  factor).  In  a  number  of  practical  problems  (for  example ,  in 
calculating  the  blades  of  turbomachines)  it  is  possible  to  ignore  the  flexure  in 
the  plane  of  maximum  rigidity,  since  A>-^v 

In  ^his  case,  from  equations  (1.1)  v/e  obtain 


where 


At. 

—  e=  -A  COS  a  ; 

dt* 

iPv 

Af- 

■=  — 5.  sin  a, 

dr' 

cos  «  —  Mj,x  sin  » 


(1.10) 


is  the  bending  moment  relative  to  axis  with  a  minimal  moment  of  inertia. 


From  equalities  (1.11)  and  (1.5)  there  ensues 

K 

Af, » iif,  -  sin  « (r)j  (''i)  ^  (''i) 

r 

-.•cos.(r)r  Jc,(r,)  A 


(1.12) 


In  considering  now  relationships  (1.6)  and  (1.10),  we  shall  obtain  a  boundary 
integral  equation  relative  to  the  bending  moment  • 

Mr„  (  1  •  13  ) 

where 

R 


X= 


/fAf,  —  sin  a  (r) 


Mi: 


(rt) 


- - Sin  a  (r,)  dr.  dti  + 

(''f) 


R  r, 


+  C08a 


J  J 


cos  a  (r,)  dr^dr^. 


(1.14) 


For  a  blade,  secured  to  cylindrical  hinge,  axis  of  which  coincides  with  axis  vj 
for  the  root  section,  the  integral  equation  will  be  such: 


(1.15) 


1  I  7 


(1.16) 


\ 


whare  Vw  —  •*;  . 


r 


where  KM^  and  >Mi|  have  th^  former  meaning,  and  ^ 

B  (r)  (aj)  sin  a  (r)  j*  C^  (r,)  rfr,  +  cos  a  (r)  rr^  J 

'  f 

A  further  simplification  can  be  attained,  if  it  is  assumed 

tan  9"^ar 

■» 

(vane  or  blade  of  constant  screw  nitch  ). 

Coefficient  a  may  be  selected  equal  to: 

ten  ac 

a  =  - 


»C11 


’Zf 


(1.1?) 


where  r*,  ^  is  the  average  radius;  «cp  =  ®(^p). 

The  calculations  showed  that  replacement  of  real  angle  of  ir.sta.Ua l  \  «n  of 


profile  in  vane  of  a  turboroachine  by  an  angle,  determinable  from  equality  (l.i'^j, 
does  not  introduce  a  noticeable  error. 

Wa  shall  have 

lift  dr* 


By  integrating  both  sides  of  equality  from  >■»  to  r  ,  we  obtain  for  a  rigidly 
fixed  blade 

The  latter  expression  is  valid  also  for  a  blade  fastened  on  hinges  by  virtue  of  the 
equality 

.J  (rO  =^(»’,)^"»“(»’o). 

f 

Now  frop  relationships  (1.12)  and  (1.1?)  it  follows 
,  )  M,  (r)  -  M,  (r)  -  u.«  sin  •  (0  J  (l  +  -^)  C,  (rO  ^  (r.)  dr, « 

t 

—  (r)— •••ln«(r) 

*0.  Yu.  Panov,  Calculation  of  a  Propeller  for  Strength,  Transactions  of  Central 
Aero-Hydrodynamic  institute.  No.  288,  1937. 
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J  ito*<i(r,)  dr  '  ”  * 


(1.18) 


In  vi«w  of  the  dependence 

^9  ' 


7-JtWrfr, +&(,,).  (1.19) 

'•  r% 

relatlonehip  (1.18)  can  be  presented  in  the  fom  of  a  boundary  Integral  equation. 
By  introducing  equality  (1.19)  for  rigidly  fixed  blade  in  i*elation- 

ship  (1.18),  we  obtain 

n 

Af,  (r)  *  Af,  (r)  —  m*  sin  o  (r) 

j  fin-«(r|| 


f  -£i-(OL  ^ sin  «  (r,) dr^dr^.  (1.20) 


Integral  equation  for  the  curvature  'HO  is  obtained  from  relationship  (1.18), 

if  both  sides  of  the  equality  are  divided  by  £/,  — —  and  dependence  (1.19)  i3 

•in* 


taken  into  account: 


(0  "j*  (0  •"  ”*  — —  f  — -  dr^dti.  (1.21) 

^  m  m. 


Integral  equation  relative  to  the  angle  of  rotation  has  the  form 

(1.22) 

dr  J  £J  J  sln’a(r,)  dr  ' 

'.  ’  rt 

Equations  (1.20-3.22)  correspond  to  one  and  the  same  problem,  however,  they  have 
their  own  peculiarities  from  the  point  of  view  of  practical  use.  The  difference  is 
found  to  be  also  the  norms  of  operators,  entering  into  these  equations. 

Influence  of  the  centrifugal  forces  on  the  flexure  depends  on  the  dimension¬ 
less  parameter  of  the  flexibility  of  rod 

»* 

where  «u*,=  X,  — first  eigenvalue  of  the  homogeneous  equation  (•'r  is 

the  value  of  the  angular  velocity  of  rotation,  with  which  centrifugal  force?, 
redirected  to  compression,  cause  a  loss  of  stability  of  the  vane). 

'.Vith  small  /alues  of  the  parameter  of  flexibility  0  the  effect  of 

elastic  dtiformai'' '■•'.0  the  rod  can  be  ignored  by  assuning 


If  the  bending  moment  Mn  from  a  lateral  load  (and  initial  displacements  of  axis 


* 


\ 

I 


whtr«  M^tii  ttia  distribution  of  ths  bonding  momants  with  the  first  form  of  tho 


rod's  loss  of  stsblllty;  o  Is  s  eosfflclsnt,  thsni 


l.e.,  with  Isrgo  flexibility  paramstarsj  the  effect  of  elastic  deformations  in  a 
field  of  centrifugal  forces  may  be  very  great. 

The  method  of  solving  the  boundary  integral  equations  was  diacusised  earlier. 

Ws  present  results  of  the  calculation,  relating  to  Lho  pri.blom  on  tha  flexure 
of  a  rod  of  uniform  eection  under  the  effect  of  diatributod  transverse  and  axiai 
loads  of  constant  intensity  (Table  6). 

Table  6.  Values  for  a  Rod  with  a  Flexibility  raramater  ^  ?.'3H 


(al  Meraji  fwuitHMi 

- 

•Foptiyia 

i^t)  lIpHOAinKrHHe  |j)  Tomiof 

1 

2 

3 

peuieHiie 

(4) 

CjOWHIII  NTtpiUHR  no  BHipXI 

(S  M) 

0.436 

0,300 

0.426 

0,412 

fV> 

CmwHIII  NTCpaUHl  C  nepeMtHHMM 
napfMTpoM 

(S.43) 

0,280 

0,399 

0,409 

0.412 

neaoSNai  NTtpiuHi  no  paieHCTay 
saoioRAel 

(5.52) 

0.463 

0,407 

0,416 

0,412 

(M 

noaoSNii  NTtpaaNf  no  Ki«ApaTNRHo< 

(i.a) 

0.432 

0,410 

0,413 

0,412 

noaoSNat  mepiuHa  no  paatncTty 
^nkuhA 

(5.S5) 

0.331 

0.302 

0,394 

0.412 

KEY:  (a)  Method  of  solution;  (b)  Formula;  (c)  Approximation;  (d)  Accurate  oolntion; 
(e)  Coinnlex  iteration  by  Viarda;  (f)  Complex  iteration  with  a  variable  narameter; 

(r)  oimilar  iteration  by  equality  of  areas;  (h)  Similar  iteration  by  quadrat. ic 
deviation;  (i)  Similar  iteration  by  equality  of  functions. 

Note:  Numbers  of  Formulas  Indicated  are  for  Chajter  3. 


2.  Vibrations  of  Rods 

Problem  has  had  numerous  engineering  apnlicatlons,  especially  in  t.iorbonachirios 
in  calculating  the  vibration  of  vanes.  Following  presentation  reiera  mainly  to 
problems  on  vibration  of  vanes,  blades  of  propellers  at  cetera. 

We  shall  consider  at  first  the  natural  vibration.  The  boundary  int.egral 
equation  of  flexure  vibrations  of  a  cantilever  rod  (vanes)  relat.ivo  t.o  omnlitude 


(2.1) 


flMCUTM  has  tha  form  ^Sas  Chaptar  3t  equation  (2.22)^ 

'  ^  iV  i? 

•  "''1 J  ir^  S 5 W  *  (^)  O'.  Or, dr,. 

0*  f,  f , 

where  P  it  the  angular  frequency  of  the  vibrations. 

In  deriving  this  equation,  there  is  used  the  equality 

Af,  (r)  (r,)  i  (r,)  </r,  dr^. 


dr* 
It  R 


/  •. 


(2.2) 

(2.3) 


Equation  (2.1)  is  valid  at  .(0  -  const  or  in  a  general  eaae  in  Ignoring 
the  influence  of  natural  torsion  which  is  admissible  in  determining  the  first 
frequency. 

From  relationship  (2.2) 


n-i 


f  A  (rt)  ^  • 


(2.4) 


for  a  rigid  fixing  of  the  root  section 


From  equality  (2.3)  we  obtain  a 


moment 


-^(r,)=0. 

dr 

boundary  integral  equation  for  amplitude  bending 

(2.5) 

'  'a  'a  '  ' 


or  in  a  bridged  form 

(2.6) 

With  a  hinged  fastening  of  the  rcxi  from  equalities  (2.4)  and  (2.3)  we  obtain 


where 


R  R 

^  ]  J  P^'  (''i)  (''I  -  r„)  </r,  dr^. 


Methods  of  solving  homogeneous  boundau’y  integral  equatior.s  are  reviewed  in 
Chanter  3. 


1Z5 


l-te  shall  nrasent  results  of  the  calculation  nertaininp,  to  rod  of  o. instant 


section.  The  frequency  is  determined  by  the  fomula 


The  k  values  are  given  in  Table  7.  The  values  fe,“3,52;  fej  =  22,03;  —  ar'^ 

accurate. 


Table  7.  Calculation  of  First  and  Second  Frequencies 
by  Heans  of  Boundary  Intsgial  equations 


1  iAdepaad  aacroia 

1 

Biop 

]>4  SlCTOTa 

Cft)MeToa  ptcwi  fgVnpHi. 

.wcHrie 

'*’\oqnoe 

l(^)  npHfiJutiKeNHe  { 

Till  not 

1 

9 

peiueNit«| 

1 

1  ^  1 

3 

lief  It  M 

(d)  CpawieHNe  opaiwir 

S.SI 

3,83 

3,82 

24,88 

22,8? 

52,90 

22.03 

UD  PaeeNcno  cMiapwel 
Mopew 

{-/.)  MHMMyM  vujipaTHSHe* 

°  ro  emaoHeNHi 

3.M 

3,82 

3.52 

23,07 

22.76 

22,73 

22,03 

3.82 

3,82 

3,52 

1 

22.64! 

i 

22,91 

22,92 

22,03 

K3Y:  (a)  Method  of  calculation;  (b)  First  frequency;  (c)  '.econd  frequency;  (d) 
Comoarieon  of  ordinates;  (e)  Ekjuality  of  scalir  norro;  (  f )  Minim’.iin  of  Square  deviation- 
(p)  Anproxiisation;  (h)  Accurate  solution. 


The  initial  aoproxination  was  selected  in  the  form 

» I  —  C. 

iiriiere  "7  •  where,  1  is  length  of  rod. 

In  the  calculation  the  rod  was  subdivided  into  10  sectors  and  the  intefrat ion 
was  made  by  the  traoezoidal  rule.  Sufficient  accuracy  is  obtained  also  in  diviainp 
the  rod  into  5  sectors. 

Condition  of  orthogonality  for  equation  (2.1)  has  the  form 


for  equation  (2.5)  corresnondingly 


(2.«) 


(2.9) 


The  determination  of  the  second  frequency  and  form  of  vibr  :■>,  ia  pointed 
out  in  Sec.  3,  Chao.  3.  Thue,  for  example,  equation  relative  to  bending  moments 


2Z^' 


has  the  form 


a 


(2.10) 


Results  of  the  calculation  are  given  in  Table  7. 

From  the  table  it  is  evident  that  the  error  due  to  use  of  trapezoidal  rule 
(vdth  10  sectors)  is  larger  than  error  from  an  "incomplete"  convergence  of  the 
nrocess . 

Anplication  of  normal  integral  equations 


e +  5  (/?)/.  +  ^  (/?)/. 


has  been  considered  in  Sec.  k.  Chanter  3. 

For  a  rod  of  constant  section  we  obtain 


(2.11) 


5(C)=/»*-^^ 


till 


i  + 


+:(i)i-^(i)/(i-;);  ('=V^)- 


The  functions  ®oG)  and  4^1  (^)  are  equal 

.... 

where 


1  /‘p/' 

X~p*  -  ■  - 
^  £J 


)• 


(2.12) 


The  frequency  equation  for  the  first  approximation  /"in  the  series  (2.12)  are 
retained  bj-^  two  of  the  first  terms_7: 


I  -  0.08334Z  +  0,000347/-  -=  0. 


Tr.e  fir»t  fonroxi-mation  for  the  first  fraqutnc.’’ 


"1  'fera  fron  the  accurate  by  1.7'.'. 

In  th«  aecord  annroxiinat ion  for  I'urot  th.-ee  frequo-iona ,  vie  obL<i..n  vajucs  of 
th#  coefficient  k  /"equality  (2,7)_7: 

4i(j)=332;  "33,6. 

In  the  th*rd  annroximation 

4io)"3.52:  *j,„»66.I3. 

We  shell  coneider  now  the  vibrations  of  a  natinrally  twisted  ri-d. 

If  “(')  and  v{r)  are  the  amolitude  disnlaceme.nts  of  nojrts  on  the  vatu's 
axis  alo-.,'-  the  axes  x  and  y,  then  the  bendin{^  moments 

(f)  —P'  J  J  (rt)  V  (r,)  </r,  dr, 

r  e 

ar  d 

M,Ar)~rJlpri  r%)u(rt)dr,  dr^. 


In  accordance  -rfith  equality  (1,11)  wa  obtain 


Af,(f)  JJ(>F {r,)u{r,) dr,  dr^  f 

+  llaa  j  Jp/^(r,)  v{r^dr,  dr,  ^ . 
r  r,  / 


B’  inteprating  equations  (I.IO)  and  (1.11)  and  intpodici-p  u.o  i:,to 

ecc’ifty  (2,13),  we  shall  have 

(as  5 

cos  *  J  j  pr(r,)  J  j  cos  a  (>,)  dr,  dr,  dr.,  dr, 


r  'I 


•mne 


(k. 14) 


n\nr,d.'.,  i  , 


If  anrl«  *('’)  -  comt,  then  equatif)n3  (-^.14)  and  (<'.')  ''incide. 

At  an  angle  of  natural  torsion  of  rod  f)f  an  order  of  30-^0“,  the  torsion 
insignificantly  increases  the  first  frequency  (to  1:^)  and  considerably  lowers  the 
second  (15  to  ^0^);  this  is  confirmed  by  exoerimental  data. 

The  determination  of  the  second  frequency  is  made  by  equation  (2.10)  with  the 
operator  KM«i,  ,  corresponding  to  ecjuation  (2.i4). 

In  calculating  both  rigidities  for  flexiire  the  problem  reduces  to  a  homogeneous, 
boundary .matrix  integral  equation. 

We  shall  present  for  this  case  the  •  cid ; cjus  ^f  ortho. jonality 


'^"SL\dr^0 

I 

(l*J.  ij~ 


1  '>  1 
I  I  *-1 


.) 


The  calculation  of  torsional  vibrations  of  rods  also  reduces  to  solution  of 

homogeneous  boundary  or  normal  integral  equations. 

The  boundary  integral  equation  for  amplitude  angles  of  rotation  has  the  form 

r  1  r 

J  'cr(.,r  j  (^.) « (n) rf..  (2.15) 


where  GT{f)  is  strength  of  section  rf  rod  to  torsion;  /^(r) 
moment  of  inertia  of  section  relative  to  trie  center  of  rigidity. 
Corresponding  equation  for  amnlitiide  torques  is  such. 

'  '% 


is  the  polar 


(2.15) 


Fig.  11.  F  .n  ueriving 
eqia’ions  of  flexme- 
torsional  vilrat  ins 


r 


Conditions  of  orthogonAllty  written  out  In  the  form 


*Ul**^t  </r*=a 

or 


The  normal  aquation  for  torsional  vibrations 

•  W  - I f  («■.)•  W  *1  <<>•■+•  (»)• 


Ws  turn  now  to  the  calculation  of  flexure-torsional  vibrations  of  a  rod  (vane). 

The  origin  of  the  local  system  of  coordinates  Is  placed  at  the  center  of 
rigidity  of  section;  coordinates  of  center  of  gravity  in  this  cnorainate  ayntam 
will  be  (Fig.  11). 

In  designating  the  amplitude  disolacement  of  canter  of  rigidity  nt  and  ft, 
angle  of  rotation  9,  we  obtain  a  system  of  three  differential  equations 

^  j-(fA-fy,)ilaeeoi‘e^^^  4-  (fyicoi*e+£y,sin*e)  * 
^j(fyism»e+£y,eoi»f)  ^-(fyt-^y,)iin«cos 

(2.17) 


This  system  is  equivalent  to  the  homogeneous  matrix  boundary  integral  equation 


where 


(2.18) 


Equation  (2.18)  is  obtained  by  previoualy  indicated  methods,  and  the  value  of 
operators  is  not  written  out  hare. 

In  calculating  vanes  it  is  possible  usually  to  ignore  the  flexure  in  plane  of 
the  maximum  rigidity,  whereas  the  flexure-torsional  vibrations  are  described  by  tho 


following  equation  relative  to  the  bending  and  turning  moments: 


where 


Ik\  ^  [k„  a..  |[mJ’ 


(2. 


a  a 


KuM^  -  cos  «  j*  J  pA  (r,)  r  r 

r  ft  ’ 


(^) 


cos  a  (r«)  rfr4  dr,  dr,  dr,  + 


't 


+  **"  •  j*  j*  ('■«)  j  J 

aa  ^ 

A„Af.  =  »«n  a  j  J (r,)  pF(r,)  |  dr,  dr,  dr,  - 

f  ft  f, 

-CO.,  j i c„ (r.) pf  (r.)  J 

f  ft  tt 

A„yM,  =  j  e,,  (r,)9F  (r,)  |  j  sin  a  (r,)  dr,  dr, dr,  - 


'<  'j 


-  j’  t„  ('■.)  ff  w  f  I  “  ('■■>  *■• 


The  calculation  is  made  by  the  equation  (See  Chap.  3,  Sec.  3) 


f>w„oT 

=pi 

Am 

1  ^«(0  J 

r(/)  jL. 

*'n 

|[Af,(,_,)  J 

Under  the  condition  of  normalization  by  a  scalar  norm 


Second  frequency  of  fXexure-torelon*l  vibretiane  is  determined  frcro  the  eciuation 


where 


We  now  consider  the  fleiwri  vibrstions  of  a  naturaii^  tvfiated  rod  in  a  fial.i 
centrifugal  forces.  They  tre  dfsoribed  by  the  following  system  of  differential 
equations  with  respect  to  VMilitude  flexures : 

J2^.(e-ytcos»s  +  £y,sin*e) 

~  [(^A  ••••  S + £/,  foi»s)g  -  (f y,  -  ify,)iin  •  coi « -  ( -  • 

where  C(0  is  the  terielon  in  the  section  of  the  rod. 

In  the  particular  case,  ws  will  have  vibrations  of  an  untwisted  vane  in  axial 

plana 


(2.21) 


by  intefhrating  this  equality,  we  find 

In  considering  relationship  of  fom  (1.4),  ws  shall  obtain  a  boundary  integral 
equation  for  a  rigidly  secured  red 


(2.22) 


where 


ff 

f  r,  r,  f. 


rfr.dr,. 


liquation  (2.22)  is  a  two-parametric  integral  equation. 

It  is  nossible  to  show  that  if  equation  of  flexure  of  vane 

where  is  the  bending  moment  of  the  transverse  load ,  and  the  equation 


corresponds  to  the  problem  on  vibrations  of  non-rotating  rod,  then  integral  eqxjatio 
of  vibration  of  rod  in  a  field  of  centrifugal  forces  has  the  form  (2.22). 

In  composing  equation  (2.22)  there  may  be  made  assumptions  of  a  different 
nhysical  nature  relative  to  the  operators  K  and  For  example,  in 

determining  the  first  frequency,  operator  can  be  taken  without  considering 

the  natural  torsion  (2.5),  and  the  operator  may  be  taken  the  same,  as  for  a 

vane  of  constant  screw  pitch  (1.20). 

In  solving  equation  (2.22)  by  the  method  of  successive  approximations  we 
obtain 

By  using  the  norm  of  function  on  basis  of  maximum  (method  of  comparing  ordinates) 
we  find 


The  orocess  of  successive  annroximations  is  unconditionally  convergent  with  the 
oarameter  of  flexibility  »<I.  For  nore  flexible  vanes  one  should  apply 
corresDonding  iterative  processes.  Thus,  by  applying  methods,  similar  iterations, 
we  obtain  an  equation  of  ordinary  structure 

(2.24) 


1  hi 


whart  by  alollar  itaratlon  on  baals  of  aquaUty  of  functiona 


by  similar  iteration  on  baaia  of  equality  of  areas 


by  similar  iteration  on  basis  of  minimum  of  square  oeviation 


(j_u  =  (i-i)  —  •' 


n(i- 

't 


In  solving  equation  (2.24)  there  are  used  methods,  indicated  for  the  vioration  ov 
non-rotating  rods. 

Above  ware  considered  the  natural  vibrations  of  rods.  Let  us  turn  now  to 
nroblem  about  forced  vibrations,  at  first  without  taking  into  account  the  forces  of 
damping . 

As  example  \n  shall  take  the  flexure  vibrations  of  a  rod. 

Suppose  onto  the  rod  is  applied  an  external  excitation  load. 

g**q{r)cos  >/. 

The  equation  for  amplitude  flexures  of  forced  vibrations  has  the  form 

KH 


*('>’■' a  s !  ^ 

ft  ft  '  '•  't 


or  in  abridged  form 


I-'^/CE+Z. 


(2.25) 


HZ 


where  f  is  the  flexure  of  vane's  axis  under  action  of  distributed  load  q(r). 
Process  of  simple  iteration 

=  -1)  +  / 

converges,  if  ''  where  p^  is  the  first  natural  frequency. 

At  one  should  apnly  the  previously  indicated  iterative  processes. 

Raoidity  of  convergence  deoends  on  form  function  f  and  at  v<4pi  it  is 
obtained  usually  entirely  satisfactory. 

The  very  best  results  in  a  number  of  nractical  examoles  wore  given  by  the 
method  of  similar  iteration  on  oasis  of  minimum  of  square  deviation. 

Iffe  now  consider  question  of  determining  the  coefficients  of  dynamic  rigidity 
of  rod  of  variable  section. 

Suppose  onto  the  root  section  of  rod  there  are  applied  dynamic  (Fi^.  12). 

Q=Q«cos  'tt, 

A!=AfoCos  V  t. 

This  will  cause  a  vibration  of  entire  vane,  with  which  in  root  section  there  will 

be  a  flexure  I  and  angle  of  rotation  " 

{  =  l«cos 

There  exist  the  linear  relationships 

Q«  =  “ii’o  +  ®it'0*  (2.26) 

Afo  =  *ii5o  +  *ji ’O' 

The  coefficients  are  called  coefficients  of  dynamic  rigidity.  They  possess  the 
property  of  reciprocity 


Fig.  12.  Determination  of  Coefficients  of  Dynamic  Rigidity. 
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Wt  sheXl  wit*  out  tha  norm*!  integral  equation  for  considered  case  in  following 
form: 


where 


I  -  >*/V;  +  ij,  +  +  Af„/,  +  gj/„ 

^  *  J  J  fJ,  ('ti  J  J  ^ 

'i '%  ’  r,  r, 

/o“l:  fi  =  r-r,^ 


(2.2?) 


The  solution  of  equation  (2.27)  will  be  such; 


where 


«  =  5o'I‘o  +  5;'l>j  +  /M,<l*,  +  Qo«V  (2.23) 

«I»,-/,  +  vW/,  +  n*A^/,+  .  . 

(i  =  0,  I.  2.  3). 


^  intrv./duclng  the  values  (2.28)  into  the  boundary  conditions 

Q(/?)=0;  ilf(/?)=0,  ' 

we  arrive  at  equalities  (2.26),  where  the  coefficients  a„  become  knowr.. 

lat  us  consider  now  forced  vibratiwis  with  a  consideration  of  linear  damping. 
The  differential  equatio.i  of  the  problem  is  written  out  as : 

^  1^)  ^  ^  ^ 

By  nutting  the  solution  in  the  form 

{-■  w  (r)  cos  »/+z  (r)  sin  »r, 

We  arrive  at  a  system  of  equations 

^  (£y,.^P)_  ,tp/r (,)  ^  ^ 

5  “  **  ~  '’*P^  ('■)  *  ('')  =  0- 


i3^ 


For  a  rigidly  secured  vane  this  system  is  equivalent  to  a  inhomogeneous  boundary 


matrix  equation 


where 


f  r,  HR 

KiiW  =  J  J  (^)  w  ('  4)  dr,  dr^  dr,  dr^, 

r*  't 

'^••^“-’11  -sTm  1 1  *  <"•> '  <'•> 

r,  r,  r,  r, 

f  r,  R  R 

“  'J  J  ttst  1 1  *  <'•>  ” 

f» '»  f*  't 

r  r,  R  R 

J  J  .f  i *• 

r.  %  ’  r,  r, 

r  ^  a  /? 

/t“0. 


(2.30) 


Equation  (2.30)  at  '*<^x  is  solved  by  the  method  of  simple  iteration,  at  »>Pi 
there  can  be  applied  the  method  of  similar  iteration  (See  Sec.  k.  Chap  3)« 

3.  Critical  Speed  of  Shafts 

The  determination  of  critical  speed  is  important  for  many  high  speed  machines, 
expecially  turbomachines. 

Let  us  consider  the  general  case  of  the  precessional  motion  of  shaft  (Fig.  13). 
Sunpose  the  plane,  containing  elastic  line  of  shaft,  revolves  with  angular  velocity 
V,  and  the  shaft  itself  is  rotated  in  reference  to  this  plane  with  an  angular 
velocity 

The  kinematic  model  of  a  similar  motion  is  shown  in  Fig.  14. 

Angular  velocity  of  shaft  is  equal  to 

(3.1) 
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Th«  »ngulAr  v«loci*ies  »  and  K  are  presented  in  the  following  form: 


(3.2) 


where  •  oerhaps,  in  ge-  ral,  is  an  arbitrary  (real)  number  . 


Fig.  13.  Processional  motion 
of  shaft. 


Fig.  14.  Kinematic  model  of  nre- 
cossional  motion  of  shaft. 


The  disk  acts  on  the  shaft  with  a  stress  and  a  moment 

(o<). 


(3.0 


Af,-.*p,y  (ft,). 

/ _ !_/ 

where  m^  is  the  mass  of  disk;  */  3  equatorial  moment  of  inertia 

of  disk. 

The  equation  of  the  stability  of  revolving  shaft  with  distributed  masses  m(x) 
and  moments  of  i.iertia  I(x)  has  the  form 


^  S)  “  “*  ^  ^  “  (2 -  •)  £(/  ( 0  v'K-t))  j  0.  (3.3) 

Equation  of  flexure  vibrations  of  shaft  with  a  calculation  of  inertia  of  turn 

£:  ( W  S)  “  <  '>  ~  “0 

At  /(x)*=0  equations  (3.5)  and  (3.6)  coincide. 

Solution  of  equation  (3.5)  can  be  used  and  as  solution  of  equation  (1.6),  if 
it  is  assumed  e=— I  (reverse  synchronous  procession)  and  to  decrease  Ari  the 
magnitude  of  l(x)  by  three  times. 

We  now  turn  to  composing  integral  equations  of  the  problem.  Lot  us  consiaar 


♦The  nrecession  is  called  synchronous,  if  |«|-|,|  In  in  identical  direction 
wand  V  the  precession  is  considered  forward,  with  a  different  direction — reverse. 
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AS  an  exAfflDlo  a  shaft  on  two  pivoting  bearings  (Fig.  13),  carrying  the  distributed 
masses  m(x).  The  soluticai  can  be  applied  also  in  the  presence  of  individual 
masses,  if  the  siass  of  the  disk  is  distributed  along  the  length  of  corresponding 
section.  The  gyroscopic  effect  we  disregard,  and  the  parameter  g=l.  Bending 
mcnent  in  section  x  will  be  equal  to 

M  (x)  =  R^x  +  <u*  f  j  m  (x,)y  (.v,)  (/x,  </x,. 

By  determining  reaction  in  the  le‘*t  support  from  condition  M(l)  =  0,  we  find 

M(x)  =  m'A,, 

vrhere 

»»  y  (-"f*)  </x,  —  -j  f  J  m  ( x,)  _y  (Xa)  rfx,  rfx,. 

0  0 


By  using  the  equation  of  flexure  dx-  EJ (x)'  we  obtain. 


y(x)' 


I*  f  (*  ■^y 


<:!^rfx,rfjc,+/(0)x+y(0).  (3.7) 


By  determining  y'(0)  from  condition  y(0—0  (magnitude  y^Q^sO), 
at  a  hanogeneous  boundary  integral  equation 


we  arrive 


■iJ 


Ay  (Xt) 

£J(x,) 


</x,rfx, - - 


I  X, 


A,(Xy) 

£J(xj) 


dx,  dxg. 


(3.8) 


Ws  now  consider  the  general  case  (Fig.  15),  *dien  there  are  considered  the 
distributed  and  concentrated  masses.  The  transverse  force  in  section  of  x 

Q  (*)  +  2,^ (•«.  *<)  +  (^.  <*i) (3.9) 

where  the  unit  function,  for  example,  S(x,  b,),  is  determined  by  the  equality 

(1  x>b,. 

In  considering  the  dependencies 

q  (x)  =■ «.’«  (x)  y  (x) ;  a  (x)  =  e*m  (  v). 

P,  =  (bi), 

Ws  shall  write  out  expression  (3.8)  in  the  following  form: 

Q  (X)  - (X)  +  S  5  (X.  a,)R,  (3.10) 


:Z37 


(3.11) 


and  for  tha  banding  maoant 

Af  (jc)  -  .M,,  (x)  +  i  S  (X.  a,)  R,  (x  -  a,). 


Mhera 


Fig.  15.  Shaft  on  two  pivoting  bearings. 

(Jf)  *=]  «  (-*^1)  y  (JCi)  rf  'fi  +  i  S  (X,  bi)  a,> 

a  ('ft)  J'  (-ft)  rfx,  rfx,  + 1?  (x,)^'  (Xi)  rfx,  + 

+  J^5(x.  (b,)  (x-b,)  +  ^y  (^)l 

(p(x)».(2-.)/(x)l. 


(3.12) 


datarainina/ 

In  /  “faction  and  R2  from  conditions 

Q(/)»0.  A1(0-=0, 


we  obtain 


whara 


(3.13) 


Af  (x)  (x). 

A,  (X)  -  (X)  +  5  (X.  a,)  [  A.,  (/)  (/- aj  -  (/)]  + 

+  S(x,  a,)-^^[A,,{l)-A„(f)(t-a,)]. 

A,  —  fli 


(3.1^*) 


In  datanaining  in  aquality  (3.7)  the  magnitudes  y  (0)  and  (O)  from  conditions 
v(ai)»0  and  >(a,)«=0,  wa  obtain  a  system  of  integral  equations 

(3.15, 


where 


X _  ?(* 


^dx.dx,- 


I— I 


l^dx.dx,- 


dx^dx^  ; 


•l  •»! 

r  r 
1) 


K^i)y  =  [  ^M!l>-  rfx,  +  — ! —  r  f  dx,  dx^ 

•  )  fy(jr,)  a,- a,  J  J  £J(x^ 


(3.16) 


•«  'I 

-IH5 


EJ(x,) 


rfXjrfjC,  . 


(3.17) 


Calculation  by  formulas  (3.16)  and  (3.17)  are  very  simple,  since  they  contain  all 


two  integral  operations 


(fr 


II  " 


£J(xm} 


Equation  (3.15)  expresses  the  matrix  integral  equation 

1><“J  L/f..  AfnJlH 

where 

In  abridged  form,  equation  (3.15)  is  written  out  as: 

[y]. 

In  solving  by  the  method  of  successive  aporoximations 


(y^ol  = ‘“(.j  (y(i-ol. 


First  line  of  this  equality 


By  method  of  comparison  of  ordinates 


,  >(1-1) 


(3.18) 


(3.19) 


(3.20) 
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wh6r0  is  ths  abscissa  of  section,  corresponding?  to  the  Tnayimiim  value  Ly(i-i)l- 

After  determining  w  find  y,„  and  y,'„  and  further 

Practice  of  calculation  showed  that  there  is  sufficient  not  more  than  two 
approximations  (second  approximation  is  for  control). 

Thus,  there  is  determined  the  mirimum  in  absolute  value  eigenvalue. 

>fe  shall  d»«ll  on  one  circumstance  associated  with  the  calculation  of  systems 
with  strong  influence  of  the  gyroscopic  effect  (for  example,  with  disks,  located 
near  supoorts). 

It  may  be  found  (in  oractical  cases  extremely  rarely)  that  “•|<0. 

This  means  that  real  angular  velocity  will  be  greater  thart  |‘•l|  and  must  be 
determined  by  taking  into  account  corresponding  condition  of  orthogonality  (see 
determination  of  second  critical  sneed). 

If  true  value  is  •>*>(), 

but  in  the  first  approximation  in  view  of  unsuccessful  selection  of  y(0)  there  is 
obtained  n 

then  one  should  continue  the  process  further  and  it  will  converge  to  a  real  angular 
velocity  “i. 

In  engineering  problems,  the  indicated  cases  may  be  encountered  as  exceptional, 
and  only  with  a  calculation  of  gyroscopic  effect. 

For  a  number  of  problems  encountered  in  practice  of  (rotor  with  large  number 
of  disks,  calculation  of  mass  proper  of  rotor  et  cetera)  the  calculations  of  the 
discussed  method  are  found  to  be  significantly  less  laborous,  than  calculation  by 
other  methods  (for  example,  requiring  the  determination  of  influence  coefficients). 

Ws  now  turn  to  determining  the  second  critical  angular  velocity.  Condition 

♦ 

of  orthogonality  in  considered  problem  has  the  form 

I  LviX/*  ('*)  ’p  (x)]dx=o 

(l  +  J,  2.  3.  ..). 

♦In  the  presence  of  concentrated  masses  and  monents  of  inertia,  integrals  are 
taken  in  the  sense  of  Stieltjes. 


whera  a  (‘if)  and  ?(x)  depend  on  distribution  of  the  masses  and  moments  of  the  inertia 


along  the  length  of  shaft 

aix)  =  e»m(x).  Mx)  =  t(2-t)/(x). 
The  calculation  is  made  by  the  equation 


where 


Pf 


Ul* 

l/C','>llyl  =  lAf'»l-Pi  (I'll. 

I  (>i«  ->i 


(3.21) 


Above  theire  has  been  considered  a  shaft  (rotor)  on  two  pivoting  bearings. 

Analogous  equations  may  be  comniled  also  for  other  cases.  Suppose,  for 
example,  the  supports  of  the  shaft  are  elastic,  then  there  exist  the  dependencies 


* 


y(ai) 

y{at) 


^ . 

K«i' 

3*. 
^  » 


(3.22) 


where  /C*i.  Kgt  is  the  rigidity  coefficient  of  the  support. 

If  the  sunport  is  a  complex  system  of  masses  and  elasticity,  then  the  magnitudes 

and  K,.  are  dynamic  rigidity. 

Relationships  (3.10)  and  (3.12)  remain  in  force,  and  from  conditions  (3-13) 
the  magnitudes  and  R2  are  determined  by  Au  and  At,. 

The  magnitude  y(0)  and  y'(0)  in  equality  (3.7)  are  found  from  condition 
(3.22)  which  results  in  a  corresponding  integral  equation.  There  are  no  great 
difficulties  in  composing  the  equations  and  in  other  calculating  cases  (shaft  on 
several  pivoting  bearings). 

Let  us  consider  as  an  examole,  a  shaft  with  large  number  of  identical  disks 
on  two  pivoting  supports  (Fig.  16).  Per  unit  of  length  of  shaft,  there  should  be 
a  mass  a  moment  of  inertia  /. 


♦In  the  presence  of  a  connection  between  the  supports,  y(«i)  and  y(a,)  are  ex- 
nressed  through  the  linear  combinations  and  R2. 


IHl 


Equations  (3 •15)  havs  the  form 


>(x)- 


(■*») 


£J(xO 


iX^dXy  — 


(■*») 

£y(jt,) 


dx^dxi 


{■*%) 

£J(xO 


dXt  dXi 


By  virtxie  of  equality  (3*  14)  and  (3.12) 

W  “  J  ]  *  ('^i)>  (-<1)  dx^  +  I  p  {x^)y’  (x,)  dx„ 


where  for  forward  synchronous  precession 

As  an  Initial  approximation  >«“»(•*)  foj.  a  shaft  on  two  pivoting  bearings. 


It  Is  Dosslble  to  take 

jf(0)  =  C(x  (x  fl|),  (3.23) 


where  the  constant  C  expediently  Is  determined  from  the  condition 


In  accordance  with  equality  (3.23) 

ri 

tt 

ttttt:' 

>io  (Jc)-Cl2x-(fl,4-fli)l 

T  ^ 

n 

UTIL.. 

tI 

Shaft  was  subdivided  into  ten  sectors  and 

- - 1  ■  - - 

the  Integration  was  made  by  the  trapezoidal  Fig.'  16.  Shaft  with  uniformly 

^  .  distributed  disks. 

rule .  The  ratio  is  • 

In  the  first  approxismtlon  there  was  obtained 


(if  there  were  made  an  accurate  integration  »(i)=  10,56 
In  second  approximation 


* 


^  10.93 


mt^ 


i^z 


The  problem  has  an  accurate  solution 


From  the  calculation  it  is  clear  that  first  approximation  gives  a  deviation  of 
an  order  of  I.65S,  and  the  inaccuracy  of  the  second  approximation  (1.8^)  is  ex¬ 
plained  by  error  in  the  anproximate  integration.  It  can  be  removed  by  the  selection 
of  a  greater  number  of  sectors,  however,  a  great  accuracy  in  the  calculation  is  not 
required . 


4.  Stability  of  Rods 

The  application  of  integral  equations  in  problems  of  stability  are  especially 

effective,  since  for  practical  purnose  there  is  required  a  seeking  of  only  the 

«■ 

minimum  eigenvalue  . 

Wb  shall  consider  the  stability  of  a  rod  of  veu'iable  section  under  the  effect 
of  concentrated  and  distributed  along  the  length. 

The  differential  equation  of  the  problem  has  the  form 

'  (4.1) 

whore  P(2)  —  compressive  force  in  section  z; 

EJ  —  minimum  strength  of  section  to  flexure. 

Equation  (4.1)  is  valid  for  any  fastened  ends  of  the  rcxl. 

Ws  shall  take  as  an  example  a 
cantilever  rod  (Fig.  17).  Since  a  trans¬ 
verse  force  in  section  ?  -  1  is  absent, 
then  order  of  equation  (4.1)  can  be 
lowered : 

~  *The  application  of  boundeuy  integral  equations  to  problems  of  stability  were 
for  the  first  time  given  by  Yu.  V.  Repman. 


Fig.  17.  Stability  of 
Cantilever  Rod. 


In  the  considered  case 


U  (-1)  +  S  -5  (z,  a,) 

t 

where  Pi  is  the  concentrated  force  in  section  z=Oi. 

The  single  function  S(z,  a<)  is  determined  now  by  equality 

'  ''10  t>a,. 

If  the  stress  Pi  is  directed  for  elongation,  then  into  the  calculation 
equation  it  should  be  introduced  with  a  minus  sign.  In  problems  of  stability, 
external  loads  contain  as  factors,  parameters  subject  to  aetenrination. 

Thus, 


f  ^  (*i)  \f\  **  (^>  ^i)  Pi' 


(4.3) 


For  the  concentrated  forces  as  parameters  X,  visually  there  are  ta^vsn  the 

magnitudes /*j.  i.  e 4  (i=I,  .  .  . /r).  Frequently  as  Xj  it  is  convenient  to  take 

„  -  \ 

the  dimensionless  parameters,  for  example  • 

For  concrete  calculation  ore  should  note  denenaence  betv.ee  ^iriAi;.et'.irs  so 

that  expression  (4.3)  contains  only  one  independent  parameter  ,V.  (for  example, 

X,<^X,  X,  aO.SX,  1,2/  et  cetera). 

.As  chief  unknown  we  shall  take  the  valte 

Equation  (4.2)  now  will  be  written  out  as: 


\  *-»  i 


B>''  integrating  both  sides  of  equality  from  z  to  1  and  considering  the  boundary 
condition 


4(0-0, 


obtain 


*  *  ^ 

^  ^  j* 9  (’i )  /u  (*i )  4  ^  \S (z,  a)  J  ^  (z,)  dz^. 


i— I  «  (4.4) 

We  note  that  on  the  right  hand-side  of  the  equality  is  the  expression  of  bending 

int.  It  can  be  determined  also  directly  from  consideration' of  Fig.  17: 

M  (z)  -  j  ^  (zj  Ly  (Z,)  - y  (z)  1  </z  +  PfS  (z.  <i()  (y  (n,)  -y  (z)  1 ;  ^  ^ 


since 


Then  equations  (4.4)  and  (4.5)  in  accuracy  agree.  After  dividing  both  sides  of 
equality  by  EJ  and  integrating  then  frcMi  0  to  ^  (?(0)*=t)),  obtain  the 

integral  equation  =  ,, 


Ooerator 


(’*)  r  ('*)  dx.  dz^  +  V  ifi  (z,  n,). 

i-1 


where 


B  (ar,  a,)  = 


(»  =  0.  1 . n). 

‘  'l 

i  a  - 

g  "I 

f— 


\)dztdzx 


z  <  a. 


One  of  the  coefficients  vi  may  be  taken  as  arbitrary  (for  example,  ^•*“0- 


Fig.  18.  Stability  of  Rod  of  Constant  Section. 

Expression  for  B{z.  Oi)  it  may  be  written  out  also  in  a  simpler  form 

t<a,  a, 

0  I. 

if  it  is  agreed  at  z'^a^  to  maintain  the  same  value  of  the  function  B(*  di). 
The  integral  equation  (4.6)  has  a  simple  structure. 


1^/f 


It  is  interesting  to  note  that  the  use  of  differential  equations  would  re- 
qu’Te  the  construction  of  solutions  in  each  sector  and  a  linking  of  the  solutions 
with  a  calculation  of  boundary  conditions. 

Let  us  consider  several  particular  cases. 

Suppose  we  have  a  rod  of  constant  section  with  force  at  end  (Fig.  13, a). 

For  this  case,  equation  (4.6)  will  be  such: 


«  ' 


or 


where 


t(^)  -  =  — j  (2,) 

0 


EJ  i 


The  accurate  solution  is  1  =  2,467. 
As  the  initial  approximation  we  shall  select 


satisfying  main  boundary  conditions  (0)  =0,  (1)  =0. 


(4.7) 


B:'  the  method  of  ccnparing  ordinates 


►(!)  = 


y<o) 


^T(0) 


2,500(1,34). 


By  the  method  of  minimum  square  deviation 

I 

f  t(0)^T<0)^ 

=  2,467  (0,Q0). 


In  the  following  aooroximation  by  the  method  of  comparing  ordinates 

X(*)  =  2,471  (0,16). 

If  even  we  select  a  rougher  initial  approximation 

?(0)-j.  (4.8) 

then  ty  the  method  of  comparing  ordinates 

^1)  =3.000  (21,6);  >,j,  =2,500  (1,34);  X,3,  =  2.471  (0,16). 


In  parentheses  is  shown  the  error  in  %. 


If  one  were  to  apply  approximate  integrating  by  trapezodial  rules,  as  is  done 
in  practica'*  calculations,  then  by  subdividir  into  ten  equal  sectors  we  obtain 

^(1)^3,008;  i(!)=  2,511;  i,3)=  2,482. 

B:-'  determining  with  the  initial  approximation  (4.8)  by  the  method  of 
minimum  square  deviation,  we  find 

^1)  =2,470. 

For  a  rod,  loaded  by  a  distributed  load  (Fig.  18,  b),  equation  (4.6)  will  be 


such: 


M  I 

»  (*)  =  -^  J  f  (^-  ^»)  ?  (^t)  dzt  dz^. 


0  t. 


since  the  compressive  force  in  section  z 

^(*)  =  >*/o(»)  =  <7(/-z) 

In  converting  to  dimensionless  form, 

(1-E,)  ?(;,)</;*  rf'.. 


where 


EJ 


After  taking  a  rough  initial  approximation  in  the  form 

?(0)  =  5, 

v/e  obtain  by  the  method  of  minimal  square  deviation 

i(i)  =  9.05, 

=  8-01, 

A<3)  =  7.93 


(4.9) 


with  an  acciorate  value 

X  -  7,837. 

For  a  rod  under  the  action  of  two  concentrated  forces  (Fig. 18b)  frcm  equation 
(4.6)  we  will  have 


?(z)  = 


EJ 


(4.10) 


*1  0 

By  setting  uo  an  initial  approximation  f'O*  ~  we  obtain  by  the  method  of 

P 

comparing  ordinates  X<i)  =  — ^—=2,67;  i(2)=2,12S. 

The  accurate  value  X=s  2,068, 

Vfe  consider  now  the  integral  equations  of  stability  of  thin-walled  rods  of 


J.H7 


constant  section. 


In  the  presence  of  distributed  longitudinal  and  transverse  stresses  the 
problem  is  described  by  the  follovdng  system  of  differential  equations,  obtained 
by  V.  Z.  Vlasov: 

-  ( N  (■;'  f-  J'  +  (,-W/)  •  0.  (4.11) 

-  (/V  (r/  -  u  V)\'  -I-  (11^6) '  =  0.  ^ ^  ^  J 

fy.e'v  _  on- -  l(r/V4-  6' I'  -I-  \q^  + 

+  a.  +  ^.Vf,V  =  0. 


where  I  and  tj-  —  are  comnonents  of  displacement  of  center  of  flexure  along  the 

nrincinal  axes  x  and  y  arising  with  the  loss  of  stability; 

•  —  is  coranlementary  angle  of  rotation  of  section  dui'ing  loss  of 

stability; 

S,  Af, —  are  the  normal  (tensile)  force  and  the  bending  moments  in  the 
section  of  rod  under  action  of  external  load; 

—  is  sectorial  moment  of  inertia; 

T  —  is  the  geometric  rigidity  to  torsion; 

(Jj,  Vj,  —  are  comnonents  of  transverse  distributed  load; 

—  are  coordinates  of  point  of  anplication  of  distributed  load  in 
plane  of  the  section; 

in’  **»'  —  coordinates  of  center  of  bend; 

'’i  and  Pjf-  —  are  geometric  characteristics  of  section. 

In  comnosing  the  boundary  integral  equation  as  the  chief  unknowns  it  Ls 
exnedient  to  take 


r< 


life  now  consider  for  example,  cantilever  rod  (Fig.  19)  with  free  upper  and 
rigidly  fixed  lower  sections.  In  this  case,  we  shall  have 


5' W=JX«,)</c,: 
0 

0 


'  ""o  /  ^ 


Fig.  19.  General  case  of 
stability  of  thin-»#alled 
rod. 


V,e  shall  vrite  out  equatr  ons  (4.11)— (4.13)  in  the  folloiiring  fora; 


\r 


\r 


A  «" 

/  i 


I 


(u.u) 


F  £  ^ 

■1- 


+  r*(A^d)' 


+  M/i  ~  ((2?^Af,  -  2?^^)  by  + 
+  (9s  -  Os)  +  9,  (e,  -  a^))  J  d  (j,)  j  ^ 


(4.16) 


In  integrating  both  sides  of  equalities  (4.14)  and  (4.15)  twice  from  z  to  1 
and  by  considering  boundary  conditions  at  z  1  we  obtain 


I*  (‘:)  *,+  /-‘d  (Z^)jdz^~ 

•“ ^ I “ J j (2»)  +  U*)  ) dZi dz,  + 1 (2p^Af^ _ 


(4.18) 


-2p^,)»(*,)^,  +f(?,(^,-«,)+?,(^,-<i,)yf 


/s' 


^(zt)dz,dz, 


-fj-a(..)... 

< 

8y  integrating  eq-oality  (4.19)  from  0  to  z,  we  find 

•  <  !  ‘\  *» 

•  -  -  J  J  ^  (^i)  j  j  + 


(4.19) 


+  (*i)  j  <fa,  rf»,  -  ^  I  -  (  f  f(  (r,)  +  (z,)  )  rfz,  dz^dz,  + 


J 


+ j  f  (2P/»,-  2p^,)  »  (z,)  rfz,  </z,  + 1 J  (z,  _  a  J  + 

+fg  (fr-«r)  )f  •  ('•)  ‘'^1  -  iij  j  » (^.)  dz,.  J 


A 

practical  problems  ,  external  loads  usually  can  be  presented  in  the 
-lowing  form: 


^(*)“  -Mi(^)+),«,(z). 

(z)  «=  (z)  +  i«m,,  (z), 

M,  (z)  =  X,m,^  (z)  +  (i), 

« J  +  g,  (e, -  a^)  ^ X,/  (z). 


(4.21) 

(4.22) 

(4.23) 

(4.24) 


-•  parameters  i,  will  be  subject  to  determination. 


6lu.tlon.  (4.17),  (4.18)  «<1  (4.20)  «•,  wrltt.n  in  Mtrix  form 


^JU  f(jta  j 
^jn 


9 

9 

d 


0  0  0 

'  9 

+ 

0  0  0 

♦ 

.0  0  K,„. 

0 

.  (4.25) 


wher«  values  of  the  ooerators  are  readily  established  after  a  cooiparison  of 
corresponding  equations . 


If  one  were  to  introduce  the  matrix-column  then  equation  (4.2!>)  can 


be  written  out  even  more  briefly: 


(4.26) 


Equation  (4.26)  is  a  five-parametric  boundary  matrix  integral  equation. 

For  a  conci^te  problem  there  should  be  known  the  relationship  between 
narameters  of  load  a  then  we  obtain  the  two  parametric  equation 


where 


IA1=^2  VlAOI. 


(4.27) 


Two-Darametric  integral  equation  already  has  been  encountered  in  nroblem  on 
vibration  of  rod  in  a  field  of  centrifugal  forces. 

If  the  equation  (^]  (Aol  |4>]  has  the  eigenvalue  then  for 

the  solution  of  equation  (4.27)  there  can  be  used  the  method  of  simple  iteration: 


The  value  No 


maximum  values. 


(•(ol“"NolA')  [0(1-1)]  +  (Afoi 

is  sought,  for  example,  by  means  of  comparison  of  the 


=  +  V  + 


for  the  (i-l)-th  and  1-th  approximation  (3ee  Chap.  3,  Sec.  3): 


Then  we  obtain 


If  lA*il>N 

For  the  equation 


then  it  is  possible  to  apply  method  of  similar  iteration. 


we  have 


y(i-i-i).»“*(o - a 


2  / 


- - 

y<i)..  —^2j  ^0.^(0./ 

*  — “  1,  2,  ...fS. 


Imo  shall  give  an  example.  Suppose 
it  is  required  to  determine  the  value  of 
force  P,  causing  the  loss  of  stability  of 
rod  (Fig.  20).  The  force  is  applied  at 
the  center  of  gravity  of  section. 

In  equality  (4.21)  we  shall  take 

PP 

«,(«)=-!: 

Then  from  equation  (4.25)  there  ensues 

Fig.  20.  Stability  of  Thin-walled  (the  subscript  ^  we  omit) 

Rod  Under  Action  of  Longitudinal 

Force.  f-X  (/C,.?  f /C.tl'  +  A'u*). 

•  — ■  X  +  K„’}f  +  +  /Cojj®. 


(4.29) 

(4.30) 


(4.31) 


where 


Km^ 

Knf 


c 


■^W"- 


(WdCgdC,; 


“Ilf 

e  I 


cm 

EJ. 


The  calculation  of  the  presented  operators  is  comparatively  simple,  since 
they  contain  only  three  different  integral  expressions.  If  section  of  rod  possesses 
an  axis  of  symmetry,  (for  example,  the  y  axis),  then  center  of  rigidity  is 
located  on  this  axis  (‘*^  =  0).  Then  equation  (4.30)  becomes  indenendent 

t*-X  |]'<l'(Ct)<«.,d!:,.  (4.32) 

and  the  two  other  will  form  a  system. 

The  minimum  eigenvalue  of  equation  (4.32)  corresponds  to  the  Suler  force. 


(A. 33) 


In  the  solution  of  equations  (/*.27)  and  (4.31)  by  the  method  of  successive 

anproximatlons  as  the  initial  aoproximations  it  is  possible  to  take 

9,o)(t)  =  l-:.  »,o)(l)  =  C, 

what  satisfies  conditions  .  ^(1)=0,  d(0)=0. 


Further,  ws  shall  have 


In  a  similar  manner  there  can  be  considered  also  more  complex  questions  on 
the  stability  of  rods. 

5.  Extension  and  Flexure  of  Round  Plates  (Disks) 

The  indicated  problem  by  virtue  of  its  practical  importance  for  calculation 

of  disks  in  turbooiachines  has  been  investigated  by  different  methods.  However, 

* 

also  in  it, the  application  of  Integral  equations  makes  it  possible  to  construct 
one  of  the  most  effective  solutions. 

We  shall  consider  an  axially  symmetric  extension  of  disk  under  action  of 
centrifugal  forces  and  nonuniform  heating  (Fig.  21).  Parameters  of  elasticity  of 
material  of  disk  (E  and  p  )  are  assumed  depending  on  radius. 

Problem  is  described  by  a  differential  second  ordor  equation  relative  to  the 


*R.  3.  Kinasoshvill,  Calculation  for  Strength  of  Turbcoachine  Disks.  Defense 
Ministry  Publ.  House,  Moscow,  1954;  I.  A.  Birger,  Integral  Methods  of  Calculation 
of  Disk,  Collection  MAP  No.  6,  Defense  Ministry  Publ.  House,  Moscow,  1950. 


IS- 3 


radial  displacement  u{r)‘. 


where 


I  —  temperature  deformation;  q(r)  is  the  intensity 

of  body  force  (for  case  of  action  of  centrifugal  forces 

^(r)  =  fx..*r,  (5  2) 

here  P  is  the  density  of  material  of  disk;  w  angular  velocity  of  rotation),  or 
system  of  two  first  order  differential  equations  (equation  of  equilibrium  and 
equation  of  congruence) 


■f  f.ui*rA  =*  0, 

(5.3) 

Pa,)+af], 

(5.4) 

where  and  are  the  radial  and  circumferential  stresses. 


Fig.  21.  Extension  of  Disk, 
equations  (5.3)  snd  (5.4)  from  a  and  to  r,  ws  obtain 


It  is  possible  to  construct  different 
integral  equations  of  the  oroblem,  where 
for  practical  application  it  is  importaiit, 
so  that  equation  does  not  contain  deri¬ 
vatives  of  the  initial  oarameters  of 
disk  (h.  E,  dO. 

By  integrating  both  sides  of 


'  'a'  ~  j*  I  • 

r 

e» — —  —  (1  —  p)  0, — £  J -^^  (»*  —  a,)  rf'i — 


— £  (e  f  —  + -7- (a*  a  —  ►‘a/J  • 


(5.5) 


(5.6) 


The  subscript  *  in  these  dependence  indicates  that  the  value  of  the  parameter 
refers  to  r  »  a.  introducing  or  from  relationship  (5.5)  into  equality  (5.6), 
we  obtain  normal  Integral  equation  relative  to 

>(/•)  =  OI  (r)  — Or(0. 

This  equation  has  the  form 

•  • 

JL^  pu,*  jr,  A  dri-£(ar-a.U  +  -^(’*»-»‘«®~)-  (5.7) 


We  shall  wite  it  out  in  following  form: 

y  (r)  -  Qi  (r)  (r.)  y  (r.)  rfr,  +  Q,  (r)  {  (r,)  >-  (r.)  Jr,  +  (5.3) 

+  /+'*«o/ll 

where 

r 

/=  J  r,h  dr,  -E{at-  2JJ. 

S 

For  a  solid  disk  (’»«“= ’•a) 

For  a  disk  with  an  aperture  (’/.“O) 


The  value  >•«  is  determined  from  the  boundary  condition 

For  the  solution  of  equation  (5.8)  it  is  exnedient  to  apply  method  of  linear 
aoproximation  (See  Chan.  2,  Sec.  4),  vdiich  turns  out  to  be  in  the  given  case  more 
effective,  than  the  method  of  successive  approximations. 

The  problem  can  be  reduced  to  the  normal  integral  equation 


•I  (r)  =  ATsi  +  f,  (r)  +  (r)  +  9,  .F,  (r). 

j  A  (r.)  :>•  (r.)  dr,  ~ 

m 

p 


(5.9) 


where 

/S/\  r  I  i*(^) 

Calculation  of  disk  in  elasto-plastic  strains  can  be  made  on  the  basis  of 
equations,  valid  for  an  elastic  disk  with  variable  parameters  of  elasticity  .  This 
remark  refers  also  to  the  calculation  of  disk  for  creep  on  the  oasis  of  theory  oi 
aging.  Problem  about  synnnetric  flexure  of  disk  (Fig.  22)  has  much  in  common  with 
problem  on  extension. 

The  differential  equation  of  flexure  of  disk  with  a  calculating  of  nonuniform 
heating  through  thickness  of  disk  and  of  forces  in  middle  plane,  has  the  form 

ft  +  i  (in  rO)  +  [f -1:  (in  rO)  +  ^  (^) - 

(5.10) 

/(,)_ft(H-rt^(lnrO)+^['.^(l+n)]- 

f 

-A  Af  (r)  (1  +  f)  +  ^  J  ^  (r.)  r,  dr,  -  ^  (Q.  +  ^na?«)- 


In  this  equation 

y(r)  angle  of  rotation  of  the  normal  to  middle  plane  of  disk; 

g** _ cylindrical  rigidity  of  disk  on  the  radius  r; 

n{\-¥f) 

—  is  the  tensile  radial  force  in  section  r; 

*1.  A.  Birger,  Certain  General  ttethods  of  Solving  of  Problems  in  Theory  of 
Plasticity,  "Applied  Mathematics  and  Mechanics,  Voi.  l‘>,  'io.  6,  1951. 


910  —  is  the  distributed  load,  perpendicular  to  middle  plane  of  disk; 


^t(r)  —  is  the  difference  of  temperature  through  thickness  of  disk. 

Temperature  at  the  point,  at  a 
distance  z  from  the  middle  plane,  is 
assumed  equal 


By  using  the  equation  of  equilibrium  and 
equation  of  congruence  relative  to  the  bending  moments,  we  arrive  at  the  normal 
integral  equation 

Lt  /.\  _  KtAjI  I  B /,\ 

(5.11) 


where 


Af,(r)  =  A^Af,+  F(r). 

f, 


1  /)(r,)(l-i^*) 

r_3- 

'X  I. 

o(.j 

r 

S 

(0  =  I [(/ (0 -/ (r,))  (l-h^)rr- dr, - 

S 

f 

— J  Q  (^i)  dr,  -J-  Ja**/  (r)  -  aN,,  In  y  ]  + 

m 

m 

f 

Q  (0  “  — ^  J  9  (r,)  r,  dr,  -r  y  aQ,. 


now  turn  to  a  consideration  of  the  general  case  of  flexure  of  disk  (roxond 

plate) . 


equati  on  of 

nroblem 

has 

the  form 

o  ,  2+1* 

d*w 

1  dw 

3  d*m 

2 

d*w  1 

^  dr  \ 

r 

dr* 

r*  dr 

r*  dV  ^ 

r* 

drdV]  '  (5.12) 

4*D\ 

dw 

4.  •*  . 

-  '  ^  /aJkr 

dw' 

\ 

dr*  \ 

[dr*  ^  r 

dr 

+  r* 

J 

r  dry"^ 

dr , 

1 

*•* 

d>v 

-V*T, 

— 

r» 

le*  “ 

-9- 

iSl 


where  w(r,  I)  is  flexure  of  middle  nlane  of  disk. 
The  temoerature  of  the  noint  of  disk 


Hr.  •. 


*(r) 


the  function 


T(r,  e)  =  (l+^)a-YD 


The  radial  and  circumferential  ®»  stresses  in  middle  olaune  of  disk 
nossess  axial  symmetry. 

« 

For  comoosing  the  integral  equation  of  nroblem,  we  shall  use  the  eqxiation  of 
equilibrium  in  integral  form 

'  ' 't 

•  * 

+  -y  j  j (d-  r)  -h 

r.  /, 

where  Af^  M*  and  Afr$  are  the  radial,  circumferential  moments  and  torques 

oer  unit  of  length;  is  the  transverse  force  on  the  contour  r^b-, 

Wb  shall  consider,  as  example  the  case,  when 

^  ^  (0  cos  6.  (5.14) 

Forces  are  absent  in  middle  plane  of  disk,  and  the  heating  is  nonuniform. 

In  accordance  with  equality  (5.14) 

w(r,  e)=w(r)cose. 

By  introducing  values  and  M,t,  expressed  by  the  derivatives  ®'* 

into  equation  (5.13),  we  shall  obtain  a  boundary  integral  equation  relative  to 


f 


r  »  » 

Ft  (»'i)  +  ^  f  ^it  (''•)  dr,  dr,  -  J  J  V  (r.)  r,  dr,  dr,  - 

*  •  '  H 


(5.15) 


If  disk  (round  plate)  has  its  external  contour  free,  and  the  internal  fixed, 


then  the  magnitude  and  are  given, 

w(a)«a,  ^(a)=0. 


Equation  (5.15) 


^  =  /C'l*  +  f 


is  solved  by  simple  iteration,  at  XJ < I  and  by  similar  iteration  at  X  >-1. 

In  a  number  of  cases,  already  the  initial  approximation 

♦(0)  =  X 

gives  a  result  with  an  accuracy  of  an  order  10  to  15%. 

Equations,  analogous  to  (5.15),  can  be  composed  also  for  more  general  cases 
of  loading. 


6.  Symmetric  Deformation  of  Shells  of  Rotation 

The  problem  has  numerous  applic  ations  in  structural  engineering  and  machine- 
building. 

The  solution  of  problem  by  means  of  finding  accurate  solutions  of  corresponding 
differential  equations  encounters  great  mathematical  difficulties,  especially  for 
shells  of  variable  thickness.  In  connection  with  this  it  is  of  interest  to 
establish  Integral  equations  of  an  axially  symmetric  deformation  of  shells  of 
rotation  and  to  apply  approximate  methods  of  their  solution. 

The  scheme  of  the  shell  is  shown  in  Fig.  23. 

The  temperature  of  the  material  of  shell  is  assumed  linearly  variable  by 
thickness 

•  A 

where  is  the  temperature  of  points  of  middle  surface; 

A/  is  the  temperature  drop  through  thickness  of  wall. 

Ordinary  variables  are  used:  angle  of  rotation  of  normal  to  middle  surface 

_  a  dw  (6.1) 


1S9 


4nd  the  magnitude  il^i^aQ. 

tfhere  Q  —  transverse  force  in  section. 


(6.2) 


f. 


Relative  to  these  variables  the 
follovdng  system  of  differential  equations 
is  obtained 


(6.3) 


Fig.  23.  For  Calculating  Shell 
of  Rotation. 


Ek  \Rt^)tlj  *  (  Ck  */ 

(6.4) 

I  • 

In  these  equations  E  and  g  are 
narameters  of  elasticity,  variable  along 
is  the  cylindrical  rigidity.  The 


the  arc  of  meridan  s;  g*) 

functions  ^f(*)t  ^/(*)  and  ♦#(*)  are  determined  by  the  equalities 


u  1 

1  Tueosf Eh  tfs  \  stay  ; 

2ic£/lsla*y  J 

Safe  sisf  \  Islay/  2K£A/?|SlD’y  \  ^  Ri  / 


*1.  (fft^  \tft 

♦,(*)  — 0(1  + y-  ^*^[^**‘"»^(* 


(6.5) 

(6.6) 

(6.7) 


where  P  is  the  resultant  of  all  forces  (concentrated  and  distributed)  anplied  to 
intercepted  nart  of  shell.  In  composing  the  integral  equations,  we  shall  select,  as 
main  unknoim  functions 


4s 

•■■’W-f . 

as 


(6.3) 


This  makes  it  oosslble  in  subsequent  calculations  to  avoid  differentiation  of 
initial  geometric  and  elastic  parameters,  which  essentially  lower  the  accuracy  of 


the  calculation. 


Further,  one  should  consider  the  equality 


S 
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(6.9) 


By  introducing  the  values  (6.8)  and  (6.9)  into  equations  (6.3)  and  (6./*)  and  by 

integrating  both  sides  of  equality  from  4.  to  s,  we  will  obtain  a  system  of  normal 

+  )V„0(»  +  ,,  (a)  /„  4-  V')  (a)  /„  + 
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+  •(<*) /«  +  »<•»  (a) /.«  +  /., 

Eh\ 


integral  equations: 


where 


(6.10) 
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The  Txinctii  /</•  entering  into  equation  (6.10)  >  ''fill  be  determined  by  the 
dependencies 
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In  matrix  form,  aquation  (6.10)  has  the  form 


(6.11) 


The  aquation  contains  four  initial  parameters. 

C|-I(a);  C,-.ii<«»(a);  C,  =  d(a);  C,  =  d(')(o). 

Wa  now  consider,  as  an  axampla,  a  conical  shell  (Fig.  24). 
shall  have 
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In  this  case,  we 
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If  thickness  of  the  conical  shall  and  oaramuters  of  elasticity  are  constant, 
than  aquation  (6.11)  is  considerably  simplified: 
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I  f  f.  Fip.  24.  Conical  Shell. 
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The  solution  of  equation  (6.11)  is  made  by  the  method  of  subsequent 
anoroximations  according  to  scheme  indicated  in  Chanters  3  and  4. 

In  a  number  of  practical  problems,  it  was  found  sufficient  to  use  three-four 
anoroximations.  Another  method  of  annroximation,  which  may  be  used  for  solving 
equations  (6.11)  —  especially  with  a  gradual  convergences  of  process  of  successive 
anoroximations,  --  method  of  linear  anproximation. 

For  the  possibility  of  applying  this  method,  the  normal  operators  must  be 
presented  in  the  form  » 

For  operators,  entering  into  equation  (6.11),  this  is  readily  attained  by 
means  of  integrating  by  parts. 

Thus,  for  example,  for  conical  shell  we  will  obtain 

H  (s)  J  V )  (t)  <f;  _  j  H  (t)  ,J(I  )  (J)  ti-.  ^ 
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^  j f <9  <«• 


+  O  0)  I*  (0  »“ » (0  j  ’  (5)  </«. 

S 

- ;  a(s)  =  f^<yrf;. 

t£(5)A(C)  '  J  t 


•0(s), 

m 

«(.)-! 


Th«  iMthod  of  linear  approxlinatinn  is  applied  in  the  form,  as  discussed  in 
Chap.  3t  Sec.  6. 
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